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SUMMARY 


A two-day workshop was held at the NASA Lewis Research Center on August 5 and 6, 1992. The intent of the 
workshop was to inform the technical community of the seal code development activity sponsored by NASA 
Lewis and to provide a forum for participants to exchange information on respective activities in seals. We 
thank the workshop presenters, the participants and the peer committee for their contributions. 

Three codes were disseminated to the technical community for beta testing. They are the ICYL (cylindrical seal 
in incompressible fluids), GCYL (cylindrical seal in compressible fluids) and SPIRALG (spiral grooved gas 
seal) code. Close to thirty companies requested these codes at the workshop. 

Several needs, shortcomings, and problem areas have been delineated from the workshop on the NASA Lewis 
seals contract with Mechanical Technology Incorporated. The deficiencies include: validation data sets, 
convection flow contributions (inlet, inertia, turbulence, body forces), surface distortions, dynamics matrices in 
the CFD code, computational efficiency (speed, storage, hardware), and awareness of current seal user 
requirements and projected needs. 

The positive aspects, however, must not be overlooked. The codes are being worked into a transportable user 
friendly format. They are available for members of the peer committee and beta-users to apply toward their 
particular area of expertise. The findings from beta-testing need to be reported to enhance program 
effectiveness. This will be an on-going activity throughout this contract. 


Workshop Chairs: 

Anita D. Liang 
Robert C. Hendricks 
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PROGRAM SCOPE 

• SEVEN YEAR EFFORT (NASA Contract NAS3-25644) 

• PARALLEL PATHS - Scientific Code and Industrial Codes 

• DEVELOPMENT OF A KNOWLEDGE-BASED SYSTEM 

• TECHNOLOGY TRANSFER VIA WORKSHOPS 

• CODE VALIDATION VIA PUBLISHED DATA, IN-HOUSE TEST WORK 
AND COOPERATIVE PROGRAMS 
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PROGRAM ORGANIZATION 



Task 1 Cylindrical Seals Code 

Task II Cylindrical Code Augmentation (Labyrinth, Honeycomb, Damper, Brush) 

Task III Face, Wave, Grooved Seals 

Task IV Tip, Contact, Non-Continuous Seals 
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MAJOR TECHNICAL ACTIVITIES 

1 . DEVELOPMENT OF ANALYTICAL CODES 

2 . DEVELOPMENT OF EXPERT SYSTEMS 

3 . INTEGRATION OF CODES AND EXPERT SYSTEMS 
INTO ONE SINGLE FRAMEWORK 


Ad 

AEROSPACE TECHNOLOGY DIRECTORATE 
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CODE DEVELOPMENT 


INDUSTRIAL CODES - TWO-DIMENSIONAL CODES FOR DIFFERENT 

TYPES OF SEALS 


SCIENTIFIC CODE - A THREE-DIMENSIONAL CFD CODE FOR 

SEAL ANALYSIS 


SEAL TYPES - Cylindrical, bushing and ring, face, labyrinth, 

tip, damping, brush, electro-fluids, and 
"smart" 
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•Common with Industrial Knowledge-Based System 
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- PORTABILITY OF THE KBS 


NASA 
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Data Base 
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NASA 

Lewis Research Center 
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JgD SPACE PROPULSION TECHNOLOGY DIVISION N^S/V. 

AEROSPACE TECHNOLOGY DIRECTORATE ^ — — 

COMPUTER SYSTEM CONSIDERATIONS 

- A DESK-TOP SYSTEM (FOR ALL USER TYPES) 

- SUFFICIENT MEMORY FOR CODE EXECUTION AND FOR 
MULTITASKING 

- I/O COMMUNICATIONS WITH THE SCIENTIFIC CODE 

- COMPUTING POWER FOR GUI 

- PORTABLE SOFTWARE DEVELOPMENT TOOLS 

- PRICING 

igO SPACE PROPULSION TECHNOLOGY DIVISION 

AEROSPACE TECHNOLOGY DIRECTORATE 

COMPUTER SYSTEMS 

PERSONAL COMPUTERS WITH OS/2 OPERATING SYSTEM 

- Intel 30386 or Higher 
Math Coprocessor 

8 Megabytes of RAM 

- 80 Megabyte hard drive 

UNIX WORKSTATIONS 

- 16 Megabytes of RAM 

MACINTOSH/IBM INTEGRATED SYSTEM 

'i 
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TECHNOLOGY DIVISION 


NASA 

Lewis Research Center 


PORTABILITY OF SOFTWARE 


Fortran and C for Codes 

GUI - User Interface and C ++ compatible with OS/2 Presentation 
Manager and OSF/MOTIF 


Expert Systems - NEXPERT or CLIPS 


An SPACE PROPULSION TECHNOLOGY DIVISION NASA 

aerospace technology DIRECTOR te Lewis Research Center 


STATUS 

• SEVERAL INDUSTRIAL CODES HAVE BEEN PRODUCED AND 
READY FOR DISSEMINATION 


• A 3D CFD CODE FOR CYLINDRICAL SEALS HAS BEEN 
COMPLETED 


• THE KNOWLEDGE-BASED SYSTEM FRAMEWORK HAS BEEN 
DEFINED 


• RESULTS AND STATUS OF THE PROGRAM HAVE BEEN 
REPORTED IN SEVERAL MAJOR CONFERENCES 


• SEVERAL COOPERATIVE PROGRAMS HAVE BEEN ESTABLISHED 


• A TECHNOLOGY TRANSFER PLAN HAS BEEN FORMULATED 
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SPACE PROPULSION TECHNOLOGY DIVISION NASA 

AEfiosPACE technology oirectomtz Lewis Research Center 

"" CODE DELIVERABLES 


COPE QR MODULE APPROXIMATE DELIVERY DATE 


1. 

CFD Cylindrical Code 

02/01/92 

(09/92) 

2. 

Augmented CFD Cylindrical Module 

02/01/94 

3. 

CFD Code, Face, Wave, Non-Continuous Module 

09/01/95 


4. 

CFD, Tip, Contact, Non-Continuous Module 

09/01/96 


5. 

GJOURN 

02/01/91 1 

| 

6. 

ICYL 

03/01/91 

\ (07/92) 

7. 

SPIRALG 

04/01/91 J 

8. 

IFACE 

02/01/92 i 


9. 

GFACE 

03 / 01/92 

V ( 10 / 92 ) 

10. 

SPIRALI 

04/01/92 J 

11. 

FACEDY 

02/01/93 


12. 

RINGDY 

03/01/93 


13. 

LABYRINTH 

04/01/93 


14. 

FACECON 

02/01/94 


15 . 

DISTORTION 

04/01/94 


16 . 

Additional Codes 

04/01/95 


17 . 

Industrial KBS 

04/01/96 


18 . 

Scientific KBS 

04/01/96 



Ad 

AEROSPACE TECHNOLOGY DIRECTORATE 


SPACE PROPULSION TECHNOLOGY DIVISION 


NASA 

Lewis Research Center 


NEAR TERM ACTIVITIES 
CONTINUE CODE DEVELOPMENT 

ADD TURBULENCE AND INERTIA IN COMPRESSIBLE CODES 
DEVELOP USER INTERFACE FOR THE CFD CODE 
DEVELOP GRAHICAL USER INTERFACE 
DEVELOP EXPERT SYSTEMS 
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92P48 


Industrial 
Code Development 


Presented by 

W. Shapiro and A. Artiles 
Mechanical Technology Incorporated 
Latham, New York 

Presented to 

NASA Seals Workshop 
Contract NAS3-25644 


□ Compile and generate sets of verified 2D and 
simplified 3D or 1 D codes 

□ Codes are intended for expeditious parametric 
studies, analysis, and design of a wide variety 
of seals 

□ Integration is accomplished by the industrial 
KBS; additional functions of the KBS are: 

— User-friendly interaction 
— Contact sensitive and hypertext help 
— Design guidance 
— Expandable data base 


5 August 1992 


OBJECTIVES 
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CODE DELIVERABLES 



Approximate 

Code Module 

Delivery Date 

GCYL (Gas Cylindrical) 

02/01/91 



ICYL (Incompressible Cylindrical) 

03/01/91 



SPIRALG (Gas Spiral Groove) 

04/01/91 

^ 

CFD Cylindrical Code 

07/31/92 



FACE (Incompressible Face) 

12/31/92 



GFACE (Gas Face) 

12/31/92 


SPIRALI (Incompressible Spiral Groove) 

12/31/92 



FACEDY (Face Dynamics) 

09/30/93 



RINGDY (Ring Dynamics) 

09/30/93 


LABYRINTH (Gas) 

09/30/93 



Augmented CFD Cylindrical Module 

02/01/94 


FACECON (Face Contact) 

02/01/94 


DISTORTION (Thermoelastic Distortion) 

04/01/94 


Additional Codes: 

04/01/95 


Brush 



Damping Seal 



CFD Code, Face, Wave, Groove Module 

09/01/95 


Industrial KBS 

04/01/96 


Scientific KBS 

04/01/96 



CFD, Tip, Contact, Noncontinuous Module 

09/01/96 
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Computer Code GCYL 
Gas Cylindrical Seals 


GCYL CAPABILITIES 


□ Varying geometries 

□ Variable or constant grid (30 x 74) 

□ Shaft eccentricity and misalignment 

□ Specified boundary pressures and periodic boundary 
conditions 

□ Symmetry in axial direction 

□ Determining load (function of displacement) or seal 
position to satisfy given load 

□ Choice of English or SI units 

□ Frequency-dependent stiffness and damping coefficients 

□ Inlet and outlet inertia* 

□ Turbulence* 


To be added. 


CYLINDRICAL GEOMETRIES 



Circumferential 

Multilobe 

(With or without grooves) 



Circumferential 
Rayleigh Step 



Direction 



Rayleigh Step in 
Direction of Flow 




Self-Energized, 

Hydrostatic 

(Inherent compensation, 
orifice compensation, 
spot orifices, recesses) 


SS2I65 


SEGMENTED RING SEAL 


Item 

1 

2 

3 

4 

5 

6 

7 

8 



Description 

Segmented Ring 

Rayleigh Step 

Spring-Radial 

Spring- Axial 

Housing 

Cover 

Stop Pin 

Seal 

Sleeve 


Material 

Carbon 

Inconel X-750 
Inconel X-750 
Stainless Steel 17-4 PH 
Stainless Steel 17-4 PH 
Stainless Steel 17-4 PH 
Teflon 
Inconel 718 
Hard Chromium Plated 



832)66 


If) 


_ 


GCYL OUTPUT 


□ Clearance distribution 

□ Pressure distribution 

□ Leakage at specified flow paths 

□ Load and load angle 

□ Righting moments 

□ Viscous dissipation 

□ Cross-coupled, frequency-dependent 
stiffness and damping coefficients 

□ Plotting routines (pressure and 
clearance distribution) 


COMPRESSIBLE REYNOLDS’ EQUATION 


_a_ 

ae 




;A im 

ae 


, a (ph) 

5T 


where: 

2 = z/R, H = h/C 0 , T = t/t 0 , P = p/p 0 , 

. 6luoR 2 . 12uR 2 

= p 0 C 2 ’ ° = p 0 C 2 
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INTEGRAL FORMAT 


/ V * QdA - Q • n dS = ~ J (PH) dA 


where: 

Q-H 3 p|£ + A (PH) 

o9 


ASSUMPTIONS 


□ Laminar flow 


p^-U/h 
r dy dy 


where: 


t = Shear stress 
H = Viscosity 
V = Fluid velocity 
U = Slider velocity 
y = Film height 


□ Inertia is small and neglected compared to viscous shear 


□ Pressure across film is constant 


□ Height of film is small compared to other geometric 
dimensions; curvature is ignored 

□ Viscosity is constant 


□ Gas is iosothermal 
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TURBULENCE 


□ Re = £^>1000 

H 

Where: 

p = Density 
U = Surface velocity 
h a Film thickness 
Re = Reynolds number 

□ The most common approach to turbulence is to use G factors 
that modify the viscosity. G factors are dependent on the 
Reynolds number and pressure gradients 

— Couette Reynolds number: Re = h Rcop/p 

— Poiseuille Reynolds number: Re* = h 3 | VP | p/p 2 

— G x = Min [G x (Re), G p (Re*)] 

— G z = Min [G 2 (Re), Gp (Re*)] 


UNWRAPED SEAL SURFACE 







FLOW BALANCE CELL 




FLOW BALANCE EQUATION 


Q„ ^5 + Q- AZ| -’ . A* ^Uo- ^fcl o + AZ| o- aZm 

Ui2 — + u 12 ~2~ + u 14 + U 14 — g- - U 34 -g- - Q 34 — g— 

“ ^23 ■ 0>3 ^2^ * (^Qj + A9j.l)(AZj + AZm) ^ 

+ ^source 


Where, for example 




[ 12 


-H? G x (P, l|t1 ■ P,,|) 
A6| 


+ A P-J2 Hi 


and 


P 12 


gy ±fy*i 

2 


NEWTON-RAPHSON LINEARIZATION 


Fi( (Pi.P 2 ,P 3 ,...Pk,...P 9 ) -o 


□ The partial derivatives for the Newton-Raphson linearization 
may be calculated numerically: 

{F|, (Pi,P 2 ....Pk +£/2....P 9 ) 

- F if (P,,P 2 ,...P k - e/2, ...P 9 )}/e 

□ Then the linear equation for determining iterated pressures is: 

Fi, (old) + £ |& (old) (P k (new) - P k (old)) - 0 


□ The equation is put into the matrix format where Pj is the 
current pressure being solved for 


9Fh 

<?Pk 


[Cj] {P|} + [EJ {Pj.-,} + [Dj] {P j+ ,} - {R,> 


The system of equations is solved by the column matrix 
method. 
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FREQUENCY-DEPENDENT SPRING AND 
DAMPING COEFFICIENTS 


/ V • QdA - <^ Q • n ds = - ^ J (1 + P) HdA 


□ Introduce the RHS and perturb clearance and pressure 
H-e* 01 , P' mpe** 


then, perturb the resulting equation with respect to 
eccentricity to obtain 

dP 

ctek 

Q The real part of the derivative with respect to eccentricity 
integrated over the area is stiffness. While the imaginary 
part divided by a and integrated over the area is damping 


FREQUENCY-DEPENDENT SPRING AND 
DAMPING COEFFICIENTS (continued) 


Jv-SdA. jQ.nds = .^r J(1 +P) HdA 

□ The final set of linear difference equations for the complex stiffness pressure 
derivatives {P ■ *] are obtained 

tC' 1 ] {Pj k } + [i']{P;.f} * tDi]{P> t } - {Ri*} - tC» k] { £ j} 

- [Di k ] {Pj +1 } 

where: 

[C*i] - [Ci] + io 
{R , k } - {& *<} - io 

cj, - HC, A, + HCjAjj + HC3A3 + HC4A4 

[ci.«] - **: W - 

n 

Q The system of equations is solved by the column method in a directly 
analogous manner to that used in solving the steady-state equation. 

The principal difference is that all the matrix operations are performed 
using complex arithmetic. 
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STIFFNESS AND DAMPING COEFFICIENTS 


Speed (rpm) 


48,000 

48,000 

Excitation Frequency (rpm) 


48,000 

0 

Stiffness Coefficients 

Principal X (Ib/in.) 

KXX = 

0.965 

5885 

Cross coupled (Ib/in.) 

KXY = 

1942 

7267 

Cross coupled (lb/rad) 

KXA = 

0.730 

-0.239 

Cross coupled (lb/rad) 

KXB = 

1.291 

1.298 

Cross coupled (Ib/in.) 

KYX = 

1040 

-7116 

Principal Y (Ib/in.) 

KYY = 

17,670 

13,050 

Cross coupled (lb/rad) 

KYA = 

1.450 

3.258 

Cross coupled (lb/rad) 

KYB = 

0.448 

1.877 

Cross coupled (in. -Ib/in.) 

KAX = 

0 

0 

Cross coupled (in. -Ib/in.) 

KAY = 

0 

0 

Principal A (in. -lb/rad) 

KAA = 

639 

421 

Cross coupled (in.-lb/rad) 

KAB = 

72 

192 

Cross coupled (in. -Ib/in.) 

KBX = 

0 

0 

Cross coupled (in. -Ib/in.) 

KBY = 

0 

0 

Cross coupled (in.-lb/rad) 

KBA = 

-194 

-294 

Principal B (in.-lb/rad) 

KBB = 

221 

134 

Damping Coefficients 

Principal X (Ib-sec/in.) 

DXX = 

1.658 

1.406 

Cross coupled (Ib-sec/in.) 

DXY = 

-0.7059 

-1.859 

Cross coupled (lb-sec/rad) 

DXA = 

0 

0 

Cross coupled (lb-sec/rad) 

DXB = 

0 

0 

Cross coupled (Ib-sec/in.) 

DYX = 

0.918 

3.012 

Principal Y (Ib-sec/in.) 

DYY = 

1.521 

1.897 

Cross coupled (lb-sec/rad) 

DYA = 

0 

0 

Cross coupled (lb-sec/rad) 

DYB = 

0 

0 

Cross coupled (in. -Ib-sec/in.) 

DAX = 

0 

0 

Cross coupled (in. -Ib-sec/in.) 

DAY = 

0 

0 

Principal A (in. -lb-sec/rad) 

DAA = 

0.090 

0.1130 

Cross coupled (in. -lb-sec/rad) 

DAB = 

-0.031 

-0.047 

Cross coupled (in. -Ib-sec/in.) 

DBX = 

0 

0 

Cross coupled (in. -Ib-sec/in.) 

DBY = 

0 

0 

Cross coupled (in. -lb-sec/rad) 

DBA = 

0.026 

0.039 

Principal B (in.-lb-sec/rad) 

DBB = 

0.067 

0.069 

360° cylindrical seal; L 
p = 3 x 10‘ 9 Ib-sec/in. 2 ; 

= 1 in.; D = 1 in.; 

C = 0.001 in.; 


0 gage pressure at both ends 
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INLET INERTIA 


where: 

C D m Coefficient of discharge 
A 0 « Orifice area, in. 2 
P s - Upstream pressure, psla 
P R m Downstream pressure, psia 
7 - Ratio of specific heats 

G c - Gas constant, sec 2 /°R 
9 - Absolute gas temperature 

Q * Orifice flow, Ib/sec 


lf(^).P CR .then(^) = P 0R 
where: 



INTERNAL HALVING ROUTINE 



24 





SPIRAL-GROOVE FACE SEAL 
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CAPABILITIES 


□ Shaft seals and face seals 

□ Compressible flow 

□ Finite eccentricity and misalignment 

□ Four degrees of freedom for shaft seals 
(three for face seals) 

□ Frequency-dependent dynamic coefficients 

□ Arbitrary end pressures 

□ Predicts load, flow, power loss, dynamic 
coefficients, shaft displacements, and 
minimum film thickness 


ASSUMPTIONS 


□ Laminar isothermal flow 

□ No fluid inertia 

□ Ideal gas law 

□ Thin film approximations are valid 

□ Narrow groove theory 

□ Linearized time dependence 

□ Ideal surfaces 
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TECHNICAL FEATURES 


□ Finite misalignment 

□ Genera! algorithm for treating frequency-dependent 
dynamic coefficients 

□ Either forces or displacements may be specified 
in ail degrees of freedom 

Q Optional implementation of automatic numerical 
damping algorithm 

□ Optional implementation of Romberg extrapolation 
algorithm 


COORDINATE REFERENCE FRAMES 



ECCENTRIC SPIRAL-GROOVED SHAFT SEAL 



Eccentricity Ratio, e x 


SPIRALG THEORY 


GLOBAL AND LOCAL PRESSURES 






GLOBAL AND LOCAL PRESSURE GRADIENTS 


Q When the number of grooves becomes large, the sawtooth 
portion of the local pressure variation may be approximated 
with linear representations 

A0 A0 A0 ^ 30 X A0 V ae j A0 


A0 g 

A0 


a. and 


A 0, 

A0 


ae 




□ The corresponding relationship in the traverse direction, 
is obtained in a similar manner. 


□ The remaining two equations required to solve for the four 
local pressure derivatives are obtained from continuity 
considerations. 


CONTINUITY CONDITIONS 


□ Continuity of pressure at groove-ridge interface 



a Continuity of flow at groove-ridge interlace 
+ ^ r ((02 - coi) sin (3 

- 4 [^ (£)•“•<£)] 

+ j r(o3i - to,) sin 0 


□ Four equations in four unknown local derivatives that can be 
solved for in terms of global derivatives 
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Q, --(1 +P)[m? (k 2 A A^ P + Ask4 R Sin p. A ( a 5 + H f )Rj 
Q, - -(1 + P) [h? (k, ± + ^ . Ajtc, R cos p] 

where: 


p = P-Po 

Po 


Q.^ Ro q 
C 3 p 0 

7 (0 t 

t * -r— t 

2A 

u h ' 
Wf ■£ 

r 

c S 

s “r; 


. 6 poo R 3 

Po G 2 

8 = !k_L!k 
c 

A 6 = A8ooa (1 - a) sin p 

a ~r 3 r 

lr + lg 


r ta 2 - a)! 

CO SS “~ 

Oi 
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k VECTOR 

□ The column matrix containing spiral-groove coefficients, 
k, (a. P, H, is: 



a (1 -a)(r 3 - 1) 2 sin 2 p + r 3 
(1 - a) r 3 + a 

a (1 - a) (r 3 - 1) 2 sin ft cos ft 
(1 - a) T 3 + a 

a (1 - a) (r 3 - 1) 2 cos 2 p + r 3 
(1 - a) T 3 + a 

(r 3 -i) 

(1 - a) r 3 + a 
(1 - a) T + a 

r 


(f- 1) sin p 
(1 - a) r 3 + a 


a(1 - a)(r 3 - 1) sin p cos p 
(1 - a) r 3 + a 


a (1 -a)(r 3 -1)(r- 1) cos p + af + (1 -a) r 3 
(1 - a) f 3 + a 



□ The column matrix containing spiral-groove coefficients, 
k, (a, p, r), is: 


a (1 - a) (r 3 - 1) 2 sin 2 p + r 3 
(1 - a) T 3 + a 

a(l -opCr 3 - 1) 2 sin p cos p 
(1 - a) T 3 + a 

a (1 -a)(r 3 -I) 2 cos 2 p + r 3 
(1 - a) T 3 + a 

(r 3 -i) 

(1 - a) f 3 + a 
(1 - a) r + a 

r 

(T ■ 1) sin p 
(1 - a) r 3 +a 

a (1 ’ a) (f 3 - 1) sin p cos P 
(t - a) T 3 + a 

a (1 -a)(r 3 - l)(f - 1) cosp + ar+(t -a) r 3 
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CELL CONTROL VOLUME 



NEWTON-RAPHSON ITERATION 

Plow Balance 

Q 12 aS ' 12 + Q^ASlg + Q * 4 AS | 4 Q * 14 AS ' 14 - O 34 AS J 4 - 

/ Q 34 A$34 . “ ^23 ^^23 * ^23 ^23 * 

* ft &M ' ^ ^ + 1/2 ■**' 

+ (otS + Hf)j + j . 1/2 ^A 4 + (oc 8 + H r ), . 1/2 j . - 1/2 AA 3 
+ (®I + Hr)*. ^2 i + ^ M 2 

• ,••• . __ '■. J-'-. •: " '*' * 

For Steady State, RHS s 0 

Fii (Hr, P 1 ^%. ^, P4 , P 5 ) If*/#; Pq , p 9 ) - 0 

Apply Newton-Raphsori 

9 ^ p 

*v* s ^ .. : 

where: 

3Pii fj| (H,, P, , ... ■ ■ P k ♦ ti P 9 ) - Fn (H r , P, , , , „ P 9 ) 

3 Pk • n 
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NEWTON-RAPHSON ITERATION 

Flow Balance 

Ql2 AS12 + Qi2^S 12 + Qi 4 AS 14 Q 14 AS u 'Q^AS^j * 

■ ^34 AS 34 ‘ ^23 AS 23 ’ ^23 AS 23 * 

* ■ jo + ^'i) [( a ® + ^r)l + 1/2, j + 1/2 AAl 
+ (aS + H r )j + 1/2 j . 1/2 AA4 + (cx 5 + H r )j . 1/2t j . 1/2 AA3 
+ (cx6 + H r )i . 1/2, j 4- 1/2 aaJ} 

For Steady State, RHS = 0 

Fjj (H r , Pi, P2.P3.P4,P 5 ,P6,P7,P8 f P9) -0 

Apply Newton-Raphson 

where: 

a Fit Fy (H r , Pi P k + ri, . . P 9 ) ■ F, (H f , P t , . ■ P,) 

a Pk r[ 

COLUMN MATRIX FORMAT 

□ Generate system of equations in the following format: 

[cq {Pj new } + [eo {p,t} + [D']{ p,:r} - m 

where: 

d J 5 l 

C| ’ itk “ 9Pi+k, | 

d __?!l 

u+k 9P|+k, 1-1 

ni . 

d, ’ ,+r afWi 

k = -1, 0, 1 ; i-2 M-1 

□ The interior elements of the column vector {R 1 } are 

1 (C|,i+k P|+k,J + E | i( . k Pj+k, j-1 + 0| kk Pj+k. j+i ) - F|| 

□ Equations solved by column or transfer matrix method 
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FLOW DIAGRAM FOR LOGIC USED IN 
SUBROUTINE SPIRAL 



Mechanical Technology Incorporated 92?73 <m) 

92P48 


34 















Computer Code GFACE 
Gas Face Seals 


GFACE CAPABILITIES 


□ Varying geometries 

□ Variable or constant grid (30 x 74) 

□ Shaft eccentricity and misalignment 

□ Specified and periodic boundary conditions 

□ Determining performance as a function of 
position or position to satisfy a given load 

□ Frequency-dependent stiffness and 
damping coefficients 

□ English or SI units 
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FACE SEAL CONFIGURATIONS 



+ 


Circumferential 
Tapered Land 



Radial Rayleigh Step Radial Tapered Land 


86265 • 1 
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POLAR GRID MESH SYSTEM 



COORDINATE SYSTEM FOR GFACE 
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GFACE THEORY 


REYNOLDS’ EQUATION 



where: 

R = r/r 0 , H = h/C 0 , T = t/t 0 , P = p/p,, 

A = 6 ^ r ° t _ 12 ^ r o 
" PoCg °" p 0 C* 


FLOW-BALANCE CELL AND ASSOCIATED 
GRID NETWORK 


9(PH) 

at 





□ The net flow through a cell can be expressed as: 



ar. 

ARj.-j 

_ A9i 

+ Q,4 

A0i 1 

Ql2 

+ ®i2 

2 

’ + <^4 ^ 

2 

-Q34 

^i.Q- 

2 U 34 

AR m 

2 

1 

O 

' ^23 

A0j-i 

2 


□ Q| 2 means the mass flow per unit length across the plus side 
of cell boundary 1 -2, etc, 

□ The Q’s are dimensionless mass flow per unit length, except 
for Q |n , which is a dimensionless source inlet flow. 

□ In the 0 direction: 


PH 3 3P AR 

r ae 2 


fAR PH 


AH 

2 


□ In the R direction: 

Q= PH 3 R M 

3R 2 
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CELL PRESSURES AND DERIVATIVES 


□ Q is defined as: 


12 ,uG c T a q 

n 2 r 3 
Po °0 


□ Pressures are taken as the average pressure across 
the boundary. For example: 

p j + Pi, 1+1 

n2 „ — l__ — 1_ 

and 


3P _ Pj, j+i ~ Pj, j 

39 12 A9j 


EXTERNAL MASS FLOW 


Q in = OFC X A 0 x P s 


r , teT * 1 " 

PsJ L lP S J J 


1/2 


where: 


OFC • 


12 pC a 

PoCo 


•\f 2 TGcTa 

V v - 1 


A 0 =itd 0 H 5 for inherent compensation 


rt d? , 

A 0 = ----- for orifice compensation 
4 Cg 
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FLOW CHOKING 


If 



^ PcR. 



PCR 



where: 



Also, if 


fa 

p s 



i 

Pr^Ps 


, andPg 


P n 


R • 


then this condition implies backflow through the orifice. 


NEWTON-RAPHSON ITERATION ON 
FLOW BALANCE EQUATION 


f (i 7> + 


5 

z 

k-1 



^ p^new) . 


0 


where the partial derivatives are explicitly determined, e.g., 

3fj, j f(Pi . P 2 1 Pk + ... P s ),. j • f(Pi ■ P 2 • — Pk + e/ 2, ... Pg),, j 

3P k e 
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COLUMN MATRIX FORMAT 


P The linearized Newton-Raphson equations may be 
written in the form: 

Cj Pj + Ej Pj.i + Qj Pj+i w Rj 


□ The column method is used to solve the new pressures 
in the set of m x n equations. 


VISCOUS POWER LOSS 



r + di/dz 



r 


b) Fluid Element 

9282 * 
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VISCOUS FRICTION TORQUE 



T, 

Po CR* 


JJH| Rd9dR .JJ 


? TT Rd6dB 


FLOW ACROSS CIRCUMFERENTIAL LINE 




FLOW ACROSS AXIAL LINE 



NONUNIFQRM CLEARANCE GEOMETRIES 



L 11 

, (— i 

* mm Taper: Type 2 Section C-C 

• Rayleigh Step 



Section P-b 
> Circumferential Taper 
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DOUBLE TAPERED-LAND CLEARANCE DISTRIBUTION 



DOUBLE TAPERED-LAND PRESSURE DISTRIBUTION 





STATUS OF INDUSTRIAL CODES 



ICYL 

GCL 

SPIRALG 

IFACE 

GFACE 

SPIRALI 

Formulation 

C 

c 

C 

C 

C 

C 

Coding 

c 

c 

c 

P 

P 

c 

Documentation 

c 

c 

c 



c 

KBS Integration 

c 

c 

c 





C = Complete 
P = In preparation 

ICYL = Incompressible Cylindrical Code 
GCYL = Gas Cylindrical Code 


SPIRALG = Gas Spiral-Groove Code 
IFACE = Incompressible Face Seal Code 
GFACE = Gas Face Seal Code 
SPIRALI = Incompressible Spiral-Groove Code 
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Preceding Page Blank 



INCOMPRESSIBLE FLOW CODES OVERVIEW 


W. Shapiro 

Mechanical Technology Incorporated 
Latham, New York 


• OVERVIEW 3 CODES: 

O Capabilities 
O Status 

• ICYL « Cylindrical Seals 

O Introduction 
O Assumptions 

O Formulation & Solution Method 

• SPIRALI « Spiral Groove Seals (Cylindrical and Face) 

O Introduction 
O Assumptions 

O Formulation & Solution Method 

• IFACE • **■ Face Seals 

O Examples 
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INCOMPRESSIBLE FLOW CODES CAPABILITIES 



ICYL 

IFACE 

SPIRALI 1 

Geometry 

Cylindrical 

Face 

Both 

Turbulent Flow 

Yes 

Yes 

Yes 


Yes 

Yes 

Yes 

■m 

No 

No 

Yes 

Film thickness 

Arbitrary 

Arbitrary 

Axis-symmetric 

Eccentricity 

Finite 

Finite 

Perturbation 

Pressure Pockets 

Yes 

Yes 

No 

Spiral Grooves 

No 

No 

Yes 


INCOMPRESSIBLE FLOW CODES «* STATUS 



ICYL 

IFACE 

SPIRALI 

Formulation 

Complete 

Complete 

Complete 

Coding 

Complete 

In Progress 

Complete 

Documentation 

Complete 


Complete 

KBS Interface 

Complete 










































Incompressible Cylindrical Program ■»* Capabilities 

ICYL by Dr. Antonio Artiles 


• 2-D cylindrical geometry. 

• Rotation of rotor and housing. 

• Roughness of rotor and housing. 

• Laminar or turbulent flow. 

• Inertia pressure drop at inlets to fluid film 

O from pressurized pockets 
O from seal ends 

• Arbitrary film thickness distribution 

O steps 

O pockets (pressurized or not) 

O tapers 
O preloaded arcs 

• Rotor position relative to housing described by four degrees of freedom 
O 2 translational 

O 2 rotational 

• User specifies: 

O Rotor lateral position -or- External forces 
O Rotor angular position -or- External moments 
O Pocket pressures -or- orifice size 

• Dynamic coefficients 
O 1 6 stiffness 

O 1 6 damping 
O critical mass 


Incompressible Cylindrical Program Assumptions 

• Isothermal and incompressible flow 

• Turbulence: bulk flow model with separate friction factors 

• Fluid inertia effects at film entrance use loss coefficients. 

• Fluid inertia effects in the film (additional to those inherent in turbulence) 
are negligible. 

• Pressurized pockets deep (constant pressure) 

• Surface roughness < < film thickness < < seal dimensions 

• Wall roughness is isotropic 
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CYLINDRICAL SEAL GEOMETRY SCHEMATIC 

(CONCENTRIC ALIGNED POSITION) 
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AXIAL CROSS SECTION OF SEAL 
WITH ECCENTRIC ROTOR 

(FILM THICKNESS EXAGGERATED) 



ICYL Program Governing Equations 


Film thickness 

H = H 0 - (e x <ZB) COS0 - (^-Zyljsine 


dH 

dt 


— +Z— )cos6 - 

dt dt) 


**-z^Le 

Si Si) 


momentum 


continuity 


(fj R ej + fb Re b) u= + ( Re /j U j + Re Jb U t) 

2 \iRdQ 2 

{fjtej + fb Re b) v _ 

2 \i dZ 


(UH) 

RdQ 


— (VH) + — = 0 
dZ dt 
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Reynolds numbers, friction factors 


Re- 


fr 


— , /tejSlOOO ( laminar ) 

■~(1 -3l[ 2 + 2s 3 ) + / , (3* 2 -2£ 3 ) 1 1 000 <3000 
Ke i 

fi , Re j 2 3000 (turbulent) 


Re r 1000 
2000 


/,* = 0.001375 1 



Vl 

, 10 6 3 

+ 4 ^ fj 


Reynolds numbers, friction factors 


Selected friction foctor (Moody) 
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ICYL Program w Boundary Conditions 
circumferential ends: 

P( Z, 0 S ) = 0 and P( Z, 0 e ) = 0, 
or periodic boundary conditions: 

P< Z, 6 s >= P( Z, 0 e ) and U( Z, 0 S ) = U( Z, 0 e ). 


left end: 

P(-L/2,0) = P, - K e y 2/ oV n 2 . 

right end: 

P(L/2,0) = P, - K' y 2 pv n 2 , 
or mid-length symmetry: 
V(O,0) = O 

pocket boundaries: 

P(Z ,0)= P p - K e y 2/ ov n 2 . 

orifice pressure drop: 

Pressurized Pockets 

p, - p, = *»«?,) |(^ C J 

orifice flow: 

a = ftf <is + f— 

i i, * 
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Discretization of continuity equation 



'Ilf ( 4 v 4 ZM X 4e '* 4e '-' ) '° 


Discretization of continuity equation at pressurized boundaries 



v « = 

Az* A z.- i AO 

+ g ( U 2 h 2 ~ U 3 h 3) + 

A 0,-i 

1 - V 2 A 2) 2 


dh tJ (Az,+ Az M ) A Oy ^ Az M A0 ; _, 


(A0^_i + A ZiJAy 


dx 4 


2 
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Evaluation of velocity components 


dp 

— ,u,v 
30 


fjRey^ Re . u + 1 2/, 2 - (Re/jA j + ReJ b \) = 0, 


^,«,v 1 - f AIA v + 12A 2 ^ = 0, 


3z 


3z 


^<«-2 M 4 *'". 




^fr= («-2A t ) 2 + v 2 , 



Evaluation of velocity components 


Eiti Ell u - 

J ^ I 


V + V 


A 0 , 


= 0 


PMJ PiJ u~ +U* 
^z i 9 2 


, v 


Pi*1J+1 Pl+'\J 


AQj 


,u 


V + V 



Pj*\r\ P'j'i «_ 

Az, 


+ u * 
2 



= 0 
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ICYL: Forces, Moments, Torque 




COS0 

z 7 . 

L 9, 

sin0 


= / / p 
-l e s 

-Z sin0 

K 

w 


Z cos0 


RdQdZ 


f = Te z = ffr x f dA 

A f 

= R ff*r X (*“«r) dA 

A f 

T - rr fA ^ _ W“- 2A /) -/^„(«- 2 A ft ) 

2C a Jj de 72 h 

Af 


} </0<*Z 


Solution of Rotor position and pocket pressures. 

/» = ~fjg 
fy( r ) = -4 

m ,( r ) = "«» 

Ps - P pi = sn 8(<lTt) Ki{<lnf> for pocket 1 , 

> 

P s ~ P P 2 = sng(q r2 ) A r2 (q r2 f, for pocket 2, etc. 
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SPIRALI (Spiral Groove Incompressible Flow) «** Introduction 


SPIRALI by Dr. Jed Walowit 


Spiral groove seals 


Provide: 


I stability 

□ increase in-phase 1 component of force 
□decrease out-of-phase component of force 
Bload support 

Isealing (pumping against pressure gradient) 


• Cylindrical seals symmetric grooves pump against each other 

• Face seals grooves pump against dam region 


'with the displacement 


SPIRALI Basic Assumptions 

• Isothermal and incompressible flow 

• Turbulence: Hirs bulk flow model generalized for separate friction factors 

• Film discontinuities use loss coefficients. Inertia effects (additional to 
those inherent in turbulence) throughout film. 

• Circumferential and transient effects use small perturbations to a steady 
state first order solution for a concentric, aligned seal. 

• surface roughness < < film thickness < < seal dimensions 

• Narrow groove theory 

O neglects edge effects and local inertia effects 2 
O valid for large number of grooves (> >2/7sin/?) 

• Machined surfaces are axisymmetric (except spiral grooves) 


2 due to groove to groove pressure variations 
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SPIRALI ■** Coordinate system 



SPIRALI *** Quadratic film variation 
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SPIRALI » INTEGRATED EQUATIONS 


MOMENTUM (0 AND s DIRECTIONS): 


fdu + v j^ + Hiy + 

V at 3s r 30 r t 


h ap 
r ae 


+ (?b" ? a) *T 


/ 




UcH 

r ae ’ r J 





CONTINUITY: 


H (rvh) + 7^ (uh) + ? =0 

r <3s r (30 dl 


SPIRALI Shear stresses, shear factors, and Reynolds numbers 
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SPIRALI 


Shear stresses, shear factors, and Reynolds numbers 


i - , ,( 2 a p!“ ~“«l ' 
'■ ■ 2 pl “'“ J ^ ^ 


(W ""a) = 


“ f a (R a ) “a) * 

4 h 


f.(R a ) = n o R » mo > f b ( R b) = n 0 R b m ° 


R a = 2h | u - u a | p/^ , R b = 2h |u - u b | p/(i 



Laminar 

Turbulent 

n o 

24.0 

0.0751 

m 0 

-1 

-0.25 


SPIRALI «r INTEGRATED DIMENSIONLESS EQUATIONS 


MOMENTUM: 


-12 = pM>(G,v,h) + p'R*f^ + v^ + ^ + l,^l , 

r ae K v ' H [$ as r ae ' r J 


-J| = p'4'(u,v,h) + p'R* 


+ Of | 

at as f ae ? 1 


CONTINUITY: 


iiiiyh) + 1 d(uh) + ah = 0 
r as f ae at 
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SPIRALI INTEGRATED DIMENSIONLESS EQUATIONS 


where 


and 


<K(u,v,h) 


(G - r<o)R a f a (R a ) + uR b f b (R b ) 


T(0,v,fi) 


R a < a (R a ) * RJb(Rb) f . 

h 2 


SPIRALI «*■ Boundary conditions 


u io , p in at seal inlet: 


U = U ln • 

P = Pin - ^P’R'O + 0 V 2 @ s = s ln 


p„ at exit: 


P = Pex @ S = S e 


S in depends on flow direction. S„ at the opposite side from S*. 
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SPIRALI Continuity conditions 


At S = Sj! 


(h + Ah)Vj = hv , 


Pj + 1 p v j 2 = P + ^P v2 ( 1 +0 


, Vj , h - Ah 


P. v, h 


-Ah 


s flow direction 


Jump in film thickness. 


C(R,h,v) 


1 - 


h - Ah 


, Ah < 0 (contraction) 


, Ah £ 0 (expansion) 


SPIRALI » Stator with inward pumping grooves 
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SPIRALI Spiral groove parameters, global and local pressures 



SPIRALI First order equations 
First order variables, H(S), P(S), 0(S) and V(S): 


^(O.V.H.I,) + R*V 


dU 

dS 


dP 

dS 


p*r(0,S/,fl,i f ) ♦ P'R'V^ 


Closed form solution: 


d(fHV) _ 0 
dS 


HnRln^in 

rH 
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SPIRALI «r Semi-Implicit Algortihm 


N differential equations: 


— = F,(S,Y 
dS ' 


,.y 2> . 


,V N ) , i = 1,2, ...,N 


Yj dependent variables 

{Y ncw } at the new position (S + AS) 

{Y} at the old position S 

{ y now ) . {YJ + AS([l]-^[k])- 1 {F(S + ^,Y 1 ,Y 2 ,...,Y N )) , 
= JjrF,(S +-y,Y 1t Y 2 ,Y, Y n ) , i,j = 1,2 N . 


SPIRALI ■*" Second order equations 

fi = H(S) + fi'(S,e,t) , u = 0(S) + u'(S,e,t) , 
v = V(S) + v'(S,0,t) , P = P(S) + p'(S,0,t) . 


i ap' 


= R* 

fp* 36 


^ + viu^ + + u ao_ + * Q/ + j v + fi/ , 

at as dS f ae 1 u v H 


.li&! = R* 
p- as 


( f*L + vi^ ♦ M., + ua^' 
at as dS ? ae ; 


+ *F* 0 u' + 4^v' + 'P* f| h / , 


vM^l + dVfh' + ?H-^ ♦ Mlv' + Q— + H— + r— =• 0 . 

as dS as dS ae ae at 
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CFD Research Corporation 

3325-DTriano Blvd. ■ Huntsville, AL 35805 ■(205)536-6576 ■ FAX: (205) 536-6590 


CFD^C 


DEVELOPMENT OF A CFD CODE FOR 
ANALYSIS OF FLUID DYNAMIC 
FORCES IN SEALS 


Presentation by 

A.J. Przekwas and M.M. Athavale 
CFD Research Corporation 
Huntsville, AL 


Presented at 

Second Seals Flow Code Development Workshop 
NASA Lewis Research Center 
Cleveland, OH 

August 5, 1992 


P-4MS 7/01 


OUTLINE 


• Objectives 

• Status Report 

• Code Capabilities 

• Rotordynamic Coefficient Methods 

• Test Results 

• Conclusions 


CFD'TC 


«»-«n5?/OJ 
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OBJECTIVE OF THE PROGRAM 


CFD^C 

• Develop Scientific 3D CFD Code for (CFDRC) 

- Prediction of Flow and Dynamics In Various Seals 

- Contribute to Data Base 

- Accuracy Standard for Industrial Design Codes 


Compile and Generate Set of Verified 
Industrial Codes 

(MTI) 

Seal Design Configuration with 

(MTI) 

- Knowledge Based System (KBS) 

- CAD/CAM and Visualization Tools 



OBJECTIVES 


• Develop Verified CFD Code for Analyzing Seals 

• Required Features Include: 

Applicability to a Wide Variety of Seal Configurations, 
such as: Cylindrical, Labyrinth, Face, Tip and Brush 
Seals 

Accuracy of Predicted Flow Fields and Dynamic Forces 
Efficiency (Economy) of Numerical Solutions 
Reliability (Verification) of Solutions 
Ease-of-Use of the Code (Documentation, Training) 
Integration with KBS 
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CFD CODE REQUIREMENTS 


CFD^C 


• Scientific Code Required Capabilities 

- 3D N-S Analysis in BFC Grids 

- Stationary and Rotating Frames 

- Incompressible and Compressible Flows 

- Variable Physical Properties 

- Steady-State and Transient Solutions 

- Rotordynamic Coefficients 

- Interfaces with Advanced Preprocessing 
(Grid-GUI) and Postprocessing Packages 


»Mllfr-7/OS 


SCIENTIFIC CFD CODE DEVELOPMENT 


CrD^C 


Task 1 : Develop a 3D CFD Code for Cylindrical Seals 

- for Annular, Tapered, Stepped 

- Verification of Code Accuracy 
Rotordynamic Coefficient Calculations 

Task 2: Augmentation of Code for Labyrinth and 

Damper Seals 

Tasks 3 & 4: Augmentation for Other Seal Configurations 


Note: Starting CFD Code = REFLEQS (developed by CFDRC 
under a contract from NASA MSFC/ED32) 


R41I4-7/0* 
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STATUS - 1991 WORKSHOP 


CFD1C 


AT 1991 WORKSHOP - 2D CAPABILITIES 

• Colocated Grid Formulation for BFC Grids 

• Strong Conservative Formulation of Momentum 
Equations with Cartesian Components 

• Choice of SIMPLE, Modified SIMPLEC 

• Higher Order Accurate Temporal Differencing 

• Higher Order (2nd, 3rd) Spatial Discretizations Available 

• Rotating Grid System for Stator-Rotor Configurations 

• Moving Grid Options for Arbitrary Rotor Whirling Analysis 


P<nt6 7«7 


CURRENT STATUS 


Accomplishments Since Last Workshop 


CFD^C 


• All Numerical Models Transferred to 3-D 


- colocated variables 

- higher-order schemes, etc. 

• Rotordynamic Coefficient Calculation Methods 

- circular whirl 

- moving grid (numerical scheme) 

• Seal Specific User Interface 

- grid generation 

- preprocessing 


•* 4 I lS I'Ot 
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CURRENT CODE CAPABILITIES 


CFD^C 


• Seals Code has: 

Finite Volume, Pressure-Based Integration Scheme 

- Colocated Variables with Strong Conservation 
Approach 

High-Order Spatial Differencing - up to Third-Order 

- Up to Second-Order Temporal Differencing 
Comprehensive Set of Boundary Conditions 

- Variety of Turbulence Models (k-e, Low Re k-e, 
multiple scale k-e) 

- Moving Grid Formulation for Arbitrary Rotor Whirl 

- Two Possible Ways to Calculate Rotordynamics: 

(i) Circular Whirl (ii) Shaker Method (moving grid) 


P41IV7XM 


CODE SYSTEM ENVIRONMENT 


CFD1C 



MO 
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FLOW SOLVER STRUCTURE 


CFD^C 



^•41 lf-7/1 1 


SEAL SPECIFIC CAPABILITIES 


CrD^C 


• GUI and Preprocessor - Geared for Seals Problems 

• Easy, Quick Geometry Definition and Grid Generation 

• Four Types of Cylindrical Seals: 

- Annular, Axial Step-Down, Axial Step-Up, 
and Tempered 

• Pull-Down Menus for Problem Parameter Specification 

• One Line Commands for 


- Automatic Grid generation 

- Integrated Quantities: Rotor Loads, Torque, etc. 

- Rotordynamic Coefficient 
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ROTORDYNAMIC COEFFICIENTS 


CFD3C 


• Relation Between Fluid Reaction Force and Rotor 
Motion 

- Small Perturbations 

- For Nominal, Centered, Rotor Position 


-Fy 

= 

K k' 

y 

+ 

C c 

y 

+ 

M 0 ' 

y 

-Fz 


-k K. 

z 


-c C- 

z 


. 0 M- 

z 


• K, k -- Direct and Cross-Coupled Stiffness Coefficients 

• C, c -- Direct and Cross-Coupled Damping Coefficients 

• M -- Lumped Mass (Direct Inertia) Coefficient 
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ROTORDYNAMIC COEFFICIENT - METHODS 


CFD^C 


-Fy 

= 

K k' 

y 

+ 

C c 

y ’ 

-f- 

M O' 

y ' 

. -Fz . 


.-k K. 

_ z 


- -c c. 

. z . 


-0 M- 

z 



Circular Whirl Orbit 
Method 


"Numerical Shaker" Method 
(with moving grid) 
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CIRCULAR WHIRL METHOD 


CFD^C 

• Rotating Frame -» Quasi-Steady Solutions 

• Calculate Quasi-Steady Solutions 

- 4 or More Whirl Frequencies, Q\ 

Integrate Pressures for each Qj to Generate Rotor 
Loads, Fy and Fz 

• Curve Fit to Calculate Rotordynamic Coefficients: 

- Fy = + K + cQ. - M Q 2 

- Fz = - k + CQ 


P«4MS-7/1S 


NUMERICAL SHAKER METHOD 


CFD^C 


Two Approaches to Calculate 
Coefficients 

Integrated Approach 

Fy = Fys sin(£2 1) + F YC cos(Q t) 
F £ — Fj£ sin(£2 1) + F cos(£2 t) 


Successive Time Step Approach 


Fyi ' 


y i 

0 

y i 

0 

Tzi 
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-y i 
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-y'i 

FyX 


y 2 
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y 2 
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Fzi - 


0 

-y 2 

0 

-y-i . 


K' = 

K + M ' & 



Horizontal Displacement Definition 

Y = A sin(Q t) + B cos(Q t) 


I 



Q. = Shaker Frequency 
(o = Rotor Speed 
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CODE ACCURACY CHECK & VALIDATION 


• Extensive Validation Effort 

- Checkout Problems 

- Benchmark Problems 

- Validation Problems 

- Field Problems 


P4IISM? 


CODE ACCURACY AND VALIDATION 


CrD^C 


• Partial List of Relevant Test Cases 

- Flows in Pipes, Channels, Very Narrow Annulii 

- 2D and 3D Driven Cavity Flows 
Laminar Flows: Wedge, Duct to 90° Bend 

- Rotating Flows: Disk in a Cavity, Stator-Rotor Config. 

Turbulent Flows in Annular and Laybrinth Seals 
(Texas A&M) 

Flow in Journal Bearings 
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DESCRIPTION OF THE FLOW PROBLEM 
• Simplified Turbine Cavity Problem: 


CFD^C 

outer wall 



NUMERICAL RESULTS 


Streamlines 


CrD=tC 





C Q = 1500 


C Q = 7200 


C Q = 13000 
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Computational results 

Cavity temperatures at three cooling flow rates 
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SAMPLE PROBLEM DEFINITION 


CFD^C 


• Annular Seal Geometry (After Nordmann, 1987) 


- RotorRadius 23.5 mm., Nominal Clearance 200(im 

- Seal Length 23.5 mm. 

Rotor Speed 1000-5000 rpm 


• Flow Conditions 


- Specified AP Across the Seal 
Specified Inlet Loss Coefficient 

Water Density 996 Kg/m^, viscosity 0.7 x 10"^ N-s/m 

• Computational Parameters 
k-e Turbulence Model 

- N x * N r * N e = 10 * 5 * 30 Computational Grid 

- 40 Time Steps per Cycle 
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FLUID REACTION FORCES 


CFD^C 


• Rotor Speed 3000 rpm., Q = 40 cps. 


Geometry 


Time Variation of 
Fluid Forces 


Z 


♦ 
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ROTORDYNAMIC COEFFICIENTS 


CFDHC 


Time Variation of Mass Parameter 

Stiffness Coefficients 
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STIFFNESS COEFFICIENTS 


CFD=K 

Qi rec t Cross-Coupled 




pa » »«PR/ta 


DAMPING COEFFICIENTS 
Direct 


CrD^C 

Cross-Coupled 
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TAPERED GAS SEAL (NELSON, 1985) 

D = 65.0 mm 
C e = 0.086 mm 

Cj = 0.172 mm 
3 L/D ratios: 0.1, 0.2, 0.4 


D 





Compressible, Turbulent, 
Transonic Flow 
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TAPERED GAS SEAL - RESULTS 
Moving Grid Method 


CrD^C 


L/D 

KN/m 

kN/m 

C N-s/m 

c N-s/m 

Exit Mach 
Number 

0.1 

948000| 

19700 

11.5 

0.016 

1.00 

1150000 

15429 

9.94 

0.043 

1.00 

0.2 

1700000 

75700 

43.9 

0.073 

0.96 

2125500 

60886 

38.62 

0.09 

0.97 

0.4 

2880000 

267000 

152 

0.27 

0.83 

3553200 

233820 

145.9 

0.57 

0.83 


| 1 Nelson, 1985 I I Present Results 
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CONCLUDING REMARKS 


CFMC 

• 3D CFD Code has been Developed for Annular Seals 

• State-of-the-Art Numerical Methods Employed 

- Colocated Grids, Pressure-Based 

- High Order Differencing 

- Moving Grids 

• Several Seal Specific Parameters 

- Rotor Loads, Torque, etc 
Rotordynamics 

• Initial Link with GUI Completed 

• Ongoing Validation Effort 


CONCLUDING REMARKS (CONTINUED) 
• Continue Code Testing & Validation 


CFD^C 


• Add Treatment for Inlet/Exit Regions 

- Geometry Definition 

- Grid Generation 


• Code Extensions 

- Augmentation of Code for Labyrinth Seals 

- Treatment of Damper Seals 

• Identified Needs 

- Multi-Domain Capability for Proper Treatment of 
Inlet/Outlet Plenums 

Two-Layer k-e Model for Near-Wall Treatment 


41 16 H?4 


84 


Seal Analysis Codes: KBS User Interface 

Presented by: 

Bharat B. Aggarwal 
Lynn Cowper 

Mechanical Technology Inc. 


Seal Analysis Codes: KBS User Interface 


Overview 


• KBS Executive 

• Access all seal analysis codes 

• Utility services like printing, browsing, plotting, etc. 

• System maintenance and configuration 

• Industrial Codes 

• Prepare input files for analysis codes 

• Run analysts code 

• Scientific Code 

• User input tailored for specific seal configurations 

• Four types of Cylindrical seals implemented. Others to be added in phases. 

• Prepare input files for preprocessor 

• Run the preprocessor 

• Communication capabilities to be added later 
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Seal Analysis Codes: KBS User Interface 


KBS System Implementation 


• Portability Requirements 

• User interest in Unix at the last workshop 

• Elected to support both OS/2 and Unix 

• Unix support for OSF/Motif environment only 


• Implementation Strategy 

• C++ and Commonview 3.0 class libraries from Glockenspiel 

• Portable to several Unix platforms with OSF/Motif interface 

• NASA will support IBM, HP, and SUN workstations 


Seal Analysis Codes: KBS User Interface 


KBS System Requirements 


• IBM PC Users 

• Intel 80386 (with 80387) or 80486 based machine 

• 8 Mb RAM 

• 80 Mb Hard Disc 

• OS/2 Version 2.0 

• Will continue supporting OS/2 version 1 .3 for another year 

• Workstation Users 

• IBM RS/ 6000, SUN, or HP 700 series workstations 

• 24 Mb RAM 

• 200 Mb Hard Disc 

• Unix with OSF/Motif user interface 
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Seal Analysis Codes: KBS User Interface 


System Architecture: KBS Executive 



Seal Analysis Codes: KBS User Interface 


System Architecture: Analysis Codes 


Industrial Codes 


Graphical U ser 
Interface 

"Z 


ASCII 
Input File 



Analysis 

Code 



ASCII 

Output File 


Scientific Code 


- 

Graphical User 
Interface 



Preprocessor 
Input File 


Pre processor 


' CFP Code 
Input File 

3 : 

CFP Code 
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Seal Analysis Codes: KBS User Interface 


KBS Executive: Main Screen 


mmmmmmw* 


File Setup Help 


CFD Seal Categories 




ping Seals, - 

Hi 1 '; BiSi-MiS. 



Seal Analysis Codes: KBS User Interface 





Program 






Ind^cat Sear^ilysJf Program’ 


mm&M: 




KBS Executive: Seal Categories Screen 


File Setup Help 


Click On A 


KK 


Seal Analysis Codes: KBS User Interface 


Industrial Code GCYL: Main Screen 


XmSXBBam 

View Help 


Data 


Analysis Options... 

Seal Geometry... 
Operating Conditions... 

s for the GCYL program 


Flow Lines... 

Grid Definition... 

Constant Pressure... 

Seal Fluid... 

Fluid Sources... 

Convergence Criteria... 

Rayleigh Steps... 


Recesses... 


Spot Recesses... 


; CFD Industrial Codes 


Seal Analysis Codes: KBS User Interface 
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Seal Analysis Codes: KBS User Interface 



Seal Analysis Codes: KBS User Interface 


Industrial Code GCYL: Help Screen 



Services Options 


Help 


Data 


Gas Lubricated Cylindrical Seal Analysis 
program (GCYL) may be used to analyze 
a wide variety of seals in a cylindrical 
coordinate reference frame. The solid 
ring seal configurations analyzed by the 
program are shown below: 


Set Title 


CFD Industrial 


Seal Analysis Codes: KBS User Interface 


Industrial Code GCYL: Plotting Program 



Seal Analysis Codes: KBS User Interface 


Scientific Code: Main Screen 



Seal 


Type of Seal 
Analysis Oplions 


Seal Geometry 
Analysis Grid 

Properties 


Initial Conditions 

4 

m 

Solution 

4 

Output Options 

4 


ies for the CFD program 


Inlet Swirl 
Exit Boundary 
Rotor Wall 

Orifii i } . 


• CFD Industrial Codes 
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Seal Analysis Codes: KBS User Interface 



Seal Analysis Codes: KBS User Interface 


Scientific Code: Analysis Options Screen 
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Seal Analysis Codes: KBS User Interface 


Scientific Code: Typical Input Screen 



Seal Analysis Codes: KBS User Interface 
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Seal Analysis Codes: KBS User Interface 



Seal Analysis Codes: KBS User Interface 

Summary and Future Plans 

• KBS implemented to support both OS/2 and Unix users 

• Interactive graphics capability to be added to ease input 

• Selection of a standard plot file format 

• Implementation of database capability 

• New codes to be added as they become available 

• Enhancements based on user feedback 
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BRUSH SEAL BRISTLE FLEXURE AND HARD-RUB CHARACTERISTICS 


Robert C. Hendricks, Julie A. Carlile, and Anita D. Liang 
National Aeronautics and Space Administration 
Lewis Research Center 
Cleveland, OH 44135 


SUMMARY 

The bristles of a 38.1-mm ( 1.5-in.) diameter brush seal were flexed by a tapered, 40-tooth rotor 
operating at 2600 rpm that provided sharp leading-edge impact of the bristles with hard rubbing of the 
rotor lands. Three separate tests were run with the same brush accumulating over 1.3xl0 9 flexure cycles 
while deteriorating 0.2 mm (0.008 in.) radially. In each, the test bristle incursion depth varied from 
0.130 to 0.025 mm (0.005 to 0.001 in.) or less (start to stop), and in the third test the rotor was set 
0.25 mm (0.010 in.) eccentric. Runout varied from 0.025 to 0.076 mm (0.001 to 0.003 in.) radially. The 
bristles wore but did not pull out, fracture, or fragment. Bristle and rotor wear debris were deposited as 
very fine, nearly amorphous, highly porous materials at the rotor groove leading edges and within the 
rotor grooves. The land leading edges showed irregular wear and the beginning of a convergent groove 
that exhibited sharp, detailed wear at the land trailing edges. Surface grooving, burnishing, “whipping,” 
and hot spots and streaks were found. With a smooth-plug rotor, post-test leakage increased 30 percent 
over pretest leakage. 


INTRODUCTION 

High-performance, lightweight engines require compliant seal configurations to accommodate flexible 
interfaces. Thus, in many aircraft gas turbine engines and other turbomachines brush seal systems are 
being proposed to replace labyrinth seals because brush seals are compliant and reliable, leak less, cost 
less, and enhance rotor stability. Brush seals have been the subject of much recent seals research (refs. 1 
to 20). 

A brush seal system consists of the brush and a hardened rub ring and can be linear, circular, or con- 
toured (see ref. 20 for a review). The bristles are oriented to make an angle of 30° to 50° with the inter- 
face, such as the rotor radius for a circular brush. This design allows the bristles to flex when rotor 
excursions occur without significant damage to either the rotor or the seal. 

A typical brush seal configuration, figure 1 (courtesy of Cross Mfg. Ltd. (ref. 1)) consists of (1) a 
backing plate (like a sealing dam), (2) a circumferential or linear set of packed wires (fibers or bristles), 
(3) a pinch plate that serves as a retainer for the brush bristles, and (4) an outside diameter surface that 
fits tightly to the housing (insert in fig. 1). The flexibility of the fibers and implicitly the performance of 
this seal are governed by many factors as expressed in terms of similitude parameters (refs. 9 and 20). 
Among these factors are fiber length and diameter, inclination to the moving surface, surface speed, inter- 
face friction, seal diameter, fluid properties, packing density, modulus of elasticity, backing plate clear- 
ance, and preload or interference fit. 

Typically for a circular brush, the wire or brush materials are superalloys and range from 0.05 to 
0.07 mm (0.002 to 0.0028 in.) in diameter. The bristles are approximately 0.96 mm (0.38 in.) long and 
are aligned at 30° to 50° to the shaft in the direction of rotation. Nominally, there are 98 bristles/mm 
(2500 bristles/in.) of circumference. The interface is characterized by a smooth (4 to 20 rms), hardened 
rub surface on the shaft (e.g., A1 2 0 3 or for short duration the uncoated shaft itself). Because brush seals 
are contact seals with radial interferences ranging from zero to more than 0.25 mm (0.01 in.), ceramic 
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coatings and superalloy materials are often used to enhance life and minimize wear at elevated surface 
speeds, temperatures, and pressure drops. 

Although brush seals show great promise for future applications, it must be acknowledged that brush 
seals are most effective as contact seals and that life and wear rates are major concerns. Whereas bristle 
blowout will cause excessive leakage, bristle loss and debris have a potential for destructive impact on the 
power plant. Thus, the issues of bristle pullout, surface rubbing, bristle wear, and debris are qualitatively 
addressed in this paper. 

The authors are aware that other data exist for ranges of interference fits and other configurations, 
but the results are proprietary. 


APPARATUS 

The “drill press” apparatus was similar to that described in reference 2. In the tests of reference 2 
the 38.1-mm (1.5-in.) diameter brush seal was fixed in a pressure vessel and the rotor was a smooth- 
surface, tapered plug turning at 400 rpm. Leakage data at various interferences and eccentricities have 
been reported (ref. 19). 

In the tests described herein the 38.1-mm (1.5-in.) diameter brush seal was again mounted in the 
pressure vessel that simulated the static housing; but for these tests flutes were machined the length of the 
plug rotor, providing a set of 40 lands and 40 grooves, or a 40-tooth rotor (fig. 2). The lands were 
1.638±0.04 mm (0.0645±0.0015 in.) just above the groove at test 3 and 1.582±0.04 mm (0.0623±0.0015 in.) 
just below the groove at test 1; see rotor sketch on table I. The groove width (fig. 3) averaged 1.397 to 
1.422 mm (0.055 to 0.056 in.) with further dimensions provided in tables I and II. Prior to testing the 
machining tool marks were clear, being axial in the grooves and circumferential on the lands. This rotor 
provided 40 impacts of the brush bristles per revolution and was rotated at 2600 rpm. 


The 38.1-mm (1.5-in.) diameter brush seal (fig. 4) was damaged in a previous series of tests related 
to reference 2. The damaged section, although quite small and having a “chewed” appearance, increased 
seal leakage. The seal could no longer be used for leakage tests but was adequate for the tests herein. 

The rotor was AISI 304 stainless steel, and the brush bristles were Haynes 25 in the annealed condi- 
tion. When stainless steel is rubbed during a machining operation, it tends to change from a “gummy” 
machining material to a surface-work-hardened material. As a result the bristles would be expected to 
rub-machine the stainless steel, and in turn the stainless steel would be expected to abrasively remove the 
bristles. 

In test 1 the interference was set at 0.025 to 0.050 mm (0.001 to 0.002 in.) with the groove depth at 
0.05 to 0.08 mm (0.002 to 0.003 in.). As expected the bristles rub-machined the rotor. 

For test 2 the plug rotor was reset to a portion of the surface that was unrubbed. At that point 
the groove depth was 0.08 to 0.13 mm (0.003 to 0.005 in.). The interference fit between the rotor lands 
and the brush was 0.05 to 0.08 mm (0.002 to 0.003 in.). 

For test 3 the plug rotor was again reset to a portion of the surface that was unrubbed. The groove 
depth was 0.178 to 0.254 mm (0.007 to 0.010 in.). For this test the rotor was initially set 0.33 mm 
(0.013 in.) eccentric, but the rotor rubbed the fence (backing washer) slightly. The fence diameter was 
39.2 mm (1.544 in.). The rotor was then reset to an estimated static eccentricity less than 0.25 mm 
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(<0.010 in.). The dynamic eccentricity was estimated to be less than 1 mm (<0.004 in.). No active 
clearance measurements were made during these tests; these estimates of eccentricities were made from 
post-test photographs and static measurements. These settings provided a test with significant rotor 
impacting and incursion at one portion of the seal and no rubbing contact diametrically opposite to that 
position. 

In tests 1 and 2 the smooth rotor and the 40-tooth rotor were assumed to be interchangeable. The 
repositioning of the stator for test 3 introduced an unaccountable bias that was estimated in order to 
correlate the measured flow rates. For test 3 the smooth-rotor initial static eccentricity was estimated to 
be less than 0.15 mm (<0.006 in.), and after test 3 the estimated static eccentricity was less than 
0.25 mm (<0.010 in.). 


RESULTS 


The results are separated into observations of (1) the brush bristle flexure cycles with associated 
interface damage to the brush seal and the rotor, and (2) the leakage or performance changes. 


Bristle Flexure and Interface Damage 

Visualization of the rotor-brush interface at creeping surface speeds (under 10 rpm) revealed little 
groove penetration. In the impact zone the brush stiffness and the low void did not permit a fully 
deflected or extended set of bristles at the interface. Instead the impact compacted the bristles in the cir 
cumferential direction into the brush and spread the bristles in the axial direction at rate of 40 times per 
revolution. 


Time of testing and brush diameters before and after testing for the three tests are presented in table I; 
additional dimensions are given in table II. The diameters were obtained by inspecting the brush on an 
optical comparator before and after each test. The chewed area and a few stray wires served as reference 
positions for measurements. 


Optical inspection of the grooves cut by the brush into the stainless steel rotor (test 1) showed that 
the cut converged from the leading edge to the trailing edge of the land as the wires (bristles) crossed the 
rotor (fig. 5). Wire grooves were clear cut and debris was evident, as is better shown in the enlargement 
(fig. 6). 


The following groove extents in millimeters (inches) were measured in test 1: 


Width at inlet 1.52 (0.06) 

Width at center 1.3 (0.051) 

Width at outlet 1.27 (0.05) 


For test 2 the inlet region at the land leading edge was extensive and not readily characterized, but a 
general convergence pattern was evident. Similar behavior at the leading edge was noted for test 3. 


In order to corroborate the optical results, profilometer results for a typical tooth of the 40-tooth 
rotor were taken. Wear area and groove depth estimates are provided in table III. Values for the extent 
of the groove in millimeters (inches) are shown in figures 7 to 9 at a resolution of 0.02 mm/division. 
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Width within 0.1 mm of inlet 0.9 (0.035) 

Width at center 0.4 (0.016) 

Width within 0.08 mm of outlet 0.27 (0.011) 


It is evident that the optical values for groove extent were much larger than those from the pro- 
filometer. The problem is the scale used in defining the groove depth. For example, at a resolution of 
0.01 mm/division the width at the center is 1.15 mm (0.045 in.), but at a resolution of 0.02 mm/division 
the width at the center is 0.4 mm (0.016 in.). At the smaller resolution the extent of the scratched inter- 
face generally agrees with the optical values, but at the larger resolution the extent of the scratched inter- 
face is not resolved (i.e., detail is lost). 

For the surface asperity resolution used herein the optical method better defined the extent (width) 
of the damaged interface; the profilometer provided the depth. At one-half the depth resolution only the 
major grooving was defined. 

The first profile, labeled “leading edge,” was taken within 0.10 mm (0.004 in.) of the leading edge 
and shows a broad damage region with deep grooving for test 2 (fig. 8(a)) and test 3 (fig. 9(a)) but a 
minor amount of material damage for test 1 (fig. 7(a)). The material buildup adjacent to the groove of 
test 3 probably occurred during rotor-fence rub. The second profile, labeled “mid section,” was taken 
midway between the tooth leading and trailing edges. The damage of test 2 (fig. 8(b)) and test 3 
(fig. 9(b)) was severe, and a twofold cut has developed in the rotor during test 2. Again moderate 
damaged was noted for test 1 (fig. 7(b)). The third profile, labeled “trailing edge,” was taken within 
5 percent of the trailing edge. The grooving seen in the midsection profile carried through with perhaps 
some sharpness of the features near the trailing edge (figs. 7(c), 8(c), and 9(c) for tests 1, 2, and 3, 
respectively). 

The brush bristle impact at the leading edge left material deposits that were magnetic (i.e., from the 
rotor) and rust color (probably Fe 3 0 4 ) with a spongelike (or cauliflower) appearance (fig. 10); these 
deposits are readily seen at higher magnification (80X) in figure 11. Debris was generated by surface 
machining grooves, “whipping” of the leading edge, burnishing, and sharp trailing edges. The deposited 
materials were fine, porous, “greasy” to the fingers, and readily removed from the rotor; removal from the 
bristles was not straightforward. Standard degreasing cleaned but not thoroughly, and ultrasonic 
cleaning was not attempted. The reasoning was to see if these deposits would inhibit the responsive 
character of the bristles. The debris can affect both the response and the leakage, but neither effect was 
observed in these tests. Further work here is warranted. 

These deposits also indicate rapid wear-in with a long oxidation period for the “machined out” 
material. The materials deposited on the groove wall at the land leading edge (fig. 11) and on the groove 
wall at the land trailing edge (fig. 12) had little or no structure; the defraction spectra were peakless. 

It is speculated that the bristles were dragged across the land, with “machined” material adhering to 
the bristle and then “impacted” off the bristle at the leading edge of the next tooth. Some of the 
materials were deposited within the groove. Black nodule-like debris tended to adhere to the groove wall 
at the land trailing edge. This black material and rust-colored materials formed in the groove at the land 
leading edge. 

Stainless steel work hardens so that the cut grooves were probably harder than the parent stainless 
steel and would wear the annealed Haynes 25 bristles. The smooth grooves, the hot spots, and the hot 
streaking may indicate that a thin layer of stainless steel flowed plastically as it was machined out 
(figs. 13 and 14). Because of the incursive impact of the toothed rotor and the heated interface, the 
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Haynes 25 bristles could lose strength, erode, fracture, or pull out as massive debris. Rut no pullouts or 
massive debris was found after any of the three tests. At higher magnification (20UX) the tips still 
appear intact without fracture, but wear is evident and oxidation debris appears to be well adhered to the 
surface (fig. 15). As further evidence of the bristle wear an examination of the bristle tip surface revealed 
tip grooving (fig. 16), and the severe impacting on the bristles is shown by erratic wear notches on the 
bristle surface (fig. 17). 

The trailing edge of the land was “cut” clean by the brush in all three tests (fig. 18) in stark con- 
trast to the erratic leading-edge surface, which was whipped by the bristles (fig. 19). Of interest is the 
contrast between the land surface cuts. Test 1 surface cuts were a simple wear scar; those of tests 2 and 3 
were multiple grooves with complex surfaces and burnishing (fig. 20). The most rotor damage appeared 
from test 2 and the most brush damage from test 3, where the rotor was set eccentric. 

These tests, although preliminary and only qualitative, begin to mitigate the fear of brush seal dis- 
integration through bristle flexing as over lxlO 9 cycles were sustained without failure, fracture, or pull- 
out. However, the required flexures are at least an order of magnitude higher with parameters such as 
surface speed, temperature, pressure, and materials to be considered. 

Total flexures - 222 hr x 60 min/hr x 40 teeth x 2600 rpm = 1.38xl0 9 

Required eccentric shaft flexures = 10 000 x 60 x 1 (flexure/rotation) X 20 000 rpm = 12xl0 9 
Required rotor disk flexures = 50 X eccentric shaft flexures 

However, bristle flexures raise an equally troubling concern over seal life, because brush seals do 
wear out. Once these seals begin to reach line-to-line contact, their leakage can be equivalent to that of 
an advanced labyrinth seal. The sealing margin and competitive edge of the worn brush seal begin to 
fade. New competitive (lower leakage) configurations for labyrinth, damper, honeycomb, feltmetal, and 
spiral-groove seals are under investigation. It is clear that long-duration testing at elevated surface 
speeds and working fluid temperatures are required. 


Correlation of Leakage Data 

Although not the primary objective of this experiment, overall changes in brush leakage were 
estimated from flow checks before and after testing. In order to determine these leakages the 40-tooth 
rotor was replaced with a smooth-surface, tapered rotor. Runout errors resulting from rotor interchange 
were unresolved as were those associated with the static eccentricity of test 3. Measurements characteriz- 
ing the rotor and brush before and after testing are given in table I. The average depth of the brush-cut 
groove as well as the estimated clearances are given in tables II and III. 


Leakage is characterized in figure 21 in terms of volumetric flow rate as a function of pressure drop 
across the brush seal before and after each of the three tests. Both pretest and post-test results are 
provided in the same figure. Because this brush seal was damaged (see APPARATUS section), absolute 
leakage measurements would require weighing, but the relative leakage should be accurate. The interfer- 
ence fits for the brush seal leakage data for pretests 1 and 2 were nearly the same, resulting in corres- 
ponding leakages. While taking data it was found that the brush would stiffen and the pressure drops 
would increase. Data points illustrating hystersis (typical in brush seals) are shown. After correcting the 
post-test 2 data for clearance, these test results agreed with those of post-test 1. 
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Setting the rotor eccentric in test 3 proved a major problem in cross correlating the leakage results. 
The estimated initial static eccentricity for the smooth rotor was less than 0.15 mm (0.006 in.), and the 
smooth rotor and the 40-tooth rotor were assumed to be interchangeable at the same spindle loading. 
However, the rotor rubbed the fence (backing washer) slightly, requiring an initial static eccentricity of 

0.36 mm (0.014 in.) and implying a difference in spindle loading. The rotor was reset to an estimated 
eccentricity of less than 0.025 mm (0.010 in.). Post-test photographs indicated that the dynamic eccen- 
tricity was <0.10 mm (<0.004 in.) and clearly illustrated the fence rub (fig. 22). 

For test 3 the repositioning of the stator and the differential spindle loading introduced an un- 
accountable bias that was difficult to estimate in correlating the leakage results. From the data of 
reference 19 a relation was found for the change in pressure as a function of eccentricity at a fixed volu- 
metric flow rate. Using this relation and corrections for clearance and assuming a pretest static eccentri- 
city of 0.15 mm (0.006 in.) and a post-test static eccentricity of 0.25 mm (0.010 in.) show that the results 
of test 3 were overcorrected by 20 percent with respect to the results of tests 1 and 2. Future testing 
requires instrumentation to overcome these positioning errors. Nevertheless, these leakage data indicate 
that under conditions of severe brush and rotor wear the brush seal leakage increased 30 percent. And, 
although brush seal performance degraded, the brush seal did not fail. 


SUMMARY OF RESULTS 

In three separate tests with a 40-tooth tapered stainless steel rotor operating at 2600 rpm and a 
38.1-mm (1.5-in.) diameter brush seal with 0.07-mm (0.0028-in.) diameter annealed Haynes 25 bristles set 
at a nominal 0.076-mm (0.003-in.) radial interference for each test, the following results were obtained: 

1. The bristles withstood over lxlO 9 cycles without pullout, fracture, or massive debris generation. 

2. Rotor grooving up to 0.076 mm (0.003 in.) in depth radially with erratic “whipped” leading-edge 
surfaces followed by convergent grooving to a clean-cut trailing edge was commonplace for each of the 
three tests. 

3. Most of the debris generated was a fine black material that appeared amorphous, but the rust- 
colored materials were iron rich and magnetic, implying Fe 3 0 4 . The debris was “cauliflower” in form and 
highly porous with low adhesion, except for that which was fine enough to adhere to the bristles. Those 
fines were not readily dislodged. Nonuniform fines (or oxidation) adhering to the bristles tended to sepa- 
rate the bristles, increasing porosity, and would enhance leakage paths. 

4. Radial bristle losses up to 0.2 mm (0.008 in.) were demonstrated, which if left uncorrected would 
lead to equivalent or higher leakages than those of labyrinth seals. Bristle loss at elevated surface speeds 
and temperatures requires further study. 

5. Generated debris can impair bristle motion and alter leakage, but within the limitations of this 
experiment these considerations were not a problem. They remain as issues to be resolved. 

6. Under conditions of severe rotor-stator interface damage, the brush seal leakage performance 
degraded 30 percent, but the seal did not fail. 
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TABLE I.— DIAMETRAL CHANGES AND TEST CYCLE TIMES FOR 
38.1-mm (1.5-in.) DIAMETER BRUSH SEAL 
[Seal fence inside diameter, 1.543±0.0005 in.] 



Position 

Average 


Vertical 

Horizontal 

Vertical 

Horizontal 



Rotation, deg 



0 

45 



Brush diameter (from optical comparator inspection), in. 

Before test 1 

1.494 

1.4918 

1.4926 

1.4945 

1.4932 

After test 1* 

1.497 

1.496 

1.495 

1.495 

1.4958 

Change, in. 

0.003 

0.0042 

0.0024 

0.0005 

0.00255 

Before test 2 

1.495 

1.495 

1.497 

1.496 

1.4958 

After test 2 b 

1.5058 

1.5061 

1.5041 

1.5037 

1.5049 

Change, in. 

0.0108 

0.0111 

0.0071 

0.0077 

0.00918 

Before test 3 

1.5058 

1.5061 

1.5041 

1.5037 

1.5049 

After test 3 C 

1.504 

1.514 

1.512 

1.5064 

1.5091 

Change, in. 

-0.0018 

0.0079 

0.0079 

0.0027 

0.00418 


*Total test time, 70.3 hr. 
b Total test time, 43 hr. 
c Total test time, 109 hr. 
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TABLE II.— ADDITIONAL TEST AND 
GEOMETRY INFORMATION 
[Seal fence inside diameter, I.543±0.0005 in. 
Dimensions are in inches.] 

(a) 40-Tooth rotor 



— Tooth-to-tooth diameter 



Tooth gap, 0.05510.002 


Position 

Tooth-to-tooth diameter 

Gap-to-gap diamter 

Top of rotor 
Bottom of rotor 

1.523±0.0005 

1.491±0.003 

1.485±0.0005 

1.483L0.003 


Test 

case 

Top-of-brush 
wear surface 

Bottom-of-brush 
wear surface 

Average 

1 

0.605 

0.632 

0.618 

2 * 

.461 

.508 

.484 

3 

.325 

.303 

.344 


(b) Smooth rotor 



l 3 


l 2 


L 


Test 

L a 

Smooth-rotor 

Free brush 

Concentric 

Static 

case 


diameter 

diameter 

radial clearance 

eccentricity 

Before test 1 

LI 

1.4992 

1.4932 

-0.0030 

0 

After test 1 


1.4992 

1.4958 

-.0017 

0 

Before test 2 

L2 

1.5046 

1.4958 

-0.0044 

0 

After test 2 


1.5046 

1.5049 

.0001 

0 

Before test 3 

L3 

1.5105 

1.5049 

-0.0028 

0.013 

After test 3 


1.5105 

1.5091 

-.007 

.013 


Corresponds to equivalent axial positions of 40-tooth rotor. 
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TAIILE 111. RRORILOMETER RESULTS KOR 40-TOOTH 
ROTOR, TESTS 1, 2, AND 3 


Location 

Test 

Peak groove 
depth, 
pm 

Estimated 
average depth, 
H m 

Wear area, 
/im 2 

Leading edge 

1 

57 

22 

24 910 


2 

130 

95 

141 680 


3 

130 

70 

114 370 

Midsection 

1 

30 

17 

7 652 


2 

70 

35 

37 106 


3 

70 

27 

56 493 

Trailing edge 

1 

25 

12 

6 610 


2 

65 

35 

48 649 


3 

67 

30 

56 169 



SttMfav 


Figure 1 . — Circular brush seal. (Courtesy of Cross Mfg. Ltd.) 
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Figure 2. — Geometry of tapered, 40-tooth rotor. 
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Figure 3.— Rotor prior to tests (40 lands). 



Figure 4 —Brush seal prior to tests (270° mark). 
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Figure 5. — Rotor surface after test 1 . (Arrow shows direction of rotation.) 



Rotation 


C-92-031 27 


Figure 6. — Rotor land and groove after test 1 . 
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Groove depth, jim 



(c) Trailing edge. 



(b) Midsection. 



8 10 12 
Groove width, mm 

(a) Leading edge. 


Figure 7. — Profilometer traces for 40-tooth rotor (test 1). 





Groove width, mm 
(a) Leading edge. 

Figure 8.— Profilometer traces for 40-tooth rotor (test 2). 
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Groove depth, fim 





1 2 3 4 5 

Groove width, mm 

(a) Leading edge. 


Figure 9. — Profilometer traces for 40-tooth rotor (test 3). 



Figure 1 0.— Rotor land with debris and rub scars (test 2). 
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Figure 1 1 .—Rotor leading-edge debris formation (test 2). 



Figure 12. — Rotor trailing-edge surface grooving and debris (test 3). 
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Figure 1 3. — Rotor leading edge, 


surface machining, and hot spots (test 3). 
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Figure 1 5. — Brush bristle tips with debris (test 1). 



Figure 1 6. — Brush bristle tip wear patterns (test 2). 
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Figure 1 8. — Rotor trailing edge (test 3). 
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Figure 1 9. — Ftotor leading edge (test 2). 




Figure 20. — Rotor trailing edge (tests 1 , 2, and 3). 
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Figure 21 . — Pre- and posttest leakage results for smooth rotor. 


60 


70 



Figure 22. — Fence damage sustained during test 3. 
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SUMMARY 

A split-ring brush seal was fabricated, installed between two labyrinth-honeycomb shroud seals, and 
tested in the fourth-stage turbine of a T-700 engine. The annealed Haynes 25 bristles rubbed directly 
against the nonconditioned, irregular Rene 80 turbine blade shroud surface. A total of 21 hr of cyclic and 
steady-state data were taken with surface speeds to 335 m/s (1100 ft/s) and shroud temperatures to 620 °C 
(1150 °F). Wear appeared to be rapid initially, with an orange flash of hot brush fragments during the first 
engine startup, to minimal after 10 hr of operation. The brush survived the testing but experienced some 
bristle pullouts and severe bristle wear; some turbine interface wear and possible material transfer was 
noted. Future design concerns center on tribological behavior at the interface with or without lubricants. 


INTRODUCTION 

Engine testing of brush seals has been reported (e.g., Rolls Royce (ref. 1) and Allison (ref. 2)) that 
demonstrated performance increases relative to labyrinth seals. These brush seal systems had smooth 
rotor interfaces (<25 rms) and operated at moderate temperatures and surface speeds. Even though 
these tests were successful, concern over catastrophic failure of the brush, such as a loss of bristles when 
subjected to high surface speeds at elevated temperatures, has not been resolved. 

The objectives of this program were first to demonstrate that a well-designed and manufactured 
brush seal could survive the “pounding” of an irregular rotor surface without catastrophic failure, second 
to illustrate the concept of running a combined brush and labyrinth seal system, and third to acquire 
metallographic data on bristles subjected to such an environment. 


EXPERIMENTAL BRUSH CONFIGURATION AND INSTALLATION 

A cross-sectional view of the split-ring brush seal configuration is illustrated in figure 1. A major 
problem in designing retrofit seals is consistency of hardware measurements, and for this experimental 
engine the seal was crafted to fit. 

In order to fit into the existing T-700 engine fourth-stage turbine shroud, the design was required to 
fit into a radial clearance about half that normally used for brush seals. The 0.07 1-mm (0.0028-in.) 
diameter, Haynes 25 bristles were angled 43° to 50° to the interface with about 2500 per inch of circum- 
ference (98.4 per millimeter of circumference). The backing washer (or fence) was angled 19° to match 
the slope of the turbine shroud. The design clearance was —0.51 mm (—0.02 in.) but could range to 
— 1.27 mm (-0.05 in.) diametral (the uncertainty reflecting that of the engine geometry) with an outside 
diameter of 333.9 mm (13.146 in.) and an inside diameter of 322.3 mm (12.690 in.). 


Figure 2 shows several views of the split-ring brush seal, illustrating the joint and the restraining pin 
hole. Figure 3 shows the unusual brush cross section that was crafted to Fit the turbine shroud between 
the two labyrinth-honeycomb seals. Figure 4 shows a closeup view of the installed brush; and Figure 5, an 
overall view of the installed brush. Figure 6(a) shows the pressure tap locations in one of the shrouds. A 
thermocouple was installed in each shroud (Fig. 6(b)); three of the four were functional. Figures 6(c) and 
(d) show views of the shroud ring and instrumentation lines. 

Assembly of the turbine with the shroud required forcing the brush past the upstream labyrinth 
tooth without any visual or instrumented guidelines; see Figure 7 (power turbine). Forcing the brush over 
the labyrinth tooth spread the bristles axially into the upstream direction. This type of spreading alters 
the bristle packing conFiguration, but the extent of alteration and the degree of spreading are unknown. 
Post-test results indicated that perhaps two to three upstream bristle rows remained spread with possibly 
one downstream bristle row in disarray. Still the brush was resilient because the remaining rows 
appeared to be in position. The flexibility of a brush seal and the abuse it can withstand appear to be 
signiFicant. 

With the brush installed, the turbine shaft was difFicult to rotate, requiring 14.7 N-m (130 in.-lbf) of 
torque. This was a major concern because heat generation could be sufFicient to melt the materials at the 
interface. The geared tooth rotor results (ref. 3), including material smearing, cutting, and local hot 
spots, indicate that high heat loads and temperatures could be present but would be conFined to the inter- 
face, with the lowest heat-sinking element (the bristles) absorbing the energy. Therefore, the bristles 
would fail, but the effects on the power turbine should be benign. 


ENGINE OPERATIONS 

The T-700 turbine section was assembled and the brush seal test was “piggybacked” on the break-in 
of the engine. Operations consisted of the standard break-in procedures with data taken under both 
steady and cyclic conditions. The engine was operated a total of 21 hr, including break-in, steady state, 
and 10 hr of cycling between ground and flight idle (4-min ground idle and 5-min flight idle). Turbine speeds 
were 10 000 and 20 000 rpm, and average fourth-stage turbine shroud temperatures were 455 and 566 °C 
(850 and 1050 °F), respectively. Maximum shroud temperatures were limited to 621 °C (1150 °F). The 
turbine inlet temperatures were about 139 deg C (250 deg F) higher. The pressure drop measurements 
across the brush were up to 0.007 MPa (1 psia) and varied from shroud to shroud. An assessment of the 
effect of the brush seal on engine performance was inconclusive and remains to be investigated further. 
Neither radial nor axial positions of the rotor were monitored, but such position sensors should be an 
integral part of the engine dynamics. 

Because of concern over the 14.7 N-m (130-in.-lbf) installed torque that was required to rotate the 
power turbine shaft, after about 10 hr of engine operations the compressor and the power turbine were 
decoupled. The turbine shaft turned freely but not in reverse. The turbine assembly has 50 shrouded 
blades with irregularities (radial, to 0.229 mm (0.009 in.); circumferential, to 0.076 mm (0.003 in.); and 
axial, to 0.051 mm (0.002 in.)) representing protrusions into the brush and the spaces between the blade 
pairs. It is not known how many cycles were required to “free the bristles,” but at 10 000 rpm and with 
24 irregular asperities impacting each bristle (4000 impacts/s at a surface speed of 168 m/s (550 ft/s)) it 
is assumed that brush break-in was rapid. The annealed Haynes 25 bristles rubbed against hardened 
Rene 80 blades and probably wore rapidly during the initial stages of engine break-in. 

Furthermore, and of signiFicance to engine designers, a flash was noted upon initial engine ignition 
that was concluded to be expulsion of brush fragments. This is important because critical components 


must be protected against initially high levels of debris generation. Analysis of these and other fragments 
showed severe oxidation with some degree of stiffness remaining. These fragments are not passive debris; 
they can cause damage to critical components. The only debris noted in the gear-tooth rotor study 
(ref. 3) was a “lubricant powder.” Thus, surface speed, rotor roughness, and brush construction play a 
major role in determining the spectrum of debris generation. 


RESULTS AND DISCUSSION 

In any potentially destructive situation one attempts to preserve the critical components. Brush 
bristle wear would degrade performance, but failure of the turbine or a shaft (depending on seal location) 
could result in engine loss, perhaps catastrophic. So the brush seal becomes the sacrificial component, not 
only in this case but when running against a coated shaft. The brush seal would then be replaced at dis- 
crete intervals, such as during overhauls. 

For these tests the brush rubbed the turbine blade shroud asperities smooth and did provide a dis- 
tinct wear track, perhaps through transfer of material. However, no direct damage was ascribed to the 
turbine blade shrouds. The engine was immediately returned to service. 

Before proceeding to discuss post-test results, we return to the installation of the brush into the 
power turbine shroud. Figures 4, 5, and 7 show views of the installed brush seal. The brush was 
designed for an interference of 0.51 mm (0.020 in.) diametral. The actual interference could not be deter- 
mined, but estimates of the pretest brush clearance were —0.51 to —1.27 mm (—0.02 to —0.05 in.) 
diametral. Measurements of the shrouds differed as did those for the rotor. Although the differences 
were only a few thousands of an inch (mils), they represented a significant percentage of the clearance 
gap. It was also determined that individual blade sets could have a step change of 0.229 mm (0.009 in.) 
from one blade set to the next. These surface irregularities are shown in the post-test photograph 
(fig. 8). 

Upon initial engine startup, an orange burst was noted and was assumed to be expelled brush 
bristles (i.e., those that were inadequately attached, were pulled out by rotor irregularities, or were 
embedded within the blade row gaps during the blind installation and “snapped” or “yanked” or “bent” 
aside during engine break-in). Residual bristles from the exhaust were photographed (fig. 9). Although 
they appeared to be highly oxidized and stressed, they were curly and still wirey. The number of residual 
bristles decreased with operation until none were noted. It is assumed that at this time the bristles and 
the rotor were in nearly line-to-line contact (i.e., rubbed in). Precooling the bristles during the initial 
rub-in may mitigate bristle loss and the orange burst. 

Turbine speeds to 20 000 rpm and shroud temperatures to 620 °C (1150 °F) were commonplace. 
These conditions provided an interface speed of 335 m/s (1100 ft/s) or a temperature-velocity (TV) 
product of over (1100) 2 with the global target of about (1500) 2 (in U.S. customary units, feet per second 
per degree F). 

Further assessment of figure 8 shows that the leading edges of the blade sets were polished. Some 
brush wear was noted and expected because Rene 80 is hard relative to the annealed Haynes 25 and the 
heat-sinking capacity of the brush is very small relative to that of the rotating blade sets. The surface 
irregularities at the rubbing interface can be seen. Although there was some evidence of material trans- 
fer, no metallurgical samplings or rotor measurements were taken because of the tight program schedules 
for the engine operations. This was a major error; however, if there comes a time when the rotor can be 


looked at, some of the transferred materials may be found still embedded in the rotor even after it has 
undergone other program tests for General Electric Co. 

Figure 10 illustrates bristle spreading after testing, with a central core of bristles rubbed (probably 
clip cut and worn to shape). Although the environment was hostile, the brush did not disintegrate, but 
bristle pullout could be noted in a few places. Some upstream bristles (two or three rows) show spread- 
ing and perhaps one row downstream (toward the exhaust). The remainder of the bristles show wear (or 
cutting plus wear). Detailed estimates of the bristle stubble heights are given in figure 11(a) as taken 
from the set of photographs in figure 11(b) corresponding to positions (joint, J+l, J+1.5, J+2, J-f-2.5, 
J+3, J+3.5, J+4, and J+4.5) shown in figure 11(a). The last two photographs in figure 11(b) correspond 
to the minimum stubble height, where the rotor actually rubbed the fence (position J+3); the view is 
looking toward the bristle stubbles to show the rubbed fence. Figure 11(c) shows wear track and fence 
rub, and figure 11(d) shows the joint wear track. 

The following measured diametral parameters in millimeters (inches) were used to establish seal 
clearances: 


Pretest brush diameter (no taper) 
Post-test brush diameter (tapered) 

Differential 

Pretest rub interface diameter . . 

Brush pretest interference 

Labyrinth cold clearance 

Bright region diameter (fig. 8(a)) 

Differential 

Stubble height (fig. 11(a)) 

Possible material transfer 

Possible engine eccentricity .... 
Blade shroud height variation . . 


321.79 (12.669) 

324.15 (12.762) 

2.36 (0.093) 

322.30 to 322.81 (12.689 to 12.709) 

0.51 to 1.27 (0.02 to 0.04) 

2.29 to 2.46 (0.09 to 0.097) 

323.42 (12.733) 

0.74 (0.029) 

0.71 (0.028) 

0.25 to 0.5 (0.01 to 0.02) 

0 to 0.8 (0 to 0.03) 

0 to 0.23 (0 to 0.009) 


The average diametral clearance estimates in millimeters (inches) are as follows: 


Pretest brush clearance -0.51 to -1.3 (-0.02 to -0.04) 

Post-test brush clearance 0 to 0.8 (0 to 0.03) 

Post-test labyrinth clearance 2.3 to 2.5 (0.09 to 0.10) 


The engine parameters are as follows: 

335 (1100) 
620 (1150) 

. ( 1100) 2 


Turbine speeds to 20 000 rpm, m/s (ft/s) 

Maximum turbine shroud temperatures, °C (°F) 

Temperature-velocity product (approximate; in U.S. customary units ft/s-°F; 
future target, (1500)*) 


Details of the bristles (fig. 12) show an ingrained wear pattern that is characteristic of a high spot in 
the rotor which cuts a shallow groove as it wears in. Furthermore, the rotor ran eccentric with respect to 
the seal either during cycling or steady state or both. Plots of bristle measurements show rubbing of the 
seal ring at 180° to the pinned point and a step across the joint where the end of the split seal rubbed 
the rotor. 
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The blind installation and operation made it difficult to assess intermediate states of wear or the 
health of the brush or the turbine shroud/brush interface. Although accelerometer measurements were in 
bounds, no detailed information concerning the orbital dynamics was available. 

Metallographic results (fig. 13) illustrate some material migration along the bristle and material 
transfer both from and to the surface. Material smears seem to be in line with a softer material rubbing 
a harder material even to the point of melting. Thus, one would conclude that the interface became very 
hot, but the interface followed the classic Block slider problem, where the penetration depth (radial) is 
very small. Thus, the thermal effect at the interface appeared to be topical, wearing the brush bristles 
until line contact with the rotor and smear transfer of both materials occurred. Such evidence is 
presented in figure 14, for an individual bristle that was clipped from the brush. The clipped end was 
removed from the interface, and the worn end is representative of interface materials. There are material 
smears, oxidized spots, and apparent pits. It is not clear how the irregularities of the interface affected 
these results, but it is clear that materials were transferred. 

Metallographic analysis of the sectioned brush will be completed and reported later. 

Derby and England (ref. 4) reported excellent brush bristle and coating wear for Alloy A (Ni-Cr-Al- 
base superalloy solid-solution strengthened) bristles and Triboglide coating. Alloy A is being used in gas 
turbine hot spots and develops a tenacious Cr 2 0 3 and A1 2 0 3 , yittria-modified oxide layer. Triboglide is a 
chromium carbide (CrC) containing additives of 12 wt% barium and calcium fluoride solid lubricants. 
Triboglide is based on the work of Harold Sliney at NASA Lewis but has no silver additive. The tests 
were performed with 1200 °F air. 

Atkinson and Bristol (ref. 5) report better wear for a cobalt alloy rubbing against chromium carbide 
at room temperature, but their high-temperature (480 °C) result shows nearly equivalent wear for either 
the cobalt- or the nickel-base alloy. However, the Co-alloy/CrC combination had less leakage under 
dynamic conditions and better wear at room temperature. The tests were conducted to simulate a CT7-9 
compressor discharge seal. The brush was 129 mm in diameter and of standard Cross Mfg. construction. 


It is apparent that the composition of both the coating and the bristles needs to be characterized 
with respect to the working fluid, the operating conditions, and the component life requirements. The 
importance of surface conditions must be emphasized. Wear decreases after operation because the brush 
rubs a smoother surface and bristle wear decreases line loading. 

Obviously, tribological pairing is important and references 4 and 5 present a good initial look into 
these problems. Limitations on speed, temperature, and preloading have to be established. Furthermore, 
limitations on surface asperities have to be established for expedient or commercial engines even though it 
has been demonstrated that direct bristle rubbing of a smooth shaft could be acceptable for expedient 
engines (ref. 6). 

SUMMARY OF RESULTS AND CONCLUSIONS 

A split-ring brush seal installed between two labyrinth-honeycomb shroud seals was tested in the 
fourth-stage turbine of a T-700 engine. The following results were obtained and conclusions drawn: 
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1. Properly designed brush seals have sufficient integrity to withstand highly irregular surface opera- 
tions at surface speeds to 335 m/s (1100 ft/s) and shroud temperatures to 620 °C (1150 °F) with a non- 
centered turbine orbit during steady and cyclic loading. 

2. Upon initial engine startup, bristle debris can be expected mostly in the form of fines and some 
larger elements. The nature and amount of debris should depend on the construction, the surface charac- 
teristics (e.g., asperities), the temperature and velocity of the interface, and the bristle preload. Critical 
components must be protected. 

3. The post-test clearance was estimated to be line to line to 0.8 mm (0.03 in.), indicating a well- 
worn but still functional seal configuration. Accurate determinations of the rotor and stator dimensions 
along with dynamic displacement measurements are necessary for assessing bristle wear characteristics 
and seal clearances. A plot of bristle stubble height versus circumferential position revealed some charac- 
teristics of the rotor and engine operations. 

4. Wear is expected to be initially rapid, then steady, and subsequently decreasing with time of 
engine operation. Cyclic operations cause more rapid wear of the bristle/rotor interface. Material trans- 
fer, smearing, and pitting of the interface are commonplace. 

5. Installation torques can be high, but rub-in torques are low. Rotor reversals are not permitted. 

6. Tribological pairing is important, and limitations on speed, temperature, preload, and asperities 
have yet to be established even though direct rubbing of a smooth shaft of an expedient engine appears 
plausible. 

7. Although high installation torques (14.7 N-m; 130 in.-lbf) probably contributed to high bristle 
wear, the effects on the power turbine were benign even though heat generation and shear were suffi- 
cient to transfer materials at the interface. The Rene 80 is hard and has large heat capacity, and the 
Haynes 25 bristles are annealed and have small heat capacity. Therefore, the brush bristles failed first. 

8. The pressure drop measurements across the brush seal were up to 0.007 MPa (1 psia), but the 
effect of the brush seal on engine performance was inconclusive and requires further assessment. Radial 
and axial rotor sensor position monitoring is recommended. 

9. Metallographic studies of the brush and rotor are being completed. 
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Figure 1 .-Cross-sectional view of split-ring brush seal configuration. (Brush design with -0.020 in. interference Dimensions are 
in inches.) 
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(a) Overview of brush seal. 



C-92-01113 


(b) Edge view of brush seal. 



(c) Split end view. 

Figure 2. — Split-ring brush seal. 
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Figure 3. — Brush seal cross section crafted to fit turbine shroud between two labyrinth- 
honeycomb seals. 



Figure 4.— View of installed brush seal. 
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Figure 5. — Overview of installed brush seal and power turbine housing. 
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(a) Pressure tap hole. 



(b) Thermocouple mount and pressure tap lines. 



(c) Shroud ring with three seal segments and instru- 
mentation lines — downstream view. 



C-92-00774 

(d) Shroud ring with three seal segments and instru- 
mentation lines — upstream view. 

Figure 6. — Pressure tap and thermocouple locations 
on shroud seal. 
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Figure 7. — Schematic of power turbine. 



Figure 8.— Fourth-stage turbine after testing, showing polishing of leading edges. 
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(a) Typical bristle tips. 





(b) Surface features. 


Figure 9. — Typical brush seal debris found in engine exhaust duct. 
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(b) Continued. 
Figure 9. — Continued. 
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(b) Concluded. 
Figure 9. — Concluded. 



Figure 10. — Post-test view of brush seal showing bristle spreading and core of rubbed 
bristles. 
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J ♦ 2. J ♦ 2.5. J ♦ 3. 

(b) Bristles at backing washer. 

Figure 1 1 . — Post-test bristle stub height and backing washer rub. Bristle diameter, 0.0028 in. 
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(a) Bristle height versus circumferential position. 
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J ♦ 3.5. J + 4. J + 4.5. 



J + 3 (rotated view) 40x. J + 3 (rotated view) 23x. 

(b) Concluded. 

Figure 1 1 . — Continued. 
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(c) Wear track and fence rub. 



(d) Joint wear track. 
Figure 1 1 . — Concluded. 
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Figure 1 2. — Details of brush seal wear pattern. 
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(b) Analysis of bristle 1 . 

Figure 13. — Post-test metallographic results for single bristle. 
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(c) Bristle 1 at increased resolution. 
Figure 1 3. — Concluded. 
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Figure 1 4.— Post-test metallographic results tor bristle tip from wear track. 
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National Aeronautics and Space Administration 


Presented August 6, 1992 
at Brush Seals Systems Workshop 
Lewis Research Qenter 
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NASA 

Nallona? Aerforrmjttes and MSFC SOdl Activlti©8 

Space Admimatratton 

Overview 

• Highly-instrumented SSME HPOTP with hydrostatic bearing on 

Technology Test Bed 

• Cryogenic fluid-film bearing tester 

• Damping seal rotor support 

• Damping seals for turbomachinery 

• Experimental verification of rotordynamic analysis 
« In-house CFD analysis 


NASA 

National Aeronautics and MSFC Seal Activities 

Space Administration 


Highly-Instrumented SSME with HPOTP hydrostatic bearing on 
Technology Test Bed 

Objective: Demonstrate feasibility of hydrostatic bearings for high 
pressure cryogenic turbopumps 

• Joint Rocketdyne IR&D/MSFC activity 

• Retrofit HPOTP pump-end with hydrostatic bearing for ball bearings 

• Operated engine at multiple power levels with variable ramps and 

multiple LQ inlet conditions ' 

• Bearing accumulated 723 seconds in 7 tests 

• Rotordynamic stability maintained during all operational phases 

• Proximity probe data indicated shaft contacted wall during both start 

and shutdown 

• Post-test bearing wear data indicated minimum wear of 0.00109 inches 

over 7 test series 
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NASA 

National Aeronautic* and MSFC Seal Activities 

Space Administration 


Cryogenic Fluid Film Bearing Tester 


Objective: Procure externally fed fluid-film bearing tester to test 

Earth-to-Orbit sized turbopump bearings in cryogenic fluids 
allowing for evaluation of bearings in realistic turbopump 
environment 

• Should obtain repeatable and usable data 

• Fluids of LHj,, LQ , KjO, and surrogates 

• Only facility available to evaluate impartially and fully fluid film 

bearing concepts in cryogenic fluids 

• Knowledge to be gained: 

0 Important rotordynamic coefficients 
0 Material characteristics, tribology of rubbing components 
0 Basic load capacities 


NASA 

National Aeronautics and MSFC Seal Activities 

Space Administration 


Damping Seal Rotor Support 


Objective: Obtain test rig to verify internally fed damping bearing 
predictions, validate rotordynamlo codes 


• Water as working fluid 

• Experimentally determine dynamic coefficients and compare to 

in-house analytical predictions 

• Include parametric changes to mechanical dimensions of seal and 

bearings as well as to fluid inlet and exit conditions 

• Develop understanding of bearing performance as a function of 

parametric changes and analysis 

• Preliminary design complete 
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NASA 

NUkrmt A#rcn« u tim and M3FC Seal Activities 

Spw# Admmtir nH i o n 


Damping Seats for Turbomaohirrery 

Objective! Verify" dynamlo coefficient predictions of damping seals 


• Rig used to determine dynamic coefficients of axially fed damping seals 

0 8 ' seal configurations tested 

• Water as working fluid 

• In-house bulk-flow codes closely match most of experimentally 

determined dynamic coefficients 

• Performed at Wyte Laboratories 


NASA 

National Aeronautics and MSFC Seal Activates 


ExperfmentaTVeHtleatfon of rtotordynamic Analysis 


Objective: Experimentally evaluate damping rotor supports for 

Improvements In suppression ot rotor whin in turbopump 


« Tests of 2 externally' and 2 Internally fed damping bearings In L0 2 surrogate 
fluid 

• Provide head-to-head comparison between 4 damping bearing 

configurations 

• Performed under contract to Rocketdyne 

• Best combinations of bearing characteristics is Candidate for replacing 

ball bearings In turbine end of HPOTP 
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f\J/\SA 

National Aeronautics and MSFC Seal Activities 

Space Administration 


In-house CFD Analysis 

Objective: Mature CFD for fluid film applications to feed and support 
one-dimensional flow models 

• Three-dimensional analysis on 60° slice of bearing 

• Single-phase incompressible Navrer-Stokes analysis assuming 

constant-y K, 

• Rotational Reynolds number based on annulus width » 4.8 x 10* 

• Multi-block solution in progress with FDNS code (Chen) 

• K-e turbulence models 


•CBNCTFV 



l.l*S 
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NASA 

Nation* * Aeronautics tnd MSPC Seal Activities 

8 pace AdmtntatrtHon 


ln*house CPD Analysis 
Deiermfnatlon of Rotordynamlc Coefficients 


• Obtained from axisymmetrlc CFD results using REFLEQS (CFDRC) 

0 Zeroth-order solution 

• Presently using perturbation method of Dr. Erian Baskharone 

(Texas A&M University) and Dr. -Steve Hensel (former Texas 
A&M graduate student) to calculate coefficients 

• Currently being applied to ATD HPOTP balance piston cavity 
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Summary 


• technology Test Bed of HPOTP hydrostatic bearing was suooessful 

• OAST tasks support range of seal and bearingconfiguratlons 

0 Internally fed 
0 Externally fed 
0 Axially fed 

0 Emphasis on rotordynamlc coefficients and seal performances 

• CfD analysis In progress for seal-flow evaluation 
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AIR FORCE BRUSH SEAL PROGRAMS 

Capt Connie Dowler 
Aero Propulsion & Power Directorate 
Wright Laboratory 

Wright-Patterson Air Force Base Ohio 


Aggressive pursuit of increased performance in gas turbine engines 
is driving the thermodynamic cycle to higher pressure ratios, 
bypass ratios, and turbine inlet temperatures. As these parameters 
increase, internal air system and resultant thermodynamic cycle 
losses increase. This conflict of reducing internal airflows while 
increasing thermodynamic efficiency and performance is putting more 
emphasis on improvements to the internal flow system. One 
improvement that has been and continues to be pursued by the Air 
Force for both man-rated and expendable turbine engine applications 
is the brush seal. This presentation briefly describes both past 
and current brush seal research and development programs and gives 
a summary of demonstrator and developmental engine testing of brush 
seals . 


r 

OUTLINE 

• 

PAST R&D PROGRAMS 

• 

CURRENT R&D PROGRAMS 

• 

DEMONSTRATOR ENGINE BRUSH SEALS 

• 

F1 19 BRUSH SEALS 

• 

SUMMARY 


... J 
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PAST AF BRUSH SEAL R&D PROGRAMS 


• ALLISON: HIGH TEMPERATURE BRUSH SEALS 

• Inconel X-750 bristles wear better than Haynes 25 in hot dynamic 
contact with chrome-carbide Journal 

• Initial build-up Interference effects static leakage 

• Brush seal with hot running clearance has less leakage than 4-knife 
lab seal at same pressure ratio 

• Report #: WL-TR-91-2005 

• TELEDYNE: BRUSH SEAL DESIGN 

• Brush seals provide a factor of 3 or more reduction In leakage flow 
over conventional lab seals 

• Brush seals retain a significantly reduced leakage over lab seals for 
time periods of limited-life engines 

• Brush seals can survive shaft excursions of over 0.025 inches 
without any performance loss 

• Report#: WRDC-TR-90-2123 


< s 

PAST AF BRUSH SEAL R&D PROGRAMS 

(cont.) 


• TEXAS A&M: BRUSH SEAL ROTORDYNAMICS 

• Last stage of seal group develops higher pressure drop than 
previous stages 

• Increasing inlet tangential velocity increases leakage slightly 

• Cross-coupled stiffness coefficient very low and generally 
negative (stabilizing) 

• Whirl frequency ratio indicates brush seal Is extremely stable 

• Rotordynamic coefficients independent of seal spacing and inlet 
tangential velocity 

• Comparison with 8-cavity lab seal Indicate brush seal will 
generally improve rotordynamic characteristics 

• Report#: WL-TR-91-201 3 (original) 

corrected report in progress 

> 
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f V 

CURRENT BRUSH SEAL R&D PROGRAMS 


• EG&G SEALOL 

• Brush Seal Development Program 

• Advanced Brush Seal Development Program (New Start) 

• PRATT & WHITNEY 

• High Speed Brush Seal Development Program (New Start) 

• IN-HOUSE 

• Brush Seal Leakage Flow Modeling 

• Brush Seal Compressor Shroud Test 

< 


> 

BRUSH SEAL DEVELOPMENT PROGRAM 

EG&G SEALOL 


OBJECTIVE : 

Develop a comprehensive design methodology for application of brush 
seals to man-rated engines 

APPROACH : 

Conduct parametric testing of design variables for brush seals to 
define performance, fatigue, oxidation, and wear characteristics 

ACCOMPLISHMENTS : 

• Examined effect of bristle angle, bristle length, stiffness, staging, 
and packwldth on hysteresis, AP capability, leakage, wear 

• Significant hysteresis with speed, but not pressure 

• Increased Packwldth gave higher AP capability, lower leakage 

• Multiple Stage Brush Seal Performance 

- Leakage reduction 

— Unequal pressure distribution 

- May be controlled by using mixed stiffness designs 

- Staging with higher packwldth seals most effective 

• Trlbopalr Test and Evaluation 

• Evaluating Advanced Designs to Reduce Hysteresis 

• Designing Full-Scale Seal for Demonstrator Engine Test 

PROJECT ENGINEER: Capt Connie Dowler, 513-255-8210 
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ADVANCED BRUSH SEAL DEVELOPMENT PROGRAM 

EG&G SEALOL 


OBJECTIVE: 

Develop a comprehensive design methodology for application of 
advanced, high performance brush seals In man-rated engines 

APPROACH : 

Conduct experimental characterization of seal design and materials 
pairs and CFD modeling to maximize single- and multi-stage brush 
seal A P capability and axial and radial excursion accommodation. 

Investigate feasibility of non-contacting brush seal. 

ACCOMPLISHMENTS : 

New Start 

PROJECT ENGINEER: Lt Carolyn Sunderland, 513-255-8210 

V . - 4 


HIGH SPEED BRUSH SEAL DEVELOPMENT PROGRAM 

Pratt & Whitney 


OBJECTIVE: 

Provide verification of advanced brush seal technology needed to 
enable continued application of brush seals In IHPTET Phase (( 
demonstrator engines. 

APPROACH: 

Conduct an application study using IHPTET Phase II engine design 
and mission flight cycle to determine surface speed and temperature 
requirements for brush seals. Design and fabricate brush seals 
for rig testing to verify their capability to operate at IHPTET 
Phase II conditions (~ 1400fps, 1400 F). 

ACCOMPLISHMENTS: 

New Start 

PROJECT ENGINEER: Lt Carolyn Sunderland, 513-255-8210 
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BRUSH SEAL LEAKAGE FLOW MODEL 

FLOW MODEL APPROACH 


• 1-D MODEL, WITH RADIAL AND AXIAL FLOW CONSIDERED 

• FLOW CORRELATIONS FROM KNUDSEN AND KATZ 

- Laminar Flow 

- Turbulent Flow 

- Transition Flow 

• SINGLE CORRELATING PARAMETER: EFFECTIVE THICKNESS 

• MODEL 1 

- Linear 

- Square Array Bristle Bed 

• MODEL 2 

- Linear 

- Hexagonal Pack Bristle Bed 

• MODEL 3 

- Modification of Model 2 to Account for Curvature Effects 

“ Effective Thickness Constant Between Journal and Mean Diameter 

- Bristle Bed Configuration Varies 

• Increased Transverse Bristle Spacing 

• Increased Leakage Flow Area at Mean Diameter 

• Hexagonal Pack Bristle Bed 


TYPICAL MODEL PERFORMANCE 
CHARACTERISTICS 



• SHAPE OF 0 VERSUS P r IS EXPECTED CHARACTERISTIC CURVE 

• FLOW INCREASES WITH B AS THE BRISTLE PACK OPENS UP 

• TRANSITION FROM LAMINAR TO TURBULENT FLOW NEAR P r = 2 
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LEAKAGE PERFORMANCE DATA 



• HYSTERESIS EFFECT REMOVED 

• ALLISON DATA ARE AVERAGED RESULTS FOR SEVERAL SEALS 

• TELEDYNE CAE DATA ARE FOR A BASELINE BRUSH SEAL 

• CROSS DATA PROVIDED CHECK ON EFFECT OF BRISTLE DENSITY 

• NASA LEWIS DATA PROVIDED CHECK ON SEAL DIAMETER EFFECT 


BRUSH SEAL LEAKAGE FLOW MODEL 

MODEL RESULTS 




• WIDE RANGE OF B INDICATIVE OF DIFFERENCE IN ACTUAL SEAL THICKNESS 

• B/Bmin INDICATOR OF BRUSH SEAL LEAKAGE EFFECTIVENESS 

• EG&G DOUBLE PACKWIDTH SEAL LESS EFFICIENT THAN STANDARD DENSITY SEAL 

• CROSS SEAL VERY EFFICIENT 

• NASA SEAL LESS EFFICIENT THAN LARGER DIAMETER SEALS 
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BRUSH SEAL COMPRESSOR SHROUD TEST 

• PART OF ADLARF TEST PROGRAM IN COMPRESSOR RESEARCH FACILITY 

• TENTATIVELY SCHEDULED TO TEST JAN 93 

• INVESTIGATING CAPABILITY OF BRUSH SEAL SHROUD TO IMPROVE 

- BLADE VIBRATIONS (DAMPING) 

- STALL MARGIN 

- EFFICIENCY 

• PROJECT ENGINEER: Lt Carolyn Sunderland, 513-255-8210 


DEMONSTRATOR ENGINE BRUSH SEAL TESTING 


• TESTING TO DATE: 

- PRIMARILY IN TURBINE SECTION 

- HAYNES 25 BRISTLES 

- CHROME CARBIDE OR ALUMINUM OXIDE COATING 

- 20 TO 80% REDUCTION IN LEAKAGE OVER LAB SEALS 

- REDUCED HEAT GENERATION 

- GENERALLY GOOD DURABILITY 

- MAX CONDITIONS: 1275 FPS, 1130 F, 55 PSI AP 

• FUTURE TESTING: 

- PLANNED FOR ALL DEMO ENGINES (ATEGG, JTDE, JTAGG, ETEC) 

- COMPRESSOR AND IGV LOCATIONS 

- BRUSH SEAL SHROUD 

- HIGHER SURFACE SPEED AND TEMPERATURE (>1400 F, >1400 fps) 
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F1 19 BRUSH SEALS 


• BILL OF MATERIALS 

- 3 STATIC SEALS 

- HPT & LPT SHROUD 


• TECHNOLOGY TRANSITION PLANS 

- HPT (3 Locations) 

- LPT (2 Locations) 

- REPLACE LAB SEALS WITH BRUSH SEALS AT ALL 
COMPRESSOR INTERSTAGE LOCATIONS 

- 0.9% TFSC IMPROVEMENT (TOTAL) 


V 




SUMMARY 

s 


• AF COMMITTED TO DEVELOPING AND TRANSITIONING BRUSH 
SEAL TECHNOLOGY 

• EMPHASIZING BRUSH SEAL RIG TESTING (CONTRACT AND 
IN-HOUSE) 

• PUSHING FOR INCORPORATION OF BRUSH SEALS IN ALL 
DEMONSTRATOR ENGINES 

• WORKING WITH PROGRAM OFFICES TO TRANSITION BRUSH SEAL 
TECHNOLOGY TO OPERATIONAL AND DEVELOPMENT ENGINES 


) 
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Navy GTE Seal Development Activity 
by 

Carl P. Grala 

Naval Air Warfare Center 

Aircraft Division- Trenton 


Under the auspices of the Integrated High Performance Turbine 
Engine Technology Initiative, the Naval Air Warfare Center conducts 
advanced development programs for demonstration in the next generation 
of air-breathing propulsion systems. Among the target technologies 
are gas path and lube oiL seals. Two development efforts currently 
being managed by NAWCAD are the High Performance Compressor Discharge 
Film-Riding Face Seal and the Subsonic Core High Speed Air/Oil Seal. 

The High Performance Compressor Discharge Film-Riding Face Seal 
Program aims at reducing parasitic leakage through application of a 
film-riding face seal concept to the compressor discharge location of 
a Phase II IHPTET engine. An order-of-magnitude leakage reduction 
relative to current labyrinth seal configurations is expected. 
Performance goals for these seals are (i) 1200 F air temperature, (ii) 
800 feet-per-second surface velocity, and (iii) 600 PSI differential 
pressure. The two designs chosen for fabrication and rig test are a 
spiral groove and a Raleigh step seal. Rig testing is currently 
underway. 

The Subsonic Core High Speed Air/Oil Seal Program is developing 
shaft-to-ground seals for next-generation propulsion systems that will 
minimize leakage and provide full life. Significantly higher rotor 
speeds and temperatures will be experienced. Technologies being 
exploited include, hydrodynamic lift assist features, ultra light 
weight designs, and improved cooling schemes. Parametric testing has 
been completed, a final seal design is entering the endurance test 
phase. 
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MECHANICAL SYSTEMS TECHNOLOGY 
S.O.A. vs FUTURE REQUIREMENTS 



SPEED AND TEMPERATURE TRENDS HAVE AN 
ADVERSE EFFECT ON MECHANICAL SYSTEMS: 

- STRESS CYCLES ACCUMULATE FASTER, 
LIFE IN HOURS GOES DOWN. 

- INTERFACIAL HEAT GENERATION GOES 
UP WITH VELOCITY. 

- MATERIALS ALLOWABLES DROP WITH 
INCREASING TEMPERATURE. 

- REQUIRED WEIGHT REDUCTIONS MAGNIFY 
THE ENVIRONMENTAL CHALLENGE. 


MON 


BEARING SPEED 


BRG / LUBE TEMPERATURE 


fi.O 

4.0 

3.0 



aoo 

600 

DEG. F 

400 




K 


i?0 TO ?00 LBS 


2.0 

S.O.A. 1 997 2003 


200 

S.O.A. 1997 2003 



MECHANICAL SYSTEMS PROJECT 



TOKKTON. NSW JKRSCY 


CHALLENGE: 

* SIGNIFICANTLY EXPAND COMPONENT OPERABILITY 

- SPEED, TEMPERATURE, LOAD 

♦ SIGNIFICANT COMPONENT/SYSTEM WEIGHT REDUCTIONS 


♦ REACT TO CURRENT DESIGN PRACTICE DEFICIENCIES 
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MECHANICAL SYSTEMS PROJECT 

(WR22-P64) 



TREKTON. NSW JVtMFt 


TASK 1: RADIAL AND AXIAL BEARINGS 

TASK 2: LUBE OIL SEALS 

TASK 3: STATIC AND DYNAMIC GAS PATH SEALS 



MECHANICAL SYSTEMS PROJECT 

TASK 1: AXIAL AND RADIAL BEARINGS 



INCREASED ENDURANCE FOR EXTENDEO COVERAGE 
FOR ASW/ASUW/AEW/EW/C 3/DRUG MISSIONS 

> INCREASED STANDOFF RANGE FOR POWER 
PROJECTION ASHORE 

HIGH PERFORMANCE ENGINES FOR EMERGING 
SYSTEM REQUIREMENTS (AX, NATF. ASTOVL. SSF) 


IMPROVED LIFE CYCLE COST 
THROUGH IMPROVED RELIABILITY 

IMPROVED A/C AVAILABILITY 


IHPTET/S&T REUAN-CE-JPL 


ClQALS; 


MAGNETIC BEARINGS 
• THRUST COMPENSATION 

• RADIAL LOAD COMPENSATION 



INNOVATIVE DESIGN 
IMPROVED MATERIALS 
IMPROVED DAMPERS 



PH 1 

PH 2 

PH 3 


PH 1 

PH II 

PH III 

TEMPERATURE (F) 

400 

600 

800 

TIME ON STATION 

•1.8* 

♦3.0* 

♦ 6 * 

BEARING SPEED (MDN) 

2.5 

3.7 

3.0 

4.2 

3.5 

4.5 

RANGE 

• 2 * 

♦ 4 * 

♦ 7 * 

WEIGHT 

- 6 % 

- 10 % 

- 20 % 

PAYLOAD 

• 18 * 

• 36 % 

♦ 60 % 





FN / WT 

•5.2* 

♦11.1* 

•25* 
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MECHANICAL SYSTEMS PROJECT 

LUBE OIL SEALS 



■nttMTOH, MIW JCftSCY 


SUBSONIC CORE HIGH SPEED AIR/OIL SEAL 


CONTRACTOR : PRATT & WHITNEY 

COST : $ 430 K 


OBJECTIVE : DEVELOP SHAFT-TO-GROUND SUMP SEAL 

SYSTEMS FOR IHPTET PHASE II CONDITIONS 
THAT MINIMIZE LEAKAGE, GIVE FULL LIFE. 


GOALS : 


PHASE II CONDITIONS: 

♦ 600 FPS * 750 F AIR 

♦ 60 PSID ♦ 400 F OIL 

♦ FULL LIFE 
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MECHANICAL SYSTEMS PROJECT 

LUBE OIL SEALS 

APPROACH : * ANALYTICALLY ASSESS MULTIPLE SEALS 

* DETAIL DESIGN AND FABRICATE THE TWO 
BEST CANDIDATE SEALS 

* 25 HRS OPERABILITY EACH, REVISE 

* ENDURANCE TESTING 

TECHNOLOGIES : * HYDRODYNAMIC LIFT ASSIST (STEIN) 

* ULTRA-LIGHTWEIGHT DESIGN (REXNORD) 

* IMPROVED PACKAGING / COOLING 

ADVANCEMENT 

BEYOND SOA : * ORDER OF MAGNITUDE LEAKAGE REDUCTION 

RELATIVE TO LAB SEALS (.1% SFC PER) 

* 30% SPEED CAPABILITY INCREASE 


MECHANICAL SYSTEMS PROJECT 





STEIN HYDRODYNAMIC CIRCUMFERENTIAL 
NON-CONTACTING SEGMENTED SEAL 


PROS: 

- EXTENSION OF SUCCESSFUL 
SEGMENTED RING SEAL 

- LOW LEAKAGE 

- WEAR TO 6 MILS O.K. 

- LIGHTWEIGHT 

- WINDBACK ALLOWS MINIMUM 
LEAKAGE AND CONTAMINATION 


CONS: 

- MINIMAL EXPERIENCE WITH 
HYDRODYNAMIC LIFT GEOM. 
IN CARBON BORE. 

- THERMAL CONEING AND 
MISALIGNMENT CONCERNS 




m ' rnLi -i— I T1 

VW M 

JO, 
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REXNORD CARTRIDGE-TYPE CONTACTING FACE SEAL 


PROS: 

- VERY LIGHTWEIGHT LOW 
DRAG DESIGN IMPROVES 
SPEED CAPABILITY. 

- TOLERANT TO CONEING 
AND MISALIGNMENT. 

- LOW TO MOD. LEAKAGE 

- LOW OPERATING LOADS 
IMPROVE LIFE 


CONS: 

- NUMEROUS PARTS 

- NARROW CARBON NOSE 
PRESENTS HANDLING 
DAMAGE RISK 

- POSSIBLE LEAKAGE AT 
VERY LOW DELTA P's 




MECHANICAL SYSTEMS PROJECT 

LUBE OIL SEALS 



STATUS : * PHASE I OPERABILITY EVAL COMPLETE 

- MIXED RESULTS 

• STEIN HYDRODYNAMIC CIRCUMFERENTIAL SEAL: 

• STATIC CAL DONE, VERY LOW LEAKAGE 

• RAN SUCCESSFULLY TO 600 FPS, 60 PSID ! 

• LEAKGE CONSISTENTLY LOW. THEN - 

* BROKE EXTENSION SPRING - EASY FIX IN HAND 

* REXNORD CARTIDGE-TYPE FACE SEAL: 

• STATIC CAL COMPLETED (INITIAL SEC. SEAL PROB) 

• VERY SUCCESSFUL THROUGH TWO DYNAMIC TESTS TO 
600 FPS. 60 PSID ! 

** LEAKAGE TOOK OFF, SEAL FAILED 
• INVESTIGATION IN PROGRESS. 

SUMMARY : * MINOR REVISIONS TO STEIN SEAL, HIGH 

CONFIDENCE FOR ENDURANCE PHASE 

• ASSESSMET OF REXNORD SEAL IN PROGRESS, 
ENDURANCE PROSPECTS TBD 

• HAVE SHOWN STABLE LOW LEAKAGE OPERATION 
AT AGGRESSIVE GOAL CONDITIONS FOR BOTH. 
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MECHANICAL SYSTEMS PROJECT 

TASK 3: STATIC & DYNAMIC GAS PATH SEALS 



FLEET ISSUES ADDRESSED; 

INCREASED ENDURANCE FOR EXTENDED COVERAGE 
FOR ASW/ASUW/AEW/EW/C 3/DRUG MISSIONS 


TURBINE GOALS. 


• INCREASED STANDOFF RANGE FOR POWER 
PROJECTION ASHORE 

• HIGH PERFORMANCE ENGINES FOR EMERGING 
SYSTEM REQUIREMENTS (AX, NATF. ASTOVL. SSF> 

. BRUSH / 

• IMPROVED LIFE CYCLE COST 
THROUGH IMPROVED RELIABILITY 

• IMPROVED A/C AVAILABILITY 

. 









\\\ /k 


QQAL& 

KYA / * 


• NO FORMAL IHPTET GOALS, 



CONTRIBUTES TO COMPRESSOR 4 




TECHNOLOGIES; 

* FILM-RIDING FACE SEALS: 

HYDRODYNAMIC. HYDROSTATIC 


ABRADABLE SEALS 



PH l 

PH 2 

PH 3 


PH 1 

PH II 

PH III 

TEMPERATURE (F) 

800 

900 

1200 

TIME ON STATION 

♦1.8% 

♦3.0% 

* 6 % 

SPEED (FPS) 

700 

860 

1200 

RANGE 

♦ 2 % 

♦ 4 % 

♦ 7 % 

WEIGHT 

- 5 % 

- 10 % 

- 20 % 

PAYLOAD 

♦ 18 % 

♦ 35 % 

♦ 60 % 





FN / WT 

♦5.2% 

♦11.1% 

♦25.% 



MECHANICAL SYSTEMS PROJECT 

GAS PATH SEALS 



TWHTO*. MIW JKASCY 


HIGH PERFORMANCE CD FILM-RIDING FACE SEAL 


CONTRACTOR : 
CONTRACT NO. : 
COST : 
OBJECTIVE : 


ALLISON 

N00 140-39-C-2728 
$460 K 

DEVELOP/DEMO FILM RIDING CD FACE 
SEAL° FOR VERY HIGH PRESSURE RATIO 
PHASE II ENGINES. 


GOALS : 


1200 F AIR 
800 FPS 
600 PSID 
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MECHANICAL SYSTEMS PROJECT 

GAS PATH SEALS 



APPROACH : APPLY FILM-RIDING FACE SEAL CONCEPT TO 

PM II C.D. APPLICATION, USING IMPROVED 
DESIGN ANALYTICS, IMPROVED MATERIALS 
ASSESS MULTIPLE LIFT FEATURES, 

FABRICATE & TEST 


TECHNOLOGIES : ‘TRANSIENT DYNAMIC FILM ANALYSIS 

•SPIRAL GROOVE & RAYLEIGH PAD 
LIFT GEOMETRIES 
•SILICON CARBIDE PRIMARY RING 
•IMPROVED PRESSURE BALANCE 


ADVANCED 

BEYOND SOA : REPLACES MULTIPLE LABYRINTH STAGES 

AT OVER AN ORDER OF MAGNITUDE 
LESS LEAKAGE 



MECHANICAL SYSTEMS PROJECT 



TWtKTO*, N*W JW1IY 


Primary 
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MECHANICAL SYSTEMS PROJECT 
GAS PATH SEALS 


TRANSITION/ 

APPLICATION : * PHASE II JTAGG 

* APPLICABLE TO ANY HIGH PRESSURE 
RATIO MACHINE 

PAYOFFS : * SIGNIFICANT CYCLE EFFICIENCY BENEFITS 

♦ 0.5 % THROUGH REDUCED LEAKAGE 





MECHNICAL SYSTEMS PROJECT 

GAS PATH SEALS 



PROGRESS/STATUS : 


•FABRICATION OF BOTH SEALS COMPLETE 


•STATIC CAL. TESTING OF BOTH SEALS COMPLETE 
•VERY LOW LEAKAGE 


•DURING EVALUATION OF CRANE SEAL - RIG INDUCED 
RUB OCCURRED - SEAL OK; REWORKABLE 


•KAYDON SECONDARY SEAL FAILURE OCCURRED 
AWAITING REVISED HARDWARE 
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GAS SF.AL CODF. DEVELOPMENT AND ANALYSIS 


P. Basu 
EG&G 

1. Developed opening force versus film thickness curves (F-h curves) at ID (hi) at different speeds for a given 
set of OD and ID pressures. The constant closing force line intersects each one of the above curves at the 
corresponding operating point. The slope of a F-h curve at the operating point gives that average Film 
stiffness. (Fig. 1) 

2. Developed film thickness profile along the interface at the operating point for a given combination of speed 
and pressure. This curve combined with the corresponding stiffness value would enable the designer to 
know whether the seal is operating with a safe and stable film at that speed and pressure. (Fig. 2) 

3. Developed the operating film thickness (hi) versus system pressures for different speeds. This family of 
curves presents the operating film level over the entire envelope of operation in one picture. Hence, an 
unsafe operation at any off-design condition will become evident. (Fig. 3) 

4. Developed leakage versus system pressure curves for different speeds. This family of curves is the most 
useful one for the purpose of correlating the experimental data with the code prediction. (Fig. 4) 

5. Additionally, the code provides the pressure and Mach number distributions along the interface that are 
printed on a file for each operating point. 

6. Seal distortion 

Since the seal face distortions are of the same order of magnitude as the average operating film thickness, 
the code iteratively determines the distortions using a finite element module. (Ref. 1) 

7. Design 

A figner spring installed in a groove beneath the SiC inner diameter allowed for differential thermal 
expansion between retainer and seal ring while keeping the ring centered. Round smooth wire bent into 
shape to give centering force between rotating seal ring and the retainer. Also used to mount retainer on the 
shaft to accommodate centrifuged growth. 

Antimony impregnated carbon graphite stator ring is used. This provides high compressive strength and 
stability: 

The seals are designed to run with a very thin operating film (3-4 |J.m) to minimize leakage and maximize 
film stiffness. 

For more information 

1. Gardner, J„ Basu, P., Muraki, R. (1992), Analysis, Design, Laboratory Testing of a High Pressure, High 
Speed, Reduced Axial Length Gas Seal, Proceedings of the 13th International Conference on Fluid Sealing. 
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Distance from ID, mm Opening Force, kN 

19.2 1.3 2.2 3.1 4.0 4.9 5.8 6.7 


Inlet Pressure = 4.5 bar 

Outlet Pressure = 1.01 bar 

Fluid Temperature = 22 °C 


□ 

o 

A 


LEGEND 
0 RPM 
5000 RPM 
10000 RPM 
16000 RPM 
Closing .Force 



0-0 ^ 2^5 ^ 7.5 " 

Film Thickness at ID, /Ltm 


10.0 


12.5 


Figure 1. ‘F-h’ Curves 



Figure 2. Film Thickness Profile 
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Leak Rate, kg/h 


140 



Operating Film Thickness, nm 1 


Figure 3. System Pressure versus Operating Film Thickness Plots 


3.2 JL 


2.4-- 


LECEND 

A 5000RPM, TEST 
a 5000RPM, CALC. 
0 7500RPM. TEST 
+ 7500RPM, CALC. 
* 12500RPM, TEST 
+ 12S00RPM, CALC. 



35 49 63 77 

System Pressure, bar 


Figure 4. Leakage Correlation 
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Compliant Seal Development 


From notes and discussion of talk J# J. Gardner, EG&G- 
by R.C. Hendricks 

The compliant metallic seal combines the noncontact feature of the labyrinth seal 
the low leakage of a mechanical seal, and the compliant nature of the brush seal! 
It consists of several thin metallic elements or leaves mounted within a ring 
which is press fit into the housing, and in form, sort of resembles a lip seal 
sections wiping the shaft. A second set of overlapping cover leaves are placed 
on top of the shaft riding leaves which reduces leakage and provides stiffness. 
The leaves can be straight or angle cut. The shaft riding fingers are designed 
with mismatched curvature to provide lift off similar to the Rayleigh lift pads 
in mechanical seals with leading edge clearances nearly twice those of the 
trailing edge as shown by Fleming to be optimal for gas flows in convergent seal 
passages. Leading edge clearances range from 300 to 500 microinches. 

Balance pockets beneath the leaves provide fluid film feed to the "Rayleigh lift" 
surface and the proper balance ratio (mechanical seal) when combined with the 
static pressure and film pressure. The leaves flex in the radial direction and 
accomodate thermomechanical behavior as well as axial motion and angular 
mi alignment. 

In the static mode, there is a net closing force on the leaves. 

The seals were tested to 70 psi at speeds to 16,000 rpm or surface speeds to 330 
fps and temperatures from ambient to 440F. A slow cycle through the rig critical 
at 10,000 rpm induced a radial vibration response of 0.004 to 0.005 inch were 
accomodated by the seal. 

Preliminary performance data are encouraging demonstrating hydrodynamic liftoff 
and non contacting operation at pressure and speeds typical of gas turbine en- 
gines. The leakage performance data are significantly better than commerical 
labyrinth and brush seals which should be expected as this design incorporates 
the features of the low leakage face or mechanical seal along with the flexibility 
of the brush configuration. 

For more information 

1. Gardner, J., Basu, P„. , and Datta, A.: A New Compliant Seal Concept for Aero- 
space Applications, Fourth Int. Symp.. on Transport Phenomena and Dynamics of Ro- 
tating Machinery (ISROMAC-4), Vol . A, April 5-8, 1992, Honolulu, Hawaii, Sponsored 
by Pacific Center of Thermal-Fluid Engineering. 
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FI&UfiE 1 



Typical Pressure Distribution for Straight Cut Compliant SeaJ 
without Slot Cover 


P 

discharge 



dotted lines represent lines of constant 
pressure between the shaft the underside 
of the compliant metallic members 

solid lines represent the direction of the 
flow Held 
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Angle Cut Compliant Seal with Slot Cover 
Figure 3 


Straight finger Compliant Seal with Rayliegli Pads and Slot Cover 


End View 


Sectional view 



Figure 4 
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ACTIVELY CONTROLLED SHAFT SEALS 
FOR AEROSPACE APPLICATIONS 


Richard F. Salant and Paul Wolff 
School of Mechanical Engineering 
Georgia Institute of Technology 
Atlanta, Georgia 30332 


NASA Research Grant NAG 3-974 
NASA Technical Officer: M. P. Proctor 



Helium Purge Assembly 
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Face Holder 


Helium 


Housing 



Lox Side 


Hot Gas Side 


Carbon Face 


Rotating Face 


CONVENTIONAL MECHANICAL SEAL 



pressure 

force 


spring 

force 


FLOATING FACE 


pressure 


FIXED FACE 



High Pressure 



Leakage 
Low Pressure 
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NORMALIZED 

PRESSURE 



To Control Sys 

From Control 
Lox Side 



To Control Sys 

From Control Sys 
Hot Gas Side 


0 _ 

Thermocouple 


Electrode ^ Cartxm Face 


CONTROLLABLE SEAL 



CONED DEFORMABLE FACE ASSEMBLY 
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Leakage (ml/min) Film Thickness(lim) 


Piezoelectric Element 



Insulating Holder 


Thermocouple 


Carbon 

Face 


Insulating Bushing Electrode 


FACE HOLDER AND DEFORMABLE FACE ASSEMBLY 



Off Scale 



Voltoge (V) 


Leakage versus Voltoge - Static Test 
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Off Sco I© 



Film Thickness versus Voltage — Static Test 


3500 r 

3000 ■ 

^2500 • 
E 

^2000 - 
|l500 - 

JSC 

"1000 - 

500 • 

0 - 

0 



— i — i — . — * — i — i — i — . — i — . — i — ■ 

500 1000 1500 2000 

Voltoge (V) 

Leakage versus Voltage - Rotating Test 
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EXPERIMENTAL AND ANALYTICAL INVESTIGATION OF 
BRUSH SEALS 


M.J. BRAUN 
Professor 


V.V. KUDRIAVTSEV 
Visiting Scientist 


NASA grant NAG 3-969 


Project Manager: Mr. George Bobula 
Vehicle Propulsion Directorate 
US Army, MS 77-12 


OBJECTIVES OF THE PRESENTATION: 


• Report on the status of the project 


• Discussion that will broaden your research and 
and help to interpret results 
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OBJECTIVES AND ACCOMPLISHMENTS 


• Develop verified family of CFD codes for Analyzing Brush Seals 
-idealized(uncompliant) 2D configuration \/ 


-compliant 2D geometry not 

-compliant 3D no t 

• Experimental Facilities for the adequate code verification 

-stationary bristlesf cylinders) \/ 

-moving bristles not 


additional experiments were designed specifically 
for code validation and investigation of features 
to be incorporated in the computer program 

• Qualitative and quantitative analyses of the Fluid Flow in 
the Brush Seal Configuration 

-flow around one bristle, level 1 \/ 

-flow around several bristles, level 2 \/ 

-flow in the deep tube bundles, level 3 \/ 

(intermediate pitch-to-diameter ratio) 


-flow through uncompliant brush 
prototypes, level 4 \/ 

(small pitch-to-diameter ratio) 


PRESENT CODE CHARACTERISTICS 


• finite-difference primary variables (u,v,p) formulation 

• collocated grid 

• conservative formulation of the momentum equations 

• conservative formulation of the pressure equation 

• conservative fully populated of the boundary conditions 

• steady-state solutions 

• transient solutions 

• 2nd and 3rd order convective schemes are implemented 

• implicit under/over relaxation (if necessary, optional) 

• accommodates an arbitrary geometrical configuration in 
Cartesian coordinates(non-contiguous boundaries) 

• regular and nonregular grid 

• graphical post-processor capabilities with PV-Wave (Unix) 
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MODEL ASSUMPTIONS: 


• incompressible fluid, constant properties 

• isothermal model 

• leakage flow 

- M < 0.3 

- Re (low, intermediate) 

• no shaft rotation 

• stationary bristles 


DETERMINED FLOW RESISTANCES THROUGH THE 
BRUSH SEAL 
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3D VIEW OF THE LINEAR BRUSH IN THE CHANNEL 



SCOPE OF THE WORK 


• development and evaluation of a numerical procedure aimed at 
the solution of the incompressible flow problem through a densely 
packed array of cylindrical pins at low to medium Reynolds numbers 

[The method promotes the solution of the 2-D Navier-Stokes 
equations written in a conservative form in primitive 
variables(u,v,p). The numerical implementation follows an implicit 
formulation that uses the alternate direction integration(ADI) 
method applied to a collocated grid. For the pressure solution, the 
two momentum equations are joined by an elliptic Poisson pressure 
equation(Roache, 1981, Ghia et al., 1983). ] 


The algorithm is applied to : 

-the classical analysis of flow and pressure development around a 
cylinder 

- an array of cylinders located in a square channel 


• flow and pressure patterns are discussed parametrically with 
respect to the variation of the Reynolds number. 
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GOVERNING EQUATIONS 


Sf 


, *(UU) . a<UV) 

1 aX + ~1T~ 




aV 4W) «(UV) 


«P , 1 «2v .Aa 



Taking the divergence of the X- and Y-momentum equations, and 
using continuity for simplification, one obtains a pressure equation 
that appears under a Poisson type format(Roache, 1981). 



D is the dilation term, 


D= 


_ dU , aV 




During calculations, D represents the residual that has to shrink to 
zero if continuity is to be satisfied. 


Boundary and Initial Conditions for the Momentum Equation. 


★The boundary conditions on the lateral solid walls assume non-slip 
conditions and non-porous walls. 

★The inflow velocity is uniform with U=l, and V=0(angle of 
incidence is zero). A reference pressure is assigned to one point on 
the inflow boundary. 


Special Form of Open Boundary Condition(Exit Flow) 


(i) the satisfaction of the continuity equation(specifically in the 
direction of the U velocity), 


dU _ aV 

SY 


(ii) the condition for fully developed velocity gradient for the 
vertical component 



(iii) satisfaction of the momentum equation(in the direction of U 
velocity) 


a(UU) a(UV) 

iy “ - 1 + ~;v 
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Boundary Conditions for the pressure Poisson Equation 


★ The dynamic pressures P, are determined from the balance of 
the normal forces with the inertia and viscous forces. The effects of 

the terms ^y and ^y both at inflow and outflow boundaries have 

been considered negligible. The following formulation is totally 
independent of the internal boundary configuration, and thus 
applicable with no restrictions. 


On the non-contiguous internal boundaries of the pins, in addition 
to previous equation, one needs to add an expression for 

1 + 1 + Re + 

5a 

At the channel’s upper and lower walls the pressure boundary 
condition takes a simplified form due to the fact that U=V=0 



Initial condition s(T=0). The input velocities are given as U=l, 
V=Q. The pressures are set initially to an arbitrary and operator 
chosen constant (P=P ref =const). 

NUMERICAL IMPLEMENTATION : ADI METHOD 


5 * 2 Vx 


• The discretization of the system of governing equations 
introduced above, follows the use of the ADI applied to a 
collocated grid. The procedure uses the full direct approximation of 
each term within the differential equation on every half time step, 
W 2. One obtains the following system of linear algebraic 

equations. 


V ‘%!/2 - =& A ~M° +1/2 + A yy( n ) n 1+ C»' +1/2 + Cy(u)“ + Q x (p)» 

' ^‘^/r 172 = +A„(u)»+1] + +C y (u)"+ 1 +Q x (p)" 

~ ~^72 y ‘ + Cx(") n+l/J + W*+ Oy(P)* 

=ife A «W" +I/J + A yy( T ) ,,+1 l+Cx(v) ,,+,/J + W +, +Qy(p) n 
-■+l/2»n+ 1 «,n+l . . 

E — w=* — a^p r +A yy (pr +1 + q(u,v)“+> 


— *+l,D+l »+l/2^»+l 

2 — / = ^xx{?)^ /2,n+l +Ayy(pr f ‘•■+ 1 + Q(u,v ) n+1 


The spatial derivatives, with the exception of the convection terms 
and cross-derivatives, are approximated by a second order central 
finite-difference. For the convection terms the first-order 
conservative scheme proposed by Torrance(1968) is employed. 
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Solution Procedure 


The solution follows according to the time-dependent alternative- 
direction method using a tridiagonal matrix elimination. The steps 
of the solution are as follows: 

• solve in the x direction for U velocity at the n+1/2 time 
step 

• solve in the y direction for the U velocity at the n-fl 
time step 

• solve in the x direction for V velocity at the n+1/2 time 
step 

• solve in the y direction for the V velocity at the n+1 
time step 

• solve the Poisson’s pressure equation at the n+1 time 
step by means of the pseudo-transient method within the 
set of internal ADI iterations n^=l,...n^ : 

(i) in the x direction at the s+1/2 pseudo-time step 

(ii) in the y direction at the s+1 pseudo-time step 
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ACCURACY CONSIDERATIONS 


• Qualitative comparison of calculated flowfield with experimental 
results of Bouard and Coutanceau(JFM,v.l01, 1980) 




Pressure Drops, AP=P |eading edge - P traj|jng edge 


Re d 

this study 

Fornberg, 1980 

100 

0.689 

0.70 

200 

0.670 

0.63 

300 

0.640 

0.60 
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FLOW 



A 

L_±d 

l * { "UPPER* WALL 

L 

V 


-1 

_D "“TV* 

-0*^" ""rv 

O 

0 

•BOTTOM* WAIL 

1 row 
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FIG. I CHARACTERISTIC CONFIGURATION 
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Effect of Re number 



FLOW THROUGH BASIC ARRAY CONFIGURATION 7 ROWS X 11 PENS 
AIRjC^ — 100 ; H) Re, =2000: 

CJDETAEL OF FLOW IN THE QUADRANT INDICATED IN FIG. 2A 
D)DETAEL OF FLOW IN THE QUADRANT INDICATED IN FIG. 2B 
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NON-DIMENSIONAL PRESSURE NON-DIMENSIONAL PRESSURE 


PRESSURE MAPS INSIDE TUBES ARRAYS 


LEADING EDGE 



TRAILING EDGE 




^22^S^, NSIONAL PRESSURE VARIATION IN THE TRANSVERSE 
Di-KjbUTION 

A ) - £1 £^*££^55? 0F THE R0W - AT R»h=100 AND Re. =2000 

B) , D) TRAILING EDGE OF THE ROW, AT Re^lOO AND R^ h =2000 


Deep bank(72) cylinders, Level 3 
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NON DIMENSIONAL PRESSURE 



CHANNEL BOTTOM 


NON-DIMENSIONAL PRESSURE MAPS FOR PRESSURE DROPS IN THE 
BASIC ARRAY CONFIGURATION (7 ROWS X 11 PINS) 

A )R^ =100; B)Rc h =2000; 


Cascade of Two Elements ™ IC , recirculation 



ZONE A ZONE B 


i-— □ C I 



ZONE C 


FLOW ACROSS ARRAY OF PINS AT RE, =2000, PTDR t =1, PTDR~=1 

A) SIX ROWS ^ L 'T 

B) CASCADE OF TWO ELEMENTS(THREE ROWS EACH), C=3d 

C) CASCADE OF TWO ELEMENTS(THREE ROWS EACH), C=5d 
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VAR IATION OF THE FLOW FIELD STRUCTURE WITH THE CHANGE IN 
THE LONGITUDINAL PITCH 




] 
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PRELIMINARY CODE VERIFICATION RESULTS 


3 rows of pins 


Re P (psi) 

P(Pa) 

,u 2 

AP(nond) 

91.3 2.263 

15,602 

7.09 

2200. 

195exp. 4.02 

27,758 

32.16 

863.0 

num. 4.31 

(7% difference) 


926.0 

327exp. 6.82 

47,056 

91.29 

516 

num. 8.6 

(26% difference) 


650 


Oil temperature: T=70 F 

Oil density=950 kg/m 3 
Fluid velocities: 8.6 cm/sek 
18.5 cm/sek 
31.0 cm/sek 


Brush Design Applications 


Pressure Drop Estimation For the Characteristic Brush Section 


inflow 

u=l 



R^=20 

outflow 


Converged Solution, n |t =5000 



History of the Pressure Drop Convergence, Re^=20 
PTDR Z = PTDRj^m 
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Brush Design Applications 


Alternative Configuration: Pin AiT&y in the Channel 


inflow 

u=I 


outflow 



Re rf =20 

PTDR l = PTDJtj.ssll.8 


CONCLUSIONS: 


• approximate mathematical model with distributed parameters 
and without porous media assumption 

• developed and evaluated a new computational algorithm for 
the solution of N-S equations in (u,v,p) formulation 

• systematic analyses of the fluid flow in the brush sea 
components 

• typical pressure and velocity distributions were estimated 

• nonlinear behaviour of the pressure drop vs PTDR^ 

• Pressure paradox for highly recirculative flows 

• Capabilities of calculation of the pressure drop 

for a given bristle geometry and a typical brush segment 

• PTDR is extemely influential in the nature of the overall flow 
inside the array and the subsequent pressure drops 

• Pressure distribution inside the brush is non-symetrical and varies 
across the cross section 

• [> e K *Uf or JJ iat *9 n on flow structure within tube bundles. 
Butterfly flow formation. 

Near wall jets that are instrumental in the butterfly formation. 

• transients are negligible, Re<1000 
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DYNAMICS OF FACE SEALS FOR HIGH SPEED TURBOMACHINERY 

Simon Le efe 
Project Engineer 
BHR Group Ltd. 

Cranf ield, Bedford, MK43 oaj, England 


ABS TR ACT 


Face seals In rocket engine fuel and oxldlser turbopumps have been the 
subject of Intense investigation for over 25 years. Whilst advances have 
been made in the understanding of thin film lubrication between seal faces; 
valuable data has been produced on the friction and wear of material pairs 
in cryogenic environments; pioneering work has been done on the effect of 
lubricant phase change In seals; and many Improvements have been made In 
mechanical seal design, relatively superficial attention has been given to 
the vibrational dynamics of face seals in high-speed turbomachinery. 

BHR Group Ltd. (formerly BHRA) has recently completed the first stage of a 
study, commissioned by the European Space Agency, to investigate this area. 
This has involved the development of a two-dimensional adiabatic, turbulent 
lubrication model for thick gas film applications, the production of an 
integrated mathematical model of gas seal vibrational dynamics for thin 
film applications, implementation in software, the undertaking of an 
experimental programme to validate software against variations in 
operating conditions and design variables, and suggestions for improved seal 
design. 
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BHR Group Limited 



CHIEF EXECUTIVE 
Ian Cooper 



SUBSIDIARY OPERATIONS 



CAUec Limited - Oil & gas 

NIMIX - NorHntu«v® mixing equipment 

Fluid Developments limited - Abrasive jet cutting equipment 

SFK Technology limited - Software development 


CUR REN T PROJECTS 


• Applications 

rotary 

reciprocating 

static 

• Technology 

experimental 
analytical 
design studies 

• Organisation 

direct contract 

consortium + Dept, of Trade and Industry 
consortium + EEC 
independent consortium 
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FLUID SEALING TECHNOLOGY 


Independent facilities and expertise:- 

• Seal Analysis:- 

• Seal testing:- 

• Pump Loops:- 

• Site Measurements:- 

• Design Audit: - 

• Rig Manufacturer- 



thermal, mechanical, lubrication 

oil, water, gas, cryogens, contaminants 

oil, water, slurry, water/air 

fixed, portable 

analysis, critical review 

design, build, modification 






FLUID FILM 


STRUCTURAL 

DISTORTION 


SEAL RING 
DYNAMICS 



SEAL 

TESTS 


INTEGRATED 

COMPUTER 

MODEL 


Validated by 

J 


PARAMETRIC 

STUDY 

RECOMMENDATIONS 



199 


S H IL Q fcQ (li 


Leakage Buffer Gas 


H 2° + GH 2 

% 


LOX 



Bearings / Buffer Seals Turbine 


Seal 


FACE SEAL OPTIONS: 

• Plain 

• Hydrostatic 

• Self-acting 



LITERATURE SURVEY 
MAIN FINDINGS - SEALING PRACTICE 

50 


20 


10 


o 

\ 


Limit of demonstrated 
success of plain face 
seals In LOX (NASA 1979) 

/ 


\ 


50 100 

Speed, m/s 


150 


■ Sell -acting face-seals in LOX (NASA evaluation) 
o Plain face -seals in LOX (demonstrated success) 

+ Plain face seals in LH2 (Japanese test programme) 

Plain face seals in GN2 
(BHRG test programme) 
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MODELLING - FILM LUBRICATION ANALYSIS 


Laminar, isothermal - Reynolds equation 

• Liquid films - treatment of cavitation 

• Gas films - grid design (adaptive, graded, etc.) 

- algorithm design (implicit, alternating, 
multigrid, ‘interior co-location', etc.) 


Turbulent lubrication - Hirs* bulk flow 

- Constantinescu 

- Ng & Pan 


Non -isothermal (higher Mach number) compressible flow 

• 1-D (radial) adiabatic model with radial taper 
and entrance effects - Zuk 


Two-phase (boiling interface) films 

• 1-D models 

• Stability approached from consideration of equilibrium 
film thickness vs. opening force curves 

(i.e not from dynamic analysis) 



LITERATURE SURVEY 
MODELLING - DYNAMICS 


• IKI and fCj matrices from fluid film analysis 
then dynamics as a separate problem 

• Integrated analysis - fluid film forces and moments 
calculated at each timestep 

• Excitation mechanisms 

• Number of vibrational degrees -of -freedom 

• Thermal and vibrational transients 3 or 4 orders of 
magnitude difference in timescale - separate problem 


@[p [fca©©lL 

• Concentrate initially on 
DYNAMICS 

• Gas seal assumed 

(2-D transient 2-phase prohibitive 
within commercial constraints) 

• Transient structural distortions 

• Full transient lubrication analysis 

• Turbulent, adiabatic AND laminar, 
isothermal leakage flow 

• Choked exit conditions catered for 

• 4 vibrational degrees of freedom 

• Mechanical damping 



MATHEMATICAL MODEL 
TURBULENT. ADIABATIC FLOW 


Iterate round 
instantaneous 
equations 


Radial and circumferential velocities: 

«, - - It) w *' ere »«■ ' 

and 2 (optional 

u ■ - i£. ♦ OS. inertia) 

* n r»# 2 

G r and G m from Hirs‘ bulk flow turbulent 
lubrication theory (or - 1/12, laminar) 

Given in terms of Reynolds numbers "as seen by" 
rotor and stator (different) 

Shear stresses from these Reynolds numbers 

2-D adiabatic energy equation relates pressure 
to circumferential shear stress for density at 
current tlmestep 


Timestepping 


Substitute u r and u m in continuity equation as 
knowns* and find density at next tlmestep 


Use film thickness, h. at mid-timestep throughout procedure 
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MATHEMATICAL MODEL 

LAMINAR, ISOTHERMAL FLOW 


Compressible Reynolds equation with ideal gas assumption 



Time discretisation: 



No energy equation required 

Velocities and shear stresses from pressure gradients 


• Check exit Mach number distribution for condition 
(isothermal flor) 


M < 


/? 


Check Reynolds number distribution for condition 

(laminar flor) 


Re * Re 
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MATHEMATlCAL-MQDEIJLINg 

FACE CONING 


Sources: 

• Pre -lapped taper and clamping forces 


Bellows load 





• Rotor centrifugal inertia 




• Change in ambient temperature for seal ring assemblies 



• Sealed pressure differential 


• Interfacial heat generation 



• Interfacial pressure distribution 




MATHEMATICAL MO DELLING 

INTERFACIAL PRESSURE CONING 



— 
3. Use influence coefficients to 
calculate face coning 


S 
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©IP 

IF® (Mil® ®l®fl 


< 4 > 




MAT HE MATI C AL i>lQi) a LIN£ 

DYNAMICS 


Equations of motion time-discretised by Newmark's method: 

l “ r ‘ ,F, .. 1 ' £ (I >.) - ~ «>. 

"~r 

«u 

|F)..,= 1^1(11.., + ft^U 

\ 

linear stiffness all nonlinear forces 

(bellons) (fluid film and 


mechanical dampers) 


Collect terms In (X) and solve: 


U U= «*•'*« (ii.) ♦ (X).) 


Inertia and stiffness matrices, [M] and IkJ, are diagonal 
in the absence of lateral stiffness, so that inverses are trivial 
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MATHEMATICAL MODELLING 
TIMESTEPPING 


COMPUTATION 

USES 

TO PREDICT 

Fluid film shape 

V* V* 

< x >, <*). f. 


Fluid film pressure 
distribution, p 

p * \*v 2 


Velocity and shear 
stress distributions 

p... ^ 

,r »*i P-..J 

Coning 



Fluid film forces 
and moments 

Mechanical damper 
forces and moments 

p - v. 

(i Taylor 

txpansioi 

^ F n«id ^ 41 

file 

<V»„, — 


Dynamics 

»>. <*». 



□ 


Output of fundamental 
importance 



MATHEMATICAL MODEM IMG 
APPROACH TO THERMAL TRANSIENT CONING 

PROBLEM 



Thermal diffusion timescale 

2-4 orders of magnitude slower than 
timescale of vibrational dynamics 


SQL7J1QN 

^ "" ■ - ■ ■ ■ — . — — 

• Off-line thermal transient ¥£. analysis 

provides coning as a function of time 

• Coning -time cuve sampled* at user-specified 
intervals to provide quasi -steady coning 
for vibrational analysis 
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dil^ilKiQ© Si&lL 
W (PI 


Bellows Damper 



Stator Assembly Mating Ring 


• Plain-faced balanced mechanical 
seals 

• Bellows-type flexibly-mounted 
stator ring 

• Rigidly-clamped rotor ring 

• Coulomb friction vibration dampers 


hH /Vm " v 


I 




I 


1 

I 


I 

I 
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TEST PROGRAMME 

SCOPE 

• 42 tests covering. 

- face coning 

- rotor eccentricity 

- rotor out-of -squareness 

- degree of damping 

• Effects of high and low temperature investigated 

• Tests at high speed (60.000 rpm) 

• All tests on typical plain face seal, modified to 
suit required conditions (face diameter = 30 mm) 


EAC1 UTT 

— 

• Hot gas (to around 220°C) supply 

• Cold gas supply (Polling to room temp) 

• Liquid cryogen supply possible 

• High pressure (rated to 20 bar) up to 14000 rpm 

• High speed (60.000 rpm) at lower pressure 


room temp 
14000 rpm 
5.10.15 bar 
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I£SI-PRQ£RAMM£ 

MISALIGNMENT AND ECCENTRICITY 



IBH 

fof/rGroup 




Accural* tracking of runout 



Stability raached quickly 
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MIAS SMUTTAIBIL 




m 


It® IPffifli 

M7 


• Liquid lubricant film 

• Cavitating film 

• Mechanical contact 

• Circumferential EHD 

• Different spring and secondary 
seal types 

• Floating rotor types 


• Ring seal geometries 
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Dynamic Coefficients for Multiple Brush Seals 

From notes and discussion of talk by D. Childs Texas A&M Univ. 
by R.C. Hendricks 

The cross sectional view of the Dressman-Chi Ids apparatus is given as figure 1. 
This is the same apparatus as described by Childs et al., ref. 1, and used to 
determine the rotordynamic coefficients for a variety of seal configurations and 
test conditions. The working fluid is compressed air. 

The experimental test sections for four, five, and six brush seal configurations 
are illustrated in figures 2, 3, and 4 respectively, see ref. 2. 

While the dynamics of the 4-stage system was reported by Connors et al . ref. 3, 
with good reliability, the results of the extensions to other brush seal sets has 
not proven satisfactory. The data are inconsistent, e.g., the leakage for 
4-brushes is less than for 6-brushes, however for all cases, the whirl frequency 
ratio was near or less than zero. But again the performance of the 4-brush con- 
figuration was dynamically better than that of the other configurations. 

These results are not yet resolved and as such are given to provide the Seals 
Workshop with some information about multiple brush configurations - at least be 
aware that multiple brush seals require attention even though dual brush seals 
are operating with high resolve, there are only a few results for more than two 
brushes in the literature and NO other dynamic data to determine or assess 
rotordynamic coefficients. 


For more information 


1. Childs, D.W., Nelson, C.C., Nicks, C., Scharrer, J., Elrod, D., and Hale, K., 
1986, "Theory Versus Experiment for th Rotrodynamic Coefficients of Annular Gas 
Seals: Part 1-Test Facility and Apparatus, "Trans. ASME J. of Triboloqy, V. 108, 
pp. 426-432. 

2. Griffin, M., Kelynhans, G., Alexander, C., Pierce, T. and Childs, D.W. 
Experimental Rotordynamic Coefficient Results for a 4-Stage Brush Seal 

TL-SEAL-17-92 #363, May 92. 

Experimental Rotordynamic Coefficient Results for a 5-Stage Brush Seal 
TRC-SEAL-7-92 #353, May 92. 

Experimental Rotordynamic Coefficient Results for a 6-Stage Brush Seal 
TRC-SEAL-8-92 #354, May 92. 

Turbomachine Laboratories, Texas A&M Univ., College Station Texas 77843 

3. Conner, K.J., and Childs, D.W.: AIAA 90-2139, "Brush Seal Rotordynamic 

Damping Characteristics" AIAA/SAE/ASME/ASEE 26th Joint Propulsion Conference 
July 16-18 1990 / Orlando, FL. 
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Table 1 Test Points 


Rotor 

Speed 

CO 

(rpm) 

Inlet 

Pressure 

P, 

(bar) 

Pressure 

Ratio 

K 

(-) 

Inlet 
Preswirl 
in the 
Direction 
of 

Rotor Rotation 

1 - 5000 

1 -7.90 

1 - 0.55 

1 - None 

2 - 12000 

2- Brl 

O a An 
X v.'n; 

2 * Intermediate 

3 - 16000 

3 - 18.3 

3-0.25 

3 - High 



4-0.14 




4— Stage Brush Seal 

Stator 


Brushes 

Iine-to-line contact 



Rotor 

0 129.4 mm 


quills 


3.175 mm 


Figure 8 Geometry and dimensions of seal insert, all dimensions in millimeters 
Fig. 8 illustrates the brush seal insert geometry used in this study. 
Uncertainty Analysis 


An experimental uncertainty analysis based on Holman (1978) for the measured 
experimental values is used in this study. The expression for the total uncertainty of a measured 
value is given below in equation (8). 
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5-Stage Brush Seal 



Figure 8 Geometry and dimensions of seal insen. all dimensions in millimeters 
Fig. 8 illustrates the brush seal insen geometry used in this study. 

Uncertainty Analysis 

An experimental uncertainty analysis based on Holman (1978) for the measured 
experimental values is used in this study. The expression for the total uncenainty of a measured 
value is given below in equation (8). 


6-Stage Brush Seal 



Figure 8 Geometry and dimensions of seal insert, all dimensions in millimeters 
Fig. 8 illustrates the brush seal insen geometry used in this study. 

Uncertainty Analysis 

An experimental uncenainty analysis based on Holman (1978) for the measured 
experimental values is used in this study. The expression for the total uncertainty of a measured 
value is given below in equation (8). 
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RESULTS OF CRYOGENIC BRUSH SEAL TESTING 


Presented by J. Scharrer, Rocketdyne 


and 


R.C. Hendricks, National Aeronautics and Space Administration 


Five-Brush Seal Configuration Leakage as a Function of 
Pressure for Selected Rotor Speeds 



Schematic and Section View of Multifunction Tester 


Bearing 

inlet 


Labyrinth seal 
resistance 
(replaced 
brush seal) 


(dp) -=:-|mPOT ii 


in 


-H/S bearing 
orifice resistance 



H/S bearing 
land resistance 


ATM 


MDOT 1/2 1/2 

labyrinth MDOT in MDOT in 

( dp) -: :-|MD0T ( dp) 


-Labyrinth 

seal 

resistance 


j- Hydrostatic 
bearing 


ATM ATM 

MDOT labyrinth = 1/2MDOT in - MDOT drain 
(Note: Some probes are not shown for clarity.) 



Labyrinth 

seals 
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Modeling bristle lift-off in Idealized Brush Seal Configurations 


Vijay Modi 

department of mechanical engineering 
Columbia university 
NEW YORK, NEW YORK 10027 


1. Introduction 

In the last decade, brush seals have emerged to be one of the most promising 
technologies for the reduction of leakage flow in gas turbine engines. Recent bench 
tests indicate a possibility of an order of magnitude reduction in leakage flow over 
multiknife labyrinth seal, Holle and Krishnan (1990). Additional potential for 
performance benefit arises from mechanical and maintenance considerations. The 
efficiency of a labyrinth type seal depends upon the clearance between the tip of the 
knife and the bore but this radial clearance may be difficult to control due to thermal 
and dynamic conditions. A brush seal on the other hand is compliant, and hence has 
the ability to recover after excursions. Flower (1990). An important question that 

remains to be answered is the relationship between brush configuration/operation 
parameters and some measure of its compliance. One measure of compliance is the 

clearance that develops between the bristle tips and the rotating element due to the 

pressure differential across the seal and due to the aerodynamic drag. 

We attempt in this paper to develop a model for the flow through brush seals and 

determine their elastic behavior in order to predict the dependence of brush/joumal 

clearance on geometry and operating conditions. Several idealizations regarding 

brush seal configuration, flow conditions and elastic behavior are made in the 
analysis in order to determine closed form parametric dependence. This formulation 
assumes that there is no initial interference between the bristle tip and the rotor. 
Also interbristle, bristle-backing plate and bristle-rotor friction is neglected. The 
bristle bundle or the brush seal as it is alternately called is assumed homogeneous 
and isotropic on a macroscopic scale so that a physical property like permeability is 
uniform. The fluid is assumed to be homogeneous, incompressible, viscous and is 
flowing under steady conditions. 

A schematic of a brush seal is shown in figure 1. If the nominal bristle-shaft 
interference is absent then under static conditions the bristles may deflect axially 
due to the imposed pressure differential. This axial deflection may create a clearance 
permitting leakage flow in excess of that which occurs through the porous matrix 
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formed by the bristle bundles. Under dynamic conditions the Couette flow created by 
shaft motion could be strong enough to cause bristle deflection and once again a 
clearance may develop. 

The paper proposes a means to determine this clearance (or at least describe its 
parametric dependence on geometry and operating conditions) under static as well as 
dynamic conditions. The study can be thought of as consisting of three separate 
modeling efforts. First a flow model that describes the coupled parallel flow through 
the porous medium made of bristle bundles and the clearance region. This 
development follows the earlier work of Beavers and Joseph (1967), Beavers et al. 
(1970), Williams (1978) and Rudraiah (1985). This model provides the macroscopic 

description of the flow field, i.e., a filter velocity in the porous medium. Second, a 

model to relate this macroscopic flow field to a local flow field and its associated drag 
on the bristle is developed. This model also permits us to determine an expression for 
the permeability in order to characterize the physical property of the bristle bundle 
in absence of an experimentally determined value. The forces on a single bristle due 
to this macroscopic flow field are then estimated assuming idealized microscopic flow 
behavior and a phenomenological description of drag on a bristle. Third, the elastic 

behavior is modeled to estimate the deflection of the bristle tip due to the flow field or 
the impressed axial pressure differential. A description of the clearance would in 

principle require a simultaneous solution of the three models. It is, however, 

assumed here that the physical property permeability remains unchanged both due 
to the presence of leakage flow through and around the bristle bundle as well as due 
to the deformation of the bristles themselves. Given this, it is then necessary to solve 
only the flow and elastic models simultaneously to determine the clearance. 

2. Circumferential Deflection 

To determine the circumferential deflection of the bristles we must first estimate 

the forces on each bristle. The forces in the circumferential direction are due to 

drag caused by the flow around the bristle bundle. This flow is in turn driven by the 
Couette flow in the bristle-shaft clearance as a result of the circumferential slip 
velocity of the shaft itself. In the absence of a clearance, the shaft motion imposes a 
circumferential velocity boundary condition directly at the bristle-shaft interface. 
A model to determine the flow field in both these circumstances is introduced next. 

This development assumes that the clearance h 0 if any is known. Its actual value will 
be determined later along with considerations of the elastic behavior. 
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2.1 Flow Model 


Let us consider the geometry shown in figure 2, with a porous medium of height 
h underlying a channel formed by a clearance of height h G . The clearance is 
bounded above by an impermeable wall moving to the right at Uo and the porous 
region is bounded below by an impermeable stationary wall. The formulation 
presented here follows that in Rudraiah (1985). The basic equations describing the 
flow are obtained after the following approximations are made. 


i) The fluid is homogeneous and incompressible. 

ii) The flow in the channel and in the porous medium is driven by a shear produced 
due to the motion of the upper plate. This flow is steady, laminar and fully 
developed. 

iii) The porous medium formed of bristle bundles is homogeneous and isotropic on a 
macroscopic scale. 

iv) The flow in the porous medium is adequately described by the Brinkman 
equations and this flow is coupled to the channel flow by a boundary condition 
given by Williams (1978). 

A. 

Following Rudraiah (1985), we write equations for u, velocity in the gap and u, 
filter velocity in the porous medium as: 



d u 1 

and ^-XE“ = ° 


( 1 ) 


where X is a positive constant and k is the permeability of the porous medium. The 

A 

velocity u in the porous medium is related to u by 


u = (1— <|>) u (2) 

where ( 1 -<J> ) is the porosity. Following Williams (1978), at the clearance-porous 
medium interface we assume that: 


du „ du 

ar <*-♦*$ 


(3) 


The remaining boundary conditions are the usual no-slip conditions at 
impermeable walls and are: 
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u = u 0 at y = ho 

u = 0 at y = -h (4) 

Solving (1) subject to (2), (3) and (4) we obtain the velocity distributions in the 
clearance and the porous regions to be 


M (h Q -y) 1 

(tanhSh + J 


y >0 


* f . , o cosh Sy ~| 

U ~ (1 — c}))(tanh6h + X8h 0 ) l n y + coth 8h J 


y < 0 


(5) 


where 5 = (Xk)' 1 / 2 . Here 8h = h (^.k)' 1 / 2 » 1 is used. This relies on X to be of order 

unity and h to be a macroscopic length scale assumed to be several times greater than 
k l /2 (which is typically of order d 0 , the characteristic dimension of the porous 

A 

matrix, i.e. bristle diameter). Hence we may approximate u by 


u = — e 5y y<0 

^ ( 6 ) 

where T\ = (1-<(>)(1+XS h 0 ). The exponential behavior of the filter velocity implies that 

it decays to the Darcy value (which in this case is zero) within a boundary layer of 
length scale 1/8, where 1/8 is of the order of k^ 2 . 

A 

Here u represents a filter velocity, a macroscopic quantity defined in order to 
avoid the more difficult question of what is the true velocity of the fluid in the 
porous region between the bristles. We now make certain idealizations about the 
bristle bundle geometry and subsequently model the flow through the interbristle 
pores in order to determine the viscous drag force directly as a function of the filter 
velocity. 

2.2. Permeability Model 

The drag force on a bundle of cylindrical bristles clearly depends upon the flow 
through the pores which in turn depends upon the particular geometric 
configuration. Let us first examine the situation for two particular geometric 
arrangements. Let us assume that the cross-section of the bundle remains same in 
the direction along the bristle axis and that the flow is normal to these axes. Let <f> be 

the solidity or here the area fraction and Z be the number of nearest neighbors. 
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Then <|> and Z are respectively Jt/4 = 0.79 and 4, for a square array and tc/2 V3~ — 0.907 


and 6 for a hexagonal array. For these arrays there is no possibility of any 
transverse flow since all neighboring cylinders are in contact. Typical solidities <}> 

for brush seals are between 0.7 and 0.8, indicating a fairly close packed geometry. It 
is very likely that the manufacture of bristle bundles and their relative movement in 
the presence of leakage flow lead to bristle configurations that are close to random. 
A particular means to generate a closely packed random array is by the following two 
step process, described by Sangani and Yao (1988). In the first step, the process of 
dropping a large number of equal-diameter cylinders in a container is simulated. 
Note that in the configuration so generated any cylinder is in contact with its 
nearest neighbors. Let <}) t be the solidity of such a closely packed random array. 
Berryman (1983) summarizes the results of several simulations and reports <> t in the 
range of 0.81-0.89, with most studies quoting a value of approximately 0.82. Sangani 
and Yao (1988) have also simulated such arrays recently with upto 1600 cylinders in a 
container and they report a value of <j>t = 0.824 and Z = 4.2. In the second step, the 
diameters of all the cylinder centers are shrunk by a constant amount 2er, while 

keeping the location of all cylinder centers fixed. Here r is the bristle radius. The 
value of E is chosen so that the new configuration has the measured solidity <{>. For 

this special kind of array the gap between neighboring cylinders is uniform and is 
given by 2£r where 


The resistance offered to the flow by the gap between pairs of nearly touching 
cylinders determine the overall drag and hence the effective permeability of the 
medium. The analysis below is due to Sangani (1990). The pressure drops as fluid 
squeezes through a gap of width 26 r at a volume flow rate Q, as shown in figure 3. 

Here Q is the two-dimensional volume flow rate. The profile of the cylinder surface 
can be approximated by 



(7) 



( 8 ) 
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Assume that inertia terms are negligible and that the viscous forces in the 
streamise direction are small compared to those in the normal direction. Hence the 

A 

equation governing the x- direction velocity component u reduces to 


5 2 u dp 


(9) 


A A 

Integrating the above equation with u = 0 at y = f(x) and du/dy = o at y = o we 
obtain 

zp. cue (io) 


Integrating the above velocity profile over the gap and identifying it with 
volume flow rate Q, we obtain 


(Ap) gap = 


9ttp.Q 

8V2 


£ -5/2 


( 11 ) 


Thus the force exerted per unit length and width in direction of flow is Ap. The 
direction of flow however is normal to the line segment joining the centers of the 
cylinder pair in question. If such a line segment is oriented at an angle 0j as shown 

in figure 4, then the force vector Fj per unit length arising over width 2r due to each 
gap for volume flow rate Qj is given by 



97tpQ: E;-5/2 

1W ~ <2r) 


( 12 ) 


The component of the force (per unit length) in the direction of mean flow 
(assumed to be along 0 i=0) is then obtained after noting that Qi = (2 usin0j) r 


971LLE ^ 

F h = IFjl sin0j = — ~= - 2usin 2 0j 

8V2 (13) 

Here u is the filter velocity or the superficial velocity. Hence, the mean force <F]> is 
given by 


<Fl> = 2^e- 5/2 

8V2^ 


(14) 
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where Z is the average number of nearest neighbors and the mean value of sin 2 0, 
averaged over all is 1/2. If the number of cylinders per unit cross-section area is 
4<t>/7td 0 2 where do is the diameter then 


IVp,l = <F 1 >4=^e- 5 ' 2 
1 *<i 0 2V2 


do 2 


(15) 


If one was interested in the permeability k of the medium defined as k = pu/IVp I then 


k=— <£e 5/2 
9Z<j> °° E 


(16) 


which is identical to the expression in Sangani and Yao (1988) except for <]) in the 
denominator. 


2.3 Elastic Behavior Model 

The expression in (15) permits us to calculate iVpl, as a function of the filter 
velocity u which may vary along the bristle axis as given by (6). Hence, we are in a 
position to determine the loading on the bristle due to the inter-bristle flow driven 
by the Couette flow in the clearance region. We define a co-ordinate system shown 
in figure 5. The deflection of the bristle tip, At, can be obtained from 

straightforward application of linear elasticity theory. The radial component of the 
deflection ho = A(CosP is then the clearance, and is given by 


h -Acosp-— c~ 5/24U ° COSp 1 1 h3 

° * ^ 2V2 £ t, (sin p) 3 El 8 3 


(17) 


Recall that T| = (l-<t>) (1+X5h 0 ), 8 = i = 7ldo^/64. Since the above equation is 

implicit in ho, consider ho » k^ 2 first, so that A,5ho»l providing an explicit form. 


2 



-5/2 4U 0 k cos (3 

(l - ^) (sin|3) 


h 3 
3 3EI 


(18) 


If permeability k is experimentally determined for the medium then that 
measured value can be used in the above expression. The recommended values of <|> t 
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and Z arc 0.82 and 4.2 for a random array of cylinders. Hence the expression is valid 
for <J) < 82). If, however, k is not known, permeability can be estimated from (16) 

allowing us to determine h 0 from 

^ cos|J h L _ /t* 

oJ '«»-*> ® ^ (19) 

In the absence of any clearance h Q vanishes and any initial increase in h 0 must 
occur with h G « k 1 / 2 . Under circumstances that permit us to assume X8h 0 «1 we 


obtain 


_ 9Z p 1 ^7 , — ~]~5/2 |4u 0 VXk cosP h 3 

£=^L'-VW,J _____ 


( 20 ) 


Once again k can either be measured or evaluated from equation (16). In the above 
expression in addition to Z and 4>t we also need an empirical estimate for the constant 

X 1/2 . The constant \ 1/2 can be identified with a in Beavers and Joseph (1967) and is a 

dimensionless quantity depending on the material parameters which characterize 

the structure of the porous medium within the boundary region where the filter 
velocity decays to zero, its Darcy value. The value of a reported by Beavers and 

Joseph are for Foametal and Aloxite, with effective pore sizes varying between 0.013 
inches and .045 inches. The value of OC for these materials was found to vary between 

0.1 and 4 with lower values observed at lower pore sizes. 

Preliminary estimates of the clearance due to tangential loading can be made 
from either (19) or (20) and if it turns out that h 0 is of the order of k 1 ^ 2 then the 
quadratic equation (17) can be solved for h 0 . Before discussing the results for the 

tangential deflection we will develop the analysis for axial deflection. This will 
permit us to determine their relative magnitudes and establish conditions under 
which, it may be possible to neglect the clearance due to deflection in one of the 
directions. 


3. Axial Deflection 

The axial deflection of the bristles is due to the pressure differential along the 
axis of the rotating element. It is the purpose of the brush seal to minimize what 
would otherwise be a leakage flow due to this pressure differential. The axial loading 
on the bristle is straightforward to estimate since the pressure differential impressed 
upon the bristle bundle can be assumed to remain unaltered in the presence of 
leakage flow. If the pressure differential is Ap = Ph-Pe ovcr the width w of the 
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bristle bundle (in the axial direction) then the force per unit length, q, acting on the 


length a of the bristle, i.e. exposed portion between the backing ring of inner 

diameter Db and the shaft, diameter D s . Since the bristle bundle is clamped at the 
retaining plate inner diameter D r , the bristle behavior can be modeled as a bar 
clamped at origin and simply supported at Db with an overhanging distributed load 

between Db and D s . If we define L to be the bristle length between the retaining 

ring ID and the backing ring ID, L = (D r -Db)/2 sin J3 then the bristle geometry with 

its axial loading diagram is as shown in figure 6. We wish to determine the 

displacement. Ah, due to the axial deflection of the unloaded member, as shown in 

figure 7. This quantity is the clearance produced due to the axial loading and would 

be observed in a static leakage test. In dynamic tests the clearance would be 

produced by a combination of effects, the tangential as well as the axial loading. 

Linear elasticity theory assumes small angles of rotation for the beam and thus 
would only permit the calculation of the vertical displacement A v while Ah would 

remain zero. Typical axial loading and the bristle length to diameter ratios are 

however, such that it becomes necessary to use large deflection theory in order to 
model the problem. Thus the differential equation of the deflection curve becomes 


The exact shape of the elastic deflection curve given by the solution to this 
equation, is called the elastica. The mathematical solution to the problem of 
determining the elastica has been obtained for many different types of beams and 
loading conditions, see Frisch-Fay (1962). The solutions to the specific loading of 
interest here was not available in the existing literature. A particular difficulty is 
that the problem is statically indeterminate and the reaction force where the bristle 
is simply supported is unknown. While a linear theory may provide a value of the 
reaction, it will be approximate at best. This value, however, could be of use as an 
initial guess in an iterative determination of the reaction force into vertical 
component R and horizontal component Rtan 0, where tan 9 is the bristle slope at 
the support, see figure 7. To determine Ah we need to solve for the elastica and 


bristle is given by q = Ap d Q 2/w. This force acts uniformly over the overhanging 
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terminate the curve at a point where the length of the curve is (L + a). Here we 

assume that any elongation of the bristle is negligible or of an order smaller than 

that under consideration. In the deformed position of the bristle the axial loading is 
no longer applied over the overhang initial length 'a' but over the overhang length 
a' = a - Ah in the deformed position. The direction of the axial force, however, 

continues to remain vertical in the deformed position. 

Solution to the problem described above requires us to integrate (21), with a 
moment distribution function M(x) given by 

M(x) = qa' (a’/2 + L-x) -R(L-x) -RtanBy 0 < x < L 

M(x) = (q/2) (L + a' - x) 2 L <x < L + a (22) 

where a' and tan9 are unknown and are determined as part of the solution y(x). The 
reaction force, however, is still unknown because the problem is statically 

indeterminate. This difficulty is easily overcome, especially since the solution is 

perhaps most easily obtained on a computer. An initial guess of R is made from 

linear theory. The preliminary solution thus obtained however, will in general fail 

to pass through the support at (L,0). The magnitude of the reaction force is gradually 

changed until the elastica does pass through, (L,0). It is convenient to normalize x 

and y with length L, i.e. x* =x/L and y* — y/L. Dropping asterisk now equations (21) 
and (22) become 

y7(l+y' 2 ) 3/2 = -M(x) with y(o) = o y’(o) = 0, where 
M(x) = C e ff { 1/2 + (1/a’) (l-x)} - a (l-x)- atan9y 0<x<l 

M(x) = (C e ff /2) (1/a') 2 (1 + a’-x) 2 i< x < 1 + a/L (23) 

where C e ff = C (a’/a) 2 , C = qa 2 L/EI and a = RL 2 /EI 

The problem shown above was solved numerically for various values of the 
loading parameter C and the overhang ratio a/L. The elastica curves for a/L = 0.2 are 
shown in figure 8. Note that the axis in the y-direction is stretched considerably for 

clarity. The quantity of interest here is the clearance Ajj the variation of which with 

loading parameter C is shown in figure 9 for several overhang ratios a/L. Note that 

the small deflection theory value of Ajj is identically zero for all C. We observe that 

the departure from this value as loading is increased depends upon the overhang 
parameter. It turns out that for typical gas turbine applications brush seals may be 
operating in a parameter range where Ah is rapidly increasing with C and a/L. 
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The actual value of the reaction force is not of immediate interest once the axial 
deflection is known, however, it may be of use in the following manner. Earlier 
analysis to determine the circumferential deflection assumed that the friction at the 
bristle-backing ring interface was negligible. The bristle orientation under axial 
loading shown in figure 8 indicates that the only point of contact may be at the inner 
edge of the backing ring where the normal force acting on the bristle is given by 
R/cos 0. An estimated value of the static Coulomb friction coefficient would permit 

an approximate evaluation of the tangential restraining force. This force would act 
on the bristle layer in contact with the backing ring. While in this paper we do not 

account for this friction force in the calculation of the circumferential deflection 

the value of R and 9 may prove useful for future work. 

4. Discussion 

Under dynamic conditions the actual deflection would be due to both the 

circumferential and axial loads and a vector sum of forces has to be used to solve the 
three-dimensional elastica problem. Additional simplification is possible if we can 
demonstrate that the deflection in one of the directions is small compared to the 

other. An explicit expression for their relative magnitude is not available because of 

the lack of a closed form formula for the axial deflection. 

The numerical results of figure 9 establish that the behavior of Ah (C, a/L) is 

indeed nonlinear, and may be of importance in the design of brush seal systems. In 
addition to the expected non linearity in the response to loading parameters C the 

displacement h is also sensitive to a/L ratios. Thus it is not just the overhang length 
that is relevant even though it is the length exposed to a pressure differential but 

also the size of the retaining plate. The analysis carried out here while making 
several idealizations may provide insight into the dependence of the clearance on 

geometry and operating conditions. 
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Packing Thickness 



Diameter D, D 

r b 8 

Figure 1 Schematic of Brush Seal System. 



Figure 2 

Physical model of clearance channel 
bounded by shaft on one side and 
porous bristle bundle on other. 
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Figure 3 Schematic of Volume Flow Q 
through a gap of size 2tir 



Figure 4 A pair of neighbouring cylinders 
in a random array. 



Figure 5 Schematic diagram of clearance region and the porous 
bristle bundle matrix. Velocity distributions in each 
are shown at left. Bristle displacement is normal to its 
axis. 


h 0 = 4 t cosp 
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Figure 6 Axial loading diagram for a single bristle. 



the length of the bristle is unchanged under loading. 
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Figure 9 Clearance generated due to axial loading for 
various values of C. a/L» « 0.2,0. 4,0. 6 




Figure 10 Variation of the (formalized reaction force 8* 
where »* -(RL. 2 /£I)/ C eff 
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fiEAE Brush Seal Program Presentation Summ ary 


Bru s h Sea l s for Turbine En aine Fu el Conservation 

GeneralEteehte^r^ft Engines (GEAE) has an ongoing program sponsored by NAWC 
(Contract No. N001 40-90-C-31 99) for the evaluation of brush seals. This is a summary 
of the program status as presented during the Brush Seal System Workshop at the 
NASA Lewis Research Center in Cleveland, Ohio on August 6, 1992 by Ron Korzun. 

The objective of the program is to design, prooire and test brush seals intended to be 
suitable for replacing the inner and outer balance piston seals in the T407 turboprop 
engine. An existing T407 low pressure turbine rig was modified to conduct the brush 
seal testing. The rig which includes the power turbine shaft accurately simulates the 
dynamics of the engine. A picture and cross section of the rig are shown in Figures 1 
and 2. ® 

Eia ure 1 . T4Q7 R i g Assemb led In last Platform 



GE Aircraft Engines 

Figure 2. T407 IBP and QBP Brush Seal D ynamic Test Rip 

OBP Seal Discharge Cavtty 



Boure 3. Seal Desion and Fit -I Ip 


Pulw Grooves To Faciltate 
Removal Ot Seals 



IVxmh Seal Inner Diameter Cold (B) 

tVueli To Runner Diametral Inlerference - Cold (0) 

Brush T 0 Runner Diametral Interference • Hoi (8) 

Brush Axial Pack Width [C] 

Booking Plata GAP (0) 

CHameirai Interference FH With Stator Suppon - Cold (A) 
•Maximum Stress - Seal Support 
'Maximum Stress -Brush Seal 
' At SS IRP, Nominal Fit Up 


one 

IBP. 

to r 

5.6- 

006 

006 

.010 

.010 

.025 

.050 

0$1 

.044 

006 

.004 

18KSI 

25 KSl 

13KSI 

17KSI 
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GE Aircraft Engines 

Fiflura 4. Kev_Brush Seal/Ria Features 



The test plan calls for cyclic test to evaluate the leakages at operating pressures and 
temperatures. The test plan also calls for the evaluation of the effect of seal pack width 
and rotor eccentricity. The brush seal design goals are to achieve or exceed a 50% 
leakage reduction versus labyrinth seals. This must be accomplished without affecting 
rotor life. The test will also evaluate runner coating wear, heat generation and 
corrosion. 

The project is carried out under the direction of Bill Voorhees, the NAWC Project 
Engineer with the assistance of Guy Ullman, the NAVAIR Project Engineer and Carl 
Graia, the NAWC Component Engineer. At GEAE Warren Ostergren is the program 
manager. Ron Korzun is the design engineer and Bert Campbell is the test engineer. 
EG&G Sealol's support has come from program manager, Chuck Nevola and deskm 
engineer, Bob Johnson. 

As of the August review, the build of the test rig for its initial testing was in process 
GEAE looks forward to presenting the results of the test program during a future NASA 
workshop. 


Program Objective 


Design, Procure, And Test Brush Seals Intended To 
Be Suitable For Replacing The Inner And Outer Low 
Pressure Balance Piston Labyrinth Seals In the T407 
Turboprop Engine. 
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Program Approach 


• Design And Procure Brush Seals 


• Modify And Instrument An Existing T407 Low 
Pressure Turbine Test Rig 

• Conduct Cyclic Tests To Evaluate Seal Leakage 
At Operating Pressures And Temperatures 


• Evaluate Effect Of Seal Pack Width And Rotor 
Eccentricity 


Project Organization 
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Brush Seal Design Goals 


• Achieve Or Exceed Predicted Leakage 
Reduction Versus Labyrinth Seals 

- OBP - 59% Reduction 

- IBP - 45% Reduction 

• Rotor Part Life Not Limited By Brush Contact 

- Runner Coating Wear 

- Heat Generation 

- Corrosion 


T407 IBP And OBP Brush Seal Dynamic Test Rig 

OBP Seal Discharge Cavity RIG OPERATING BOUNDARY CONDITIONS 
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Power Turbine Rotor Assembly 



Exhaust Frame Stator Assembly With 
Brush Seals Installed 
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Exhaust Frame/Outer Housing Assembly 



T407 Rig Assembled In Test Platform 
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Key Brush Seal/Rig Features 



Spacer Arms Used To Position 
Seals In 2 Locations On Runner 


Chromium Carbide Runner Coating 
Compatible With Haynes 25 Bristles 


Axial Retention Sleeves With 
Locking Pin 


INCO 625 Support With HAST-X 
Brush Seal Backing Plate^ Gives 
Best Growth Compatibility To 
Avoid Errant Rig Leakages 


Outer Balance Piston Brush Seal 
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Inner Balance Piston Brush Seal 



Seal Design and Fit • Ups 


Puller Grooves To Facilitate 
Removal Of Seals 





IBP 

Brush Seal Inner Diameter Cold (B) 

io.r 

5.6" 

Brush To Runner Diametral Interlerence - Cold (B) 

.006 

.008 

Brush To Runner Diametral Interference - Hot (B) 

.010 

.010 

Brush Axial Pack Width (C) 

.025 

.050 

Backing Plate GAP (D) 

.051 

.044 

Diametral Interference Fit With Stator Support - Cold (A) 

.006 

.004 

‘Maximum Stress - Seal Support 

18KSI 

25 KSI 

‘Maximum Stress -Brush Seal 

13KSI 

17KSI 

‘ At SS IRP, Nominal Fit Up 




• Brushes Maintain Contact With Rotor At Alt Operating Conditions 

• Backing Plate Distance Sized For Worst Case Conditions Expected In Field 

• Backing Interference Maintained At All Operating Conditions - Avoid Leakage 
■ Stress Is Acceptable - Below .2% Yield Strength 
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Brush Seal Data Used To Size 
T407 Brush Seals 


Test Conditions 
Pressure Drop 
Temperature (F) 
Surface Speed (Ft/S) 

Seal Parameters 
Diameters 

Backing Plate Height 
Pack Thickness 
Bristle Diameter 
Bristle Angle 
No. Of Stages 


GEE 

Sealol 

0-80 

15-25 

70 

70 

430 

375 


5.1, 11.3 

5.4 

.04-. 130 

.047 

.018-.027 

.025 

.002..0028 

.0028 

45-56 

36-56 

1 

1 


Technical Literature 

20 

600 

668 


5.1 

.02..03 

.027 

.0028 

45 

1.2 


• The above ranges stated are for the test data which was readily available (10/91 ) 
to use in evaluating the T407 Seals 

• Additional studies are continually in process at GEAE (Ruston) and at Sealol 


Brush Seal Test Data - Flow Parameter 
Versus Pressure Parameter 
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Test Plan Schedule 


1992 


I. Break In T esting - 1 st Set 


50 hrs 




t 

150 hrs 


III. Cyclic & Endurance Testing - 2nd Set 

IV. Alternate Seal Design Testing 

V. Eccentricity Testing - 3rd Set 


150 hrs 





50 hrs 


50 hrs 


Test Plan Includes Extensive Testing On Three Sets 
Identical Seals And One Alternate Design. 


Summary 

• OBP And IBP Brush Seal Designs Completed 
And Fabricated 

• LPT Dynamics Rig Modifications Completed 

• Test Plan Defined 

• Test Cell Checkout Completed 

• Testing To Begin In August 
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Tribopair Evaluation of Materials for Brush Seal Applications 

From notes and discussion of talk by J. Derby , EG&G 
by R.C. Hendricks 

Derby and England investigated several brush seal system tribopairs from bristle 
materials Haynes 25, Inconel 718, and proprietary materials called Alloys A, B, 
C, D, and E and coatings chromium carbide, Triboglide, tungsten carbide, chromium 
oxide, aluminum oxide, Tribaloy, Tribomet T-104C. 

Haynes 25 is a cobalt chromium nickel tungsten based superalloy that forms an 
oxide resistive 0203 film. 

Inconel 718 is a precipitation hardenable nickel chromium superalloy good to 1300 
F. 

Alloy A (Ni Cr A1 base superalloy solid-solution strengthened) is being used in 
gas turbine hot spots and develops a tenacious 0203 and A1203, yittria-modified 
oxide layer. 

Alloy B (solid-solution, and carbide strengthened Ni Cr, W based superalloy) with 
good strength at elevated temperatures developing an 0203 oxide film that is 
enhanced by lanthanum, and used in combustion cans. 

Alloy C (mechanically alloyed, oxide dispersion strengthened, Fe, Cr, A1 based 
superalloy) has excellant resistance to oxidation due the a tenacious Y203 sta- 
bilized A1203 and 0203 film, good to perhaps 2100 F. 

Alloy D (solid-solution strengthened by addition of W and precipitation of M6 
and M23C6 carbides Co, Cr, W, Ni based superalloy) is one of the first sheet alloys 
developed for aircraft engines. Adding lanthanum is reported to modify the pro- 
tective oxide layer up to 2000 F. 

Alloy E (age hardened Ni, Mo, Cr based superalloy without gamma prime precipitate 
and hardening believed based on Ni2(cr,Mo) stoichiometry) with service temper- 
atures to 1400 F. 

Although not reported Haynes 125 is known to have a higher use temperatue but more 
data are required on tribological pairing. 

Chromium carbide has superior wear and friction performance. Sprayed coatings 
contain Cr3C2 and Cr7C3 mixtures produce an increase in hardness. The CrO was 
applied using HVOF. 

Triboglide is a CrC containing additives of 12 wt. percent barium and calcium 
fluoride solid lubricants based on the work of Sliney at NASA-LeRC. Spray pa- 
rameters (standoff, angle, powder size, feed rate, temperature) were investigated 
and use of HVOF was prefered. 

Tungsten carbide contains a Co binder and exibits phases WC, W2C, and Co3W3C with 
decarbonization of some WC above 900F. Union Carbide D~Gun was preferred with 
8-percent Co binder. 
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Chromium oxide has higher friction but used within the industry has low thermal 
expansion and thermal conductivity. 

Aluminum oxide was applied using air plasma spray (APS) with a bond coat of 

0. 001.0.002 inch. 

Tribaloy T-800 is a Co based with Mo, Cr, and Si additives and APS and HV0F applied 
followed by a heat treatment. 

Tribaloy T-104C is Co based that is electrodeposited such that CrC particles are 
codeposited in the coating. Exhibits good friction and wear characteristics and 
is almost entirely metallic ( good rotor adherence). 

The tests were performed using a 2.0 inch diameter, 0.040 inch-two-1 obed-cam rotor 
at 1200F at speeds of 10, 20, and 30 krpm, providing about 15E9 fatigue cycles. 
The interface was assumed to be line to line at the mean equivalent diameter al- 
though no preload information was given. 

In a second of tests at 0-pressure drop, conducted at 2, 20, 40, and 60 krpm at 
800 F the nominal brush-rotor clearance was -0.010 inch-radial. 

Inconel 718 and Haynes 25 bristles form unstable Cr203 layers with leading edge 
alloy fractures. Alloy A is least damaged by fatigue cycles and most oxidative 
resistant. 

Oxide coatings exhibit microfracturing and grain pullout which promotes wear and 
cobalt alloy bristles perfromed poorly due to a transformation from hexagonal 
close to face centered at 400C (750F). 

Alloy A/Triboglide (15-20 microinch Ra ) provided good friction behavior vs 
sliding speed (varies from 0.14 to 0.25) with Inco 718/Triboglide constant at 0.3. 
The Haynes 25 /Triboglide is good at surface speeds of 100 to 350 ft/s while Haynes 
25/CrC increased linearly from 0.23 to 0.45 at 500 ft/s. 

As an editorial note, there does exist some transition from hexagonal close pack 
to cubic crystal structure between 350 C and 400 C. Within this temperature range 
the wear increases and is up to ten times higher in a vacuum (ref. 2). In air 
the formation of oxides tends to mitigate wear and for low sliding speeds wear 
actually decreases as temperature is increased (ref. 3). However data at both 
elevated temperature and surface speed are needed. 

For more information 

1. Derby, J., and England, R.: AIAA 92-3715, "Tribopair Evaluations of Brush Seal 
Applications,” AIAA/SAE/ASME/ASEE 28th Joint Propulsion Conference, July 6-8, 
1992, Nashville, Tennessee. 

2. Buckley, D.H.; and Johnson, R.L.: Friciton and Wear of Hexagonal Metals and 
Alloys as Related to Crystal Structures and Lattice Parameters in Vacuum. ASLE 
Trans. 9, 121-135 (1966). 

3. Murray, F.S.; and Calabrese, S.J.: Low Speed Sliding Behavior of Metal-Ceramic 
Couples at Temperatures to 800 C. Lubrication Engr., V. 49, 5, pp. 387-397 
(1993). 
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NASA LEWIS RESEARCH CENTER SEAL WORKSHOP IN CLEVELAND/OHIO, 5 - 6 AUGUST 1992 


DAMPING SEAL BEARINGS 


George L. von Pragenau, Huntsville/Alabama, (205) 536-6832 
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PRESENTATION TOPICS 


1. ABOUT TECHNOLOGY 

2. HIGH PRES. FUEL TURBOPOMP (HPPTP) 

3. HIGH PRES. OXYGEN TURBOPUMP (HPOTP) 

4. HIGH PRESSURE TURBOPUMP WHIRL DATA 

5. BASIC ROTORDYNAMIC MODE 

6. ROTOR SPEED LIMIT 

7. DAMPING SEAL 

8. SPACE SHUTTLE PUMP SEAL HISTORY 


9. HPOTP DAMPING SEAL BEARINGS 

10. HPOTP DAMPING SEAL BEARING ENGINE TESTS 

11. HYDROSTATIC BEARING 

12. DAMPING BEARINGS 

13. HP0TP/HPFTP TURBINE DAMPING BEARINGS 
14. SEALS k BEARINGS DYNAMIC & LEAKAGE DATA 
15. GENERIC BEARING COMPARISON 

16. OUTLOOK 


ABOUT TECHNOLOGY 
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HIGH PRESSURE FUEL TURBOPUMP (HPFTP) 



HIGH PRESSURE OXYGEN TURBOPUMP (HPOTP) 
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HIGH PRESSURE TURBOPUMP WHIRL DATA 



BASIC ROTORDYNAMIC MODEL 



ROTOR WITH SPEED £1 
AND MASS M r 


BEARING POINT 


r » 





Ky- (oO ) 2 #f sC+s 2 M r Q+cO{C+s2M) 
■Q~cO(C+s2M) K z -(cO) 2 M+sC+s 2 M r 


0 ROTOR SPEED IN RADIANS 
S LAPLACE OPERATOR RADIANS 


F| EXTERNAL ROTOR FORCES y, z ROTOR DISPLACEMENTS 
C SEAL INDUCED DAMPING C=V/W COUETTE FLOW RATIO 
K i BEARING+SEAL STIFFNESSES M r ROTOR MASS PLUS M 
M SEAL FLOW INDUCED MASS Q OTHER WHIRL DRIVERS 

V COUETTE FLOW VELOCITY v ROTOR SURFACE SPEED 
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ROM SPEED LIMIT 



DAMPING SEAL 
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SPACE SHUTTLE PUMP SEAL HISTORY 


EAEY IN 1976 THE HIGH PRESSURE FUEL TURBOPUMP (HPFTP) SPEED LIMIT WAS 22,000 RPM TO 
AVOID A 175 HZ SUBSYNCHRONOUS ROTOR WHIRL. THE REQUIRED OPERATIONAL SPEED IS 36,000 RPM. 

STIFFER BALL BEARING SUPPORTS AND TWO STIFFER SMOOTH SHAFT SEALS (INSTEAD OF LABY SEALS) 
RAISED THE 1st CRITICAL SPEED TO 18,000 RPM AND THE WHIE TO AN INTERMITTENT 300 Hz. 

DAMPING SEALS WERE INVENTED 1981 TO ELIMINATE THE 50 \ SUBSYNCH. WHIRLS OF SMOOTH SEALS. 
DAMPING SEALS WERE INSTALLED 1989 IN THE HPFTP. COULOMB FRICTION DEVICES WERE ABANDONED. 

ONE HIGH PRESSURE OXYGEN TURBOPUMP (HPOTP) EXPLODED 8 SEPTEMBER 1977 FROM A 95 * SUB- 
SYNCHRONOUS WHIE, CAUSED BY INTERNE ROTOR FRICTION, THE MOST VIOLENT WHIRL-DRIVER. 

SINCE THE 1st REEIGHT, DISCOVERY ON 29 SEPTEMBER 1988, EVERY HPOTP EIGHT PUMP HAS 
DAMPING SEALS IN THE PREBURNER PUMP (PBP) INLET AND DISCHARGE TO AVOID THE 95 I WHIE. 

HIGH BALL NEAR LIMITS THE HPOTP TO 2 FLIGHTS AND IGNITED A FIE 23 JUNE 1989 IN A TEST. 
TESTS W/O BALL BEARINGS CONFIRMED THAT DAMPING SEALS AE SUPERIOR HIGH SPEED BEARINGS. 


HPOTP DAMPING SEAL BEARINGS WITH AND W/O BALL BEARINGS 
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HPOTP DAMPING SEAL BEARING ENGINE TESTS WITH BALL BEARINGS 



HPOTP PUMP END BEARING TEST SAMPLE 


HPOTP 

ENGINE 

TEST NO. 

TIME 

SUM 

0307# 

NUMBER 


SEC. 

SEC. 

R3 

2105 

901-501 

519 

519 

R3 


502 

520 

1039 

R3 


503 

520 

1559 

R3 


504 

520 

2079 

R3 


505 

520 

2599 

R3 


506 

520 

3119 

R3 


* 507 

590 

3709 

R4 

2106 

902-405 

200 

3909 

R4 


406 

200 

4109 

R4 


407 

200 

4309 

R4 


408 

200 

4509 

R4 


409 

200 

4709 

R4 


410 

520 

5229 

R4 

0210 

750-283 

300 

5529 

R4 


** 284 

300 

5829 


* BALL BEARING PRELOAD AND SUPPORT WAS LOST AT 3509 SEC. OF TEST TIME. 

** PUMP END BEARING TEST TIME WAS 5829 SEC. WITH 2320 SEC. (8 STARTS AND 9 STOPS) ON THE 
DAMPING SEAL BEARING. THE TURBINE END BEARING REACHED 8303 SEC. TEST TIME. 


HPOTP DAMPING SEAL BEARING ENGINE TESTS W/O BALL BEARINGS 
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HYDROSTATIC BEARING 



DAMPING BEARINGS 



STATOR FEED VERSION 


ROTOR FEED VERSION 
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HPOTP (LEFT) / HPFTP (RIGHT) TURBINE DAMPING BEARINGS 



SEALS k BEARINGS DYNAMIC A LEAKAGE DATA 

STIFFNESS DAMPING WHIRL LEAKAGE 


SEALS AND BEARINGS MN/m kNs/m FREQUENCY kg/s 

(klb/in) (lb s/in) RATIO (lb/s) 


HPOTP PBP DUPLEX BAIL BEARING 

70 

(400) 

3.5 

(20) 

n/a 

2.7 

(6)* 

HPOTP PBP DISCH. LABYRINTH SEAL 

3.5 

(20) 

4.0 

(23) 

0.50 

2.3 

(5) 

HPOTP PBP DISCH. DAMPING SEAL 

100 

(580) 

53 

(300) 

0.19 

2.7 

(6) 

HPOTP TURBINE HYDROSTATIC BEARING 

1100 

(6000) 

400 

(2300) 

0.50 

18 

(39) 

HPOTP TURBINE DAMPING BEARING 

1200 

(6800) 

670 

(3800) 

0.27 

9.1 

(20) 

HPFTP SMOOTH SEAL 

68 

(390) 

28 

(160) 

0.50 

0.6 

(1.4) 

HPFTP DAMPING SEAL 

26 

(150) 

44 

(250) 

0.17 

0.3 

(0.7) 

HPFTP TURBINE ROLLER BEARING 

610 

(3500) 

none 

n/a 

1.1 

(2.5)* 

HPFTP TURBINE HYDROSTATIC BEARING 

510 

(2900) 

89 

(510) 

0.50 

2.3 

(5.0) 

HPFTP TURBINE DAMPING BEARING 

530 

(3000) 

160 

(920) 

0.18 

1.1 

(2.5) 


* LEAKAGE THROUGH BALL AND ROLLER BEARINGS ARE COOLANT FLOWS. HPOTP IS THE HIGH PRESSURE 
OXYGEN TURBOPUMP. PBP IS THE PREBURNER PUMP OF THE HPOTP. HPFTP IS THE HIGH PRESSURE 
FUEL TURBOPUMP. HPOTP AND HPFTP DATA ARE FOR 30,000 RPM AND 36,000 RPM, RESPECTIVELY. 
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GENERIC BEARING COMPARISON 


ITEMS\BEARING TYPES 

BALL 

ROLLER 

HYDROSTATIC 

DAMPING 

MOVING PART COUNT 

33, duplex 

17 

1 

1 

DN, BORE mi x RPM 

< 2*10“ 

< 2*10® 

< 10 7 

< 10 7 

HOOP STRESS MPa (ksi) 

172 (25) 

310 (45) 

n/a 

n/a 

HERTZIAN ST. MPa (ksi) 

2412 (350) 

2412 (350) 

41 (6) 

41 (6) 

PRECISION un (mil) 

0.5 (0.02) 

0.5 (0.02) 

25 (1) 

25 (1) 

ROLLING STABILITY 

skidding 

skewing 

n/a 

n/a 

HIGH SIDELOAD 

wear 

skewing 

ok 

ok 

WHIRL ONSET SPEED 

limit 

limit 

limit 

none 

BEARING DAMPING 

negligible 

none 

high 

high 

BEARING TILTING 

< minute 

« minute 

> minute 

> minute 

DEADBAND/TIGHT FIT 

conflict 

conflict 

relaxed 

relaxed 

FAILURE MODE 

catastrophic 

catastrophic 

benign 

benign 

BEARING LIFE LIMIT 

wear 

skewing 

none * 

none * 


* NO WEAR AT FOIL SPEED. ONLY OCCASIONALLY LIGHT BRUSHING AT LOW SPEED IN START AND STOP. 


OUTLOOK 
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CROSS-FORCES FROM LABYRINTH SEALS. 
MECHANISMS AND UPSTREAM COUPLING 


By 

Knox Millsaps 

and Manuel Martinez-Sanchez 
MIT, Gas Turbine Laboratory 


OUTLINE 

• EXPERIMENTAL DATA AND THEORETICAL ANALYSIS OF 
ROTORDYNAMIC FORCES IN SINGLE-CAVITY LABYRINTHS. 


• PARAMETERS VARIED: INLET SWIRL. PRESSURE RATIO. 
WHIRL SPEED, SPIN SPEED. SEAL LENGTH. LAND SURFACE 
(SMOOTH VS. HONEYCOMB). 

• DATA WERE SEAL CAVITY REAL-TIME PRESSURE 
DISTRIBUTIONS 


MAJOR FINDINGS 


• PRESSURE FORCES ARE SUM OF IDEAL COMPONENT, 
PROPORTIONAL TO RELATIVE SWIRL, PLUS FRICTION- 
INDUCED COMPONENT, PROPORTIONAL TO SWIRL CHANGE 

• THE IDEAL COMPONENT VANISHES WHEN FLUID SWIRL 
EXACTLY FOLLOWS TRAVELLING WAVE OF GAP VARIATION 
DURING WHIRL. THIS COMPONENT IS SOLELY RESPONSIBLE 
FOR DAMPING 

• THE VISCOUS COMPONENT ALTERS CROSS-STIFFNESS, BUT 
NOT DAMPING 

• DAMPING DATA CAN BE EXTRACTED FROM STATIC TESTS 
WITH VARYING SWIRL 


• UPSTREAM NONUNIFORMITIES DUE TO SEAL ECCENTRICITY 
STRONGLY FEED BACK INTO CAVITY PRESSURE PATTERN 
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CROSS-FORCE MECHANISMS 


(1) I N VI SQP 



• WITH ECCENTRIC ROTOR, TANGENTIAL 
FLOW IN CAVITY SEES VARYING CROSS 
SECTION 

• IN REGIONS OF INCREASING AREA, 
CONTINUITY REQUIRES ADDITIONAL FLOW 
TO BE BROUGHT FROM UPSTREAM (OR A 
REDUCED DISCHARGE DOWNSTREAM) 

• THIS THEN MEANS LOW P IN THESE AREAS, 
WHICH INTEGRATES TO FORWARD- 
WHIRLING FORCE. 


- THIS ARGUMENT ASSUMES CONSTANT VELOCITY. BERNOULLI 
EFFECTS MODIFY (REDUCE) IT, BUT BASICS REMAIN 

- THE RELEVANT TANGENTIAL VELOCITY IS WITH RESPECT 
TO THE WHIRLING FRAME. I.E., V* - OR. THIS IS THE ORIGIN 
OF DAMPING (FORCE = 0 WHEN V* = OR) 


CRQ S S-F QRCE-MBQi A MSMS 


(2) VISCOUS 



BECAUSE OF WALL FRICTION, V* (IN CAVITY) 
CAN BE SMALLER (OR LARGER, AT HIGH co ) 
THAN Vj* (INLET) 

ASSUMING V* < Vi, INCOMING LEAKAGE FLUID 
ENERGIZES CAVITY FLOW > oj 

MORE FLOW ENTERS IN WIDER GAP REGION, 
SO P PEAKS 90° AFTER IT. THIS RESULTS IN 
FORWARD-WHIRLING FORCE 

NOTE ESSENTIAL VELOCITY HERE IS Vi - V*, 
REGARDLESS OF O. SO THIS FORCE 
COMPONENT IS INDEPENDENT OF WHIRL 
(NO DAMPING) 
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Figure 2*6: Test section schematic showing the geometry and flow variables when the swirl 
vanes feed into an upstream swirl cavity. Also shown is the leakage path to the large center 
cavity. 


SIMPLIFIED THEORETICAL PREDICTIONS 

• ASSUMING UNIFORM INLET AND EXIT CONDITONS. AND 
NEGLECTING A FEW MINOR TERMS, WE CALCULATE 



WHERE K » A ib. SENSITIVITY OF CARRY-OVER TO GAP) 

n, iS 2 

a = (SEAL DIVERGENCE, EQUAL TO I FOR US) 

V. = INLET TANGENTIAL VELOCITY 
V* = CAVITY TANGENTIAL VELOCITY 

V x = Vap/p = axial flow speed in gap 

Q. - WHIRL SPEED 

dF dF 

NOTE: C. .= = - R, THIS YIELDS C X X FROM STATIC DATA 



Table 7.1 The felt column ihowi the cross-stiffness correlation for all builds and for the 
static data of Benckert (281 The next two columns show the total cross force and frictional 
component evaluated at fl* = 300 (~ ! ) and design pressure ratio. The next column 
is the measured nondimensional direct damping coefficient. The final two columns give 
the damping calculated from the cross-stiffness coefficient and the average value that was 
directly measured. 



Figure 4-3: Cross section of nested spindle Spin/Whirl producing rotating rig and test 
section assembly. 
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Figure 4-8: Layout of the core facility showing the drive systems for both the spin and 
whirl. 


liar-umi bt/.i: j ;i*q nun.* r.w j * * ; 1 1 j *u ** 

(1.0) 

0.0370 

0'plsle 

1.27, 1.38. 1.48 
1.43. 1.83 

0 

±7.52. ±20.87 
±33.17, ±48.13 

(1.0) 

0.0370 

15* plate 

r 1.25. 1.39. 1.42 
1.45, 1.84 

0 

±7.52, ±20.67 
±33.17, ±46.15 

(10) 
0. 0370 

30* 

1.27. 1.38. 1.44 
1.44. 1.82 

0 

±1.52. ±20.87 
±33.17, ±48.15 

(1.0) 

0.0370 

40* 

^ 1.25. 1.40. 1.44 

1.68. 1.83 

0 

±7.52, ±20.87 
±33.17. ±46.15 

(3.8) 

01407 

0'pUle 

h 1.13. 1.21. 1.27 
1.38. 1.48. 1.65 

0 

±7.52, ±20.87 
±33.17, ±48.13 

(*•») 

0.1407 

0 # pkt« 

1.73. 1.79, 1.94 
2.07. 2.29. 2.40 

0 

±742. ±20.87 
±33.17, ±48.15 

(3.8) 

0.1470 

O'plaie 

1.48 

±44.87. ±66.52 

±7.52. ±20.87 
±33.17. ±46.15 

(3.8) 

0.1470 

13'pUu 

1.24. 1.40. 1.55 
1.64. 1.81 

0 

±7.52, ±20.87 
±33.17. ±48.15 

(3.8) 

0.1407 

13’piaU 

1.55 

±44.87. ±48.52 

~ ±7.52. ±20.67 

±33.17. ±48.15 

(3-8) 

0.1470 

30* 

1.23. 1.38. 1.33 
1.68. 1 76 

0 

±7.52. ±20.67 
±33.17, ±48.13 

(3-8) 

0.1470 

i 30* 

1.35 

1 ±44.87, ±64.52 

±7.52, ±20.87 
±33.17. ±48.13 

(38) 

0.1407 

60* 

1 11. l.»». 1.23 
1.38. 1.53. 1 66 

0 

±7.53. ±20.87 
±33.17. ±48.15 

(3.8) 

0.1407 

60* 

| 1 79. 1.88. 2.02 

2.21. 2.30. 2.38 

0 

±743, ±20.87 
±33.17. ±46.15 

(3-8) 

0.1407 

60* 

1.34 

±44.87, ±66.52 

±7.34, ±20.87 
±33.17. ±48.15 

(t-3) 

0.3704 

0*pki« 

1.20. 1.38, 1.58 
1 67. 1.88 

0 

“ ±7.53. ±20.87 
±33.17, ±48.15 

(T3) 

0.3704 

13*pUte 

1.21. 1 36. 1.33 
167. J 93 

0 

ii.52, ±20.87 
±33.17, ±48.15 

(107) 

0.3943 

O'plaie 

1.23. 1.39. 1.52 
1.83. 1.79 

0 

±7.52. ±20.87 
±33.17, ±48.15 

(10.7) 

0.3943 

60* 

1.25. 1.40. 1.46 
1.66. 1.83 

0 

±7.52. ±30.87 
±33.17. ±48.13 

(13.1) 

0.4832 

0*plaie 

1.27. 1.38. 1.46 
1.63. 1.63 

0 

” ±7.52. ±20.87 
±33.17. ±46.15 

(13.1) 

0.4832 

13*pltl«r 

1.25. 1.39. 1.42 
1.63. 1.84 

0 

±742, ±20.67 
±33.17, ±48.15 

(13-1) 

0.4831 

30* 

1.27. 1.38. 1.44 
1.66. 1.82 

0 

±7.52, ±20.87 
±33.17. ±48.13 

(131) 

0.4831 

80* 

1.25. 1.40. 1.46 
1 66. 1.85 

0 

±742. ±20.87 
±33.17. ±48.13 


Table 5.5: Test Matrix for build #3. 


261 



VARIABLE 

ABSOLUTE 

UNCERTAINTY 

RELATIVE 

UNCERTAINTY 

COMMENTS 1 

Kulite Sensitivity, 

±7.1 Pa 

±1.4% - ±2.2% of P 

linearity, hysteresis 
and calibration 

Repeatability P 

lor„ = 0.304N 

±8.5% 

95 % 


2 ar T - 0.148N 

±4.6% 

Confidence 

tnlet 

0.0 ± 0. 75^-0* 

undefined 

mass flow, 

swirl 

28.4 ± 2.2a- is# 

7.7% 

alignment, 

Vi 

40.7 ±3.6 a- 30* 

7.2% 

repeatability, 

Build#3 design 

69.6 ± 4.4a.5o* 

6.3% 

calibration 

Pi, P\ P* 

±344 Pa Max 

0.5 % each MAX 

Total 




accuracy 

Spin speed 
u 

±2.1 RPM 

±0.1 % MAX 

O 1275 RPM 

Whirl speed 

±4.5 RPM 

±0.9 % MAX 

O 520 RPM 

n 




Temperature 

1*P 

0.2% 

To calculate 

T 



density 


1.0± 0.2 mils 


Most probable 

* 

3.8± 0.2 mils 

1.5% 

error due to 

and 

7.3± 0.1 mils 

0.6% 

non-circularity 

h 

10.7± 0.1 mils 

0.2% 

of orbit 


13.1± 0.2 mils 

0.3% 

and casing. 


17.1± 0.3 rails 

0.3% 


lXX. 

±0.001 in. MAX 

0.5% MAX 

Machining 




tolerances 




Uncontrolled 

s: 

0.012±0.008 in. 

+67% -50% 

axial 




clearance 

Digitization 

k 

0.4% MAX 



Table 5.8: Table of uncertainties used in calculating the overall confidence interval for the 
force data. 


TYPICAL RESULTS FOR ZERO INLET SWIRL 


DATA SHOWN FOR CONF. #3 (SMOOTH, LONG SEAL) 

RESTORING DIRECT FORCE, INDEPENDENT OF SWIRL, SCALING 
WITH AP 

STABILIZING CROSS FORCE, ZERO AT ZERO WHIRL RATE, 
PROPORTIONAL TO WHIRL 


, SO DATA CAN BE COLLAPSED TO ONE LINE 

a p -Jap 
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Figure 7.13’. Experimentally obtained direct force, Fff, n. tbe whirl frequency, for five 
different inlet pressures, r, =1.27, 1.38, 1.47 1.69 and 1.79. There data are from build #3 
with 0* inlet s wlrl and ut =* 0. 



Figure 7.14: Experimentally obtained crow force. Ft, ti. the whirl frequency, for fire 
different inlet presrures, w , =1.27, 1.38, 1.47, 1.69 and 1.79. These data are from build #3 
with 0* swirl and u> = 0. 


I# * 10 TO THE CCNTIMCTCN M II W CM 
•t Krur.tL » cum co «*« » .»* 


46 1510 
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R ESULTS WITH INLET SWIRL 


• NONZERO CROSS FORCE AT ZERO WHIRL (Af^ * 0) . 

• EFFECT OF WHIRL SAME AS BEFORE. 

• SINCE F r ~ AP AND BOTH, TANGENTIAL VELOCITY AND CIR 

SCALE AS VA P, PLOT OF VS. -jEL STILL COLLAPSES DATA 

AP « JAP 

• SL0PE IS THE SAME AS WITH ZERO INLET SWIRL. 

SWIRL I NCREASES K xv . DOES NOT AFFECT Cw 



Figure 7.17: Experimentally obtained direct force, F Nt y*. the whirl frequency for fire 
different presiure ratio*. There data are from build #3 with 8.6* inlet iwfrl and w = 0. 



Rgure 7.18: Experimentally obtained cron force, Fr , ▼*• the whirl frequency, for flee 
different presture ratio*. These data are from build #3 with 8.6* inlet twirl and w = 0. 
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SHAF T BQT AJION-THE JR1CT10 N AL EFFECT 


• ROTATION IN THE INLET SWIRL DIRECTION REDUCES CROSS- 
FORCE Kxv (AND VICE-VERSA) 

• REDUCTION IS THE SAME AT ALL WHIRL SPEEDS (NO CHANGE 
IN DAMPING Cxx 

• EFFECT DUE TO VARIATIONS OF THE TANGENTIAL VELOCITY 
CHANGE V* - Vi BETWEEN INLET AND GLAND 
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Figure 7.19: Experimentally obtained direct force, F*, n. the whirl frequency for flee 
different spin rates, r, = 1.53. These data are from build #3 with 8.6* inlet swirl. 



Figure 7.20: Experimentally obtained cross force, Fy, vs. the whirl frequency, for fixe 
different spin rates, t, » 1.47. These data are from build #3 with 8.6* inlet swirl 



figure 7.21: Experimentally obtained direct force, Fs , ▼*. the whirl frequency for flee 
different spin rates, w, = 1.55. These data are from build #3 with 21.4* inlet swirl 



Figure 7.22: Experimentally obtained cross force, Fy, vs. the whirl frequency, for five 
different spin rates, *, = 1.55. These data are from build #3 with 21.4* inlet swirl. 
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EFFECTS OF HONEYCOMB LAND 


SOME REDUCTION (~ 25%) OF CROSS FORCE. 

PROBABLY DUE TO CARRYOVER DISRUPTION. THEORY GIVES 
17% REDUCTION WITH K = 0. 

MAJOR REDUCTION, EVEN REVERSAL, OF DIRECT FORCE 
(BECOMES NEGATIVE SPRING). SINCE F N DUE ALMOST 
ENTIRELY TO CARRYOVER DISRUPTION. 

DATA SHOWN FOR CONF. #4 (LONG, HONEYCOMB) FOR 
COMPARISON TO #3 (LONG, SMOOTH) 

EFFECTS FOR SHORT SEAL (#5 VS. #2) ARE WEAKER. 



Figure 7.23: Experimentally obtained direct force, Fff, va. the whirl frequency for fi*e 
different pressure ratios. These data are from build #4 with 0° inlet swirl and u> = 0. 
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NO N-DIMENSIONAL ROTORDYNAMIC coefficients 


FORCES FITTED AS 

y = - /(„ - c v a + n 1 

^ = K - C Q 

r 

COEFFICIENTS MADE NONDIMENSION AL BY 

K “ = 

' t R { P x - P) 

r - _ C, 5, 

c ' ~ / **. JpM-p.) 



7.29: The effect of pressure ratio on the nondimeasionnl direct stiffness coefficient, 
K£, for builds #2(oarrow rotoT-smooth land) , #3(wide rotor-smooth land), #4(wide rotor 
- honeycomb land) and #5( narrow rotor - honeycomb land). The inlet ssrirl for all cases is 
0* and «! =0.1407. 



^* ure 7 -30: The effect of pressure ratio on the non dimensional cross stiffness coefficient, 
for # 2 t #3i #4 and #5. All were taken with the 15* swirl orifice plate and 
<1=0.1407. 
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Figure 7.31: The effect of pressure ratio on the direct damping coefficient, for builds 
#2(narTOw rotor-smooth land) , #3(wide rotor-smooth land), #4(wide rotor • honeycomb 
land) and #5( narrow rotor - honeycomb land). The inlet swirl for all cases is 0* and 
<,=0.1407. 



Figure 7.32: The effect of pressure ratio on the nondimensional cross damping coefficient, 
C” for build #2, #3, #4 and #5. AU were taken with the 15* swirl orifice plate and 
<,=0.1407. 


EFEECTS OF UPSTREAM NONUNIFORMITIES 


• PREVIOUS STUDIES HAVE ASSUMED UNIFORM Pi ,Vj 

OUR RESULTS STRONGLY SUGGEST THIS IS NOT IN GENERAL A 
GOOD ASSUMPTION 

• SINCE EFFECT NOT SUSPECTED, NO PRESSURE DATA WERE 
TAKEN UPSTREAM, AND NO SPECIAL CARE WAS TAKEN TO 
MEASURE OR MAINTAIN A CONSTANT AXIAL GAP IN 
UPSTREAM FACE SEAL (NOMINAL - 10 mil) 

• LINEARIZED THEORY IGNORING UPSTREAM EFFECTS 
PREDICTS ALL TRENDS CORRECTLY, BUT YIELDS K XY , C XX 
VALUES 2-3 TIMES TOO LOW 

• THIS SAME DISCREPANCY HAD BEEN PREVIOUSLY NOTED BY 
SEVERAL RESEARCHERS WHEN TRYING TO CALCULATE 8P IN 
THE FIRST CAVITY OF MULTICAVITY SEALS 
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UESTREAM NDNUNIFORMITIES - CONTINUED 


MODEL WAS EXTENDED BY EXPLICITLY COUPLING THE INLET 
CAVITY CONDITIONS TO THOSE IN THE SEAL CAVITY, 
SOLVING FOR BOTH SIMULTANEOUSLY. 

• 

RESULTS SHOW VERY LARGE INCREASES OF Kxy AND Cxx 
(FACTORS OF 3-5) IF INLET CAVITY DOES NOT LEAK INTO HUB 
AREA. 


• ALLOWING FOR AN AXIAL GAP INTO A UNIFORM HUB VOLUME 
REDUCES THESE FACTORS, BUT FOR OUR NOMINAL 10 mil GAP, 
THE FACTORS ARE STILL 2-3. 

• WE ALSO REDUCED FIRST-CAVITY DATA OF BENCKERT- 
WACHTER. DISCREPANCY BETWEEN DATA AND THEORY CAN 
BE REMOVED IF AN AXIAL GAP OF - 10 mil IS CONSISTENTLY 
ASSUMED. 


NQIF/ THE SAME THEORY PREDICTS PRESSURE IN CAVITIES 
OTHER THAN 1ST TO WITHIN ± 15% 



Figure 3-41: Ratio of direct damping with upitream coupling to that with no coupling v«. 
iwirl cavity to seal gland area ratio. The maximum occurs a: 1.35 and the damping ratio 
comes to within 1% at a area ratio of 75. 



Figure 3-42: Predicted direct and cross force at Q = 0 from the fully coupled model vs. 
the relative axial clearance ratio, {*. The geometry is the same as for build #3. t # =1.4, 
<*. = 15*, f\ = 0.1407 and u> = 0. 
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CONFIG. | 

Fy(*uea j.) 

f t (s; = o) 

FtW = oo) 

Matched 1] 
1 

1 

10.2 1 ( N ) 

16.28(N) 

4.99(N) 

0.008" 

2 

8.281 N) 

16.15(N) 

4.25(N) 

0.011" | 

3 

ll.9l(N) 

L5.51(N) 

4.09(N) 

0.010" | 


Tible 7.2 Comparison of the data of Benckert and Wachter [28] to the coupled model 
predictions. All cases are for standard conditions. F, = 1.58 (bar), F. = 1 (bar) and 
a v = 28.4° The first column shows the experimental value. The second column gives 
the value predicted with full coupling, that is no leakage flow. The third column has the 
predictions for constant upstream boundary conditions. The last column gives the value of 
the axial space needed for the model to match the experimentally obtained value 


THEORY EXPERIMENT 



Figure 7.34: K“ n. a for the experimental data and theory. The axial gap, is used as 
a parameter. All experimental values fall between the theory with 0.004”(0.0001m) < f* < 
0.017(0. 0004m). The top thick line is for full coupling (ie. no leakage). The bottom one is 
for uniform inlet conditions (ie. no coupling). The calculations are for build #3 geometry. 
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DISCUSSION 


• THE SENSITIVITY OF CROSS-FORCES TO UPSTREAM 
CONDITIONS WIN. NECESSITATE ACCOUNTING FOR THESE IN 
I)KSl(iN WORK. 

• THIS WAS RELATIVELY SIMPLE FOR OUR GEOMETRY (NO 
BLADES), BUT WILL NOT BE SIMPLE FOR REALISTIC TURBINE 
SHROUD GEOMETRY. 

• THESE NON UNIFORMITIES ARE THEMSELVES DUE TO THE 
KOTOR ECCENTRICITY. SO NOT A-PRIORI KNOWN 

MORE RESEARCH IS NEEDED IN THIS AREA - 


CONC LUS IONS 


LABRYRINTH CROSS-FORCE MECHANISMS CLARIFIED 

• INVISC1I) (KINEMATIC) EFFECTS RESPONSIBLE FOR DAMPING. 
DAMPING INSENSITIVE TO SPIN. 

• VISCOUS MECHANISM AFFECTS CROSS-STIFFNESS, NOT 
DAMPING. 

• DAMPING CAN BE EXTRACTED FROM STATIC DATA. 

• INDUCED UPSTREAM NON-UNIFORMITIES CAN MORE THAN 
DOUBLE CROSS FORCES. THIS NEEDS TO BE ACCOUNTED FOR. 

• DIRECT FORCES MAINLY DUE TO CARRYOVER VARIATIONS 

• HONEYCOMB LAND REDUCES DIRECT FORCES STRONGLY, 
CROSS-FORCES WEAKLY. 
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STABILITY OF TWO-PHASE FACE SEALS 


From notes and discussion 


of talk by J. Yasuna, CMU 
by R.C. Hendricks 


A variable temperature model which considers squeeze film effects and thermal 
transients that was developed for two phase flows (1, 2) was extended to examine 
the axial stability of two-phase mechanical seals that tracks the axial response 
through position and velocity of the rotor from the perturbed state (3). 

Step increases in film thickness usually engendered a response that decayed as 
the rotor asympototical ly approached the initial steady state condition. 

Step decreases in film thickness that are large enough brought about monotone 
opening (seal failure) and other forms of failure through film collapse or "poping 
open" due to the onset of instabilities. 

In certain conditions the response to small step decreases leads to stable film 
thicknesses and rotor goes to steady state or some orbital state. 

Small amounts of coning tend to stabilize the configuration, even for large dis- 
turbances. The magnitude of the thermal response parameter is important to re- 
sponse time and overshoot. 

Model is limited through 1-D heat transfer and a conjugate solution would appear 
more appropriate and must incorporate the thermomechanical behavior of the seal 
rotor and stator interface, which includes seal "waviness." The distinct 
liquid-vapor boiling interface and laminar flow assumptions are also limitations, 
but the "smeared" two-phase interface with flow turbulence and transient behavior 
is quite complex but necessary for a better understanding of the dynamics. 


For more information 

1. Hughes, W.F., Basu, P., Beatty, P.A., Beeler, R.M., and Lau, S.: Dynamics of 
Face and Annular Seals With Two-Phase Flow, NASA CR 4256, 1989. 

2. Yasyna, J.A., and Hughes, W.F.: A Continuous Boiling Model for Face Seals, 
Trans. ASME J. Tribology April 1990, Vol . 112, pp 266, 274. 

3. Yasuna, J.A., and Hughes, W.F.: Squeeze Film Dynamics of Two-Phase Seals, 
Trans. ASME J. Tribology April 1992, Vol. 114, pp 236, 247. 
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HIGH PRESSURE COMPRESSOR DELIVERY BRUSH SEAL OF THE INTERNATIONAL 
AERO ENGINES (IAE) V2500-A1 GAS TURBINE ENGINE 


Peter A. Withers 
Rolls-Royce pic 
Filton, Bristol 


High pressure compressor delivery brush seal of the international aero engines (IAE) 
V2500-A1 gas turbine engine 

The V2500-A1 is the first production aero gas turbine to be certificated with brush seals. The engine h cs 
brush seals in three positions ofwhich the HPC seal discussed in this short presentation has the most 
arduous duty. 

Foil 1 shows the Rolls-Royce brush seal development "milestones" prior and subsequent to the entry 
into service of the V2500-A1 in 1989. Prior to the 1970's the potential of brush seals was known but 
limited by materials (early tests were conducted with Nylon bristles). Extensive R&D ws conducted 
during the 1 970's, continuing in parallel with demonstrator engine testing in the early 1 980's. By the 
mid 1980's sufficient technology and confidence in brush seals had been acquired to commit to 
application and development in the civil V2500-A1 engine with certification and entry into service in 
the late 1980's. 

Foil 2 shows the HPC brush seal in relation to the HPC delivery, the seal providing control of turbine 
zone ventilation and temperatures, lower sfc, and in conjunction with the labyrinth shown, a balance 
piston, for control of bearing thrust load. 

Foil 3 and 4 show constructional details of the brush seal and its associated journal, with Foil 5 
showing the nominal operating conditions which are quite arduous. Rolls-Royce experience shows 
that to obtain good seal life it is essential to finish the journals aluminium oxide coating by grinding 
and lapping to a surface finish of better than 0,5 micrometers CLA. 

The V2500-A1 engine designed with this HPC delivery, and two other, brush seals was successfully 
developed and certificated, Foil 6. During development the seals were found to quickly run in to a 
stabilised condition and provided the journal coating surface finish was good the seal would run 
indefinitely. Whilst craze/axial cracking of the aluminium oxide coating was observed this has not so 
far proved to be a problem. 

Foils 7 and 8 show the service experience to date. No unscheduled removals of brush seals have 
occurred and the intention is not to remove seals for inspection until 5,000 engine cycles, (current 
high time engine is 3,600 cycles). A few engines have been stripped for other reasons and visual 
inspection of the seals on these engines have been shown the seals to be in good condition. 

With the successful development and encouraging service experience of brush seals in the V2500- 
A1, Rolls-Royce will continue in the 1990's to exploit the advantages of brush seals through R&D, 
demonstrator engines and further engine applications. 
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R-R BRUSH SEAL DEVELOPMENT "MILESTONES" 


PRE 1970’s 

1972 ► 

198Q’s 

1985-1988 

1989 ► 

1990’s 


- POTENTIAL OF BRUSH SEALS KNOWN 
BUT LIMITED BY MATERIALS 

- R & D OF BRUSH SEALS 

- DEMONSTRATOR ENGINE TESTS 

- DEVELOPMENT FOR A CIVIL ENGINE 
APPLICATION 

- SERVICE APPLICATION IN CIVIL ENGINE 

- CONTINUING R 4 D AND APPLICATIONS 
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V2500-A1 HPC DELIVERY BRUSH SEAL DETAILS 
1. THE JOURNAL 


• JOURNAL DIAMETER - 315 mm 

• JOURNAL COATING - PLASMA SPRAYED 

ALUMINIUM OXIDE 

• FINISHED BY GRINDING AND LAPPING TO 
A SURFACE FINISH OF BETTER THAN 
0,5 MICRO METERS CLA. 

w 
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V2500-A1 HPC DEUVERY BRUSH SEAL DETAILS 
2. THE BRUSH SEAL 

• NUMBER OF SEALS IN SERIES - 2 

• BRISTLE FIT AT BUILD - RADIAL CLEARANCE 

• BRISTLE LAY ANGLE - 45° 

• BRISTLE LENGTH - 15 mm 

• BRISTLE WIRE DIAMETER - 0,071 mm NOMINAL 

• NUMBER OF BRISTLE ROWS PER SEAL - 1 1 

• BRISTLE DENSITY - 2250 TO 2650 PER 25 mm 

• BRISTLE WIRE MATERIAL - HAYNES 25 

BACKING RING CLEARANCE - 1,2 mm (RADIAL) 

© 1992 Roni-Royc* pie 
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V2500-A1 HPC DELIVERY BRUSH SEAL DETAILS 
3. OPERATING CONDITIONS 

• SURFACE VELOCITY - 250 m/s 

• DIFFERENTIAL PRESSURE - 430 kPa 

• TEMPERATURE - > 550 °C 
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V2500-A1 HPC DELIVERY BRUSH SEAL 
DEVELOPMENT/CERTinCATION 


• THE V2500-A1 ENGINE WAS SUCCESSFULLY DEVELOPED 
AND CERTIFICATED WITH BRUSH SEALS 

• OBSERVATIONS FROM THE DEVELOPMENT PROGRAMME: 

- SEALS QUICKLY RUN IN TO A STABILISED 
CONFIGURATION AND THEREAFTER , IF THE ROTOR 
CONDITION IS GOOD , THE SEAL WILL RUN INDEFINITELY 

- CRAZE/AXIAL CRACKING OF THE ALUMINIUM OXIDE 
JOURNAL COATING AT THESE EXACTING TEMPERATURES 
AND SPEEDS IS NORMALLY SEEN 

- BRISTLES WAXED OPEN FOR ENGINE ASSEMBLY 


© 1982 Rolli-ftoycc pic 
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V2500-A1 HPC DELIVERY BRUSH SEAL IN SERVICE 
1. RUNNING EXPERIENCE (TO END MAY 19921 


• 77 Al’s IN SERVICE 

• HOURS 

TOTAL 

HIGH TIME 7200 


354,300 

* 


CYCLES 

212,100 

3600* 


♦NOT SAME ENGINE - DIFFERENT USEAGE 
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V2500-A1 HPC DELIVERY BRUSH SEAL IN SERVICE 
2. SEAL CONDITION 


• FIRST SCHEDULED SERVICE REMOVALS WILL 
NOT TAKE PLACE BEFORE 5,000 CYCLES 

• NO UNSCHEDULED REMOVALS OF BRUSH SEALS 

• ON ENGINES STRIPPED FOR OTHER REASONS 
THE CONDITION OF THE BRUSH SEALS HAS 
BEEN GOOD 

00 

/ 
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"Modulus vs temp, for 1600°C DDS Fiber 
Bulk CVD (Gulden) and bulk alpha-SiC" 



Test temperature, °C 

Figure 15. Plot Young's modulus (GPa)of fiber made at 2912°F(1600°C) from DDS 
and hydrogen/versus temperature (°C) compared with literature data (ref. 24, 25). 
Good agreement is evident. Gulden work was on CVD material. 
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Comparison of 
creep rates at 
various loads for 
CVD and single 
crystal materials. 
From 0. Conran 


Thoria sc 


Alumina sc 


Yttrfa sc 




CVD b-SIC 
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